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SECTION  I 
INTRODUCTION 

Tills  investigation  concentrates  on  the  damage  produced  in  chemically 
vapor  deposited  (CVD)  zinc  sulfide  (ZnS)  by  water  drop  impacts  on  the  order 
of  1  mm  in  diameter  impinging  over  a  velocity  range  from  22.2  to  684  ms 
(730  to  2240  fps) .  A  combined  experimental  and  analytical  study  was  pur¬ 
sued  to  determine  the  relationships  between  water  drop  impact  damage  and 
the  material  properties  of  CVD  zinc  sulfide  in  order  to  evaluate  and  guide 
efforts  to  develop  erosion  resistant  window  materials  for  use  at  8-12  ym 
wavelengths.  The  general  problem  can  be  divided  into  three  major  cate¬ 
gories  and  each  of  these  categories  involves  a  number  of  experimental  and 
analytical  tasks.  The  suggested  categorization  outlined  below  is  essen¬ 
tially  our  overview  of  the  general  problem  and  places  some  perspective  on 
the  relative  balance  between  experiment  and  analysis. 

1,  Fracture  Initiation 

Identification  of  the  material  properties  associated  with  fracture 
initiation  due  to  a  single  water  drop  Impact  (fracture  threshold) . 

Identification  of  the  material  properties  associated  with  fracture 
initiation  due  to  multiple  water  drop  impacts  below  the  single  water  drop 
fracture  threshold. 

Identification  of  the  material  properties  controlling  fracture 
initiation  and  growth  at  water  drop  impact  velocities  well  above  the  single 
water  drop  fracture  threshold. 

Evaluation  of  the  change  in  the  fracture  response  due  to  modification 
of  a  controlling  material  property. 

Development  of  a  model  for  fracture  initiation. 

2.  Crack  Growth 

Identification  of  the  material  properties  controlling  crack  growth. 

Evaluation  of  the  change  in  crack  growth  and  crack  morphology  due 
'  to  modification  of  a  controlling  material  property. 
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Evaluation  of  the  relation  between  crack  growth  and  morphology 
and  infrared  transmittance. 

Development  of  a  model  for  crack  growth. 

3.  Erosion  Pit  Nucleation  and  Growth 


Identification  of  the  erosion  mechanisms  responsible  for  mass 
removal  (erosion  pit  nucleation). 

Identification  of  the  material  properties  controlling  the  erosion 
process  (erosion  pit  growth) . 

Evaluation  of  the  change  in  erosion  rate  and  erosion  pit  morphology 
due  to  modification  of  a  controlling  material  property. 

Evaluation  of  the  relation  between  erosion  pit  growth  and  morphology 
and  infrared  transmittance. 

Development  of  a  model  for  the  general  erosion  process. 

Development  of  a  model  for  the  prediction  of  the  degradation  of 
infrared  transmittance. 

Not  all  of  the  above  tasks  are  relevant  for  the  specified  impact 
conditions  nor  have  all  that  are  relevant  been  addressed  in  this  study. 

For  example,  the  first  two  tasks  on  fracture  initiation  are  not  directly 
applicable,  since  the  fracture  threshold  for  zinc  sulfide  is  already  ex¬ 
ceeded  at  the  minimum  velocity  of  interest  (222  ms  ^) .  Erosion  pit  growth 
was  not  investigated  extensively,  since  the  objective  of  the  program  can 
be  satisfied  without  pursuing  this  complex  subject.  Analytical  and  experi¬ 
mental  investigations  of  the  events  leading  to  material  removal  (the  in¬ 
cubation  period  for  erosion  pit  nucleation)  are  sufficient  for  providing 
the  desired  material  guidance  for  the  development  of  more  erosion  resis¬ 
tant  windows.  The  primary  emphasis  of  the  program  was  on  crack  growth 
and  to  a  lesser  extent  on  the  events  leading  to  erosion  pit  nucleation. 

The  relevant  velocity  range  for  this  investigation  determined  from  an 
initial  overview  of  water  drop  impact  damage  as  a  function  of  impact  veloc¬ 
ity  reported  in  Section  III, 2  is  from  342  to  512  ms  ^  (1120  to  1680  fps). 
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».  Water  drop  impact  experiments  were  carried  out  on  CVD  zinc  sulfide 

1 

specimens  in  the  ETI  Liquid  Drop  Impact  Facility  and  are  described  In  Sec¬ 
tion  III.l.  This  facility  provides  controlled  and  well  characterized 
impact  conditions.  The  water  drops  fall  from  the  tip  of  a  hypodermic 
needle  and  have  diameters  in  the  range  from  1,5  to  2  mm.  Although  many 
of  the  qualitative  aspects  of  the  rain  erosion  response  of  brittle  mate¬ 
rials  are  known,  the  water  drop  impact  experiments  on  CVD  zinc  sulfide  are 
useful  for  they  provide  specimens,  for  detailed  materials  characterization 
in  order  to  identify  the  dominant  surface  features  and  micros tructural 
properties  which  influence  fracture  Initiation  and  growth. 

Mathematical  analyses  of  water  drop  impacts  on  an  elastic  target  with 
and  without  a  surface  crack  have  been  formulated.  The  general  analysis  is 
based  on  Green's  theorem  formulated  for  the  convolution  of  the  displacement 
vector  and  the  stress  tensor.  It  was  initially  thought  that  the  extension 
suggested  by  Hackworth  and  Kocher  (1978)  of  the  TURBAN  computer  program 
could  be  used  to  generalize  Blowers*  analytic  solution  for  the  stresses  in 
an  elastic  half-space  due  to  a  water  drop  collision  (Blowers,  1969).  It 
was  found  upon  closer  examination  that  this  extension  is  quite  restrictive 
and  in  fact  in  the  form  suggested  is  not  valid  for  the  pressure  distribu¬ 
tion  used  by  Blowers.  A  simple  generalization  of  this  work  is  easily  found, 
however,  which  is  valid  for  Blowers*  pressure  distribution  but  the  resulting 
distribution  is  still  so  restrictive  that  this  extension  of  Blowers*  work 
cannot  be  usefully  incorporated  into  the  present  program  as  originally  pro¬ 
posed,  This  led  to  an  independent  effort  to  provide  an  analytic  generali¬ 
zation  of  Blowers*  work  in  conjunction  with  consideration  of  crack  growth 
in  a  transient  stress  field. 

An  analytic  formulation  of  the  crack  growth  problem  in  materials 
subject  to  water  drop  impact  is  developed  which  does  not  involve  explicit 
numerical  integration  of  the  equations  of  motion.  The  latter  approach, 
while  straightforward,  has  a  number  of  disadvantages  not  the  least  of  which 
is  the  large  expense  for  the  required  computer  time.  Other  disadvantages 
include  the  lack  of  transparency  of  the  solution  with  regard  to  parametric 
variation  and,  at  a  more  fundamental  level,  the  problem  of  disentangling 
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subtle  physical  effects  from  numerical  artifacts  of  the  finite  grid  size 
used.  In  explicit  numerical  studies,  parametric  studies  are  done  by  gen¬ 
erating  several  solutions  for  a  range  of  values  of  the  parameters,  a  pro¬ 
cedure  which  is  usually  Impractical  for  more  than  one  or  two  parameters 
because  of  the  great  expense  involved,  for  example,  in  finite  difference 
schemes.  Nevertheless,  it  is  often  the  case  that  explicit  numerical 
schemes  are  the  only  alternative  (particularly  in  the  case  of  nonlinear 
problems  for  which  analytic  techniques  are  generally  lacking)  and,  irres¬ 
pective  of  the  expense,  they  do  have  the  distinct  advantage  that  solutions 
can  always  be  found  provided  the  singularities  of  the  solution  are 
understood. 

For  linear  or  quasilinear  problems  the  possibility  exists,  at  least 
in  principal,  for  obtaining  analytic  solutions  to  the  problems  in  question. 
Analytic  is  used  in  the  sense  that  it  includes  the  use  of  numerical  methods 
for  evaluating  integrals  and  summing  series,  but  excludes  the  use  of  large 
grids  of  values  of  the  independent  and  dependent  variables.  It  is  the  use 
of  the  latter  which  accounts  for  the  large  cost  of  the  explicit  numerical 
schemes.  The  two  approaches  are  complimentary  in  many  respects.  Parametric 
dependence  is  made  explicit  by  analytic  solutions  and  approximations  made  to 
achieve  a  solution  are  usually  better  understood  so  that  the  accuracy  of  the 
solutions  obtained  can  be  better  assessed.  For  interesting  problems,  on  the 
other  hand,  analytic  solutions  are  almost  always  more  difficult  to  obtain. 

If  found,  however,  they  provide  more  information  at  less  expense  than  do  the 
explicit  numerical  methods.  For  this  latter  reason,  it  is  always  advanta¬ 
geous  to  expend  some  effort  on  any  problem  to  determine  if  at  least  an  ap¬ 
proximate  analytic  solution  can  be  found  before  undertaking  the  development 
of  an  explicit  numerical  solution.  Even  if  the  analytical  approach  does 
not  yield  a  solution,  one  usually  benefits  from  the  effort  by  achieving  a 
better  understanding  of  the  problem,  in  particular  the  singularities,  which 
helps  to  understand  the  numerical  results  generated  by  the  explicit  schemes. 
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SECTION  II 
MATERIALS 


All  of  the  formulations  of  chemically  vapor  deposited  (CVD)  zinc 
sulfide  (ZnS)  investigated  were  produced  by  the  Research  Division  of 
Raytheon  Corporation.  A  baseline  material  was  supplied  by  AFML  which 
was  fabricated  as  a  FLIR  window.  The  piece  of  material  received  was 
designated  Plate  A,  Run  No.  MT-207,  Section  B-2,  and  was  62x104x17  mm 
thick.  Subsequently  Section  G-2  was  also  received  from  Run  No.  MT-207. 

This  material  was  selected  as  the  baseline  material  since  it  has  been 
characterized  by  the  University  of  Dayton  Research  Institute  (Graves, 
et  al.,  1977;  Wimmer  and  Graves,  1977,  1978;  Wimmer,  1979).  This  sample 
will  be  referred  to  as  material  A  in  this  report.  According  to  Wimmer  and 
Graves  (1977)  it  was  deposited  at  640°C  for  250  hours. 

Additional  formulations  of  zinc  sulfide  have  been  received  from  the 
Research  Division  of  Raytheon  Corporation  for  the  purpose  of  investigating 
the  influence  of  different  microstructures  on  water  drop  impact  damage. 

Two  of  the  additional  foraulations  deposited  at  690  and  728‘’C  (Material  B 
and  C,  respectively)  were  investigated.  No  physical  or  mechanical  property 
data  have  been  received  for  these  materials. 

1.  INITIAL  CHARACTERIZATION 

Materials  A  and  B  are  amber  in  color,  while  Material  C  has  a  milky 
yellow  appearance.  Visual  inspection  of  these  materials  reveals  bands  of 
color  change  and  differing  inclusion  concentrations  parallel  to  the  deposi¬ 
tion  plane. 

Table  1  is  a  summary  of  the  characterization  data  generated  by 
Wimmer  and  Graves  (1977)  and  Wimmer  (1979)  for  material  A.  Microhardness 
measurements  were  made  on  the  polished  surfaces  of  the  specimens  used  in 
the  water  drop  impact  tests.  The  results  are  presented  in  Table  2. 
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TABLE  1 


I 

► 
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Properties  of  CVD  ZnS,  MT-207, 

:fc 

Plate  A  (Material  A) . 


Lattice  Parameter 

Density 

Grain  Size 

4-Point  Bend  Strength 
** 

^IC 

Weibull  Parameters 
(surface  flaw  model) 


Knoop  Hardness  (100  g  load) 
Young's  Modulus 
Poisson's  Ratio 

Calculated  Dilatational  Wave 
Velocity 

Calculated  Distortional  Wave 
Velocity 

Thermal  Expansion 

Thermal  Diffusivity 
Thermal  Conductivity 
Specific  Heat 


5.409  A 
4090  kg/m^ 

8-10  ym  in  growth  plane;  up  to 
60  ym  (columnar)  in  growth  direction 

100.8  MPa  (14.6  ksi) 

0.67  MN/m^^^  (609  psi  ) 

g 

a  =  29.0  X  10  (dimensions  in  meters, 
stress  in  pascals) 

3 

0  =  14.3  X  10  (dimensions  in  inches, 

stress  in  psi) 

m  =  6.07 

2 

239  kg/mm  (growth  plane) 

74.1  GPa  (10.7  x  10^  psi) 

0.30  (in  plane  of  plate) 

0.35  (in  growth  direction) 

4,960  ms  ^ 

2,600  ms  ^ 

6.33  X  10"^/°K  (-67‘’C  to  25'’C) 

7.42  X  10"^/°K  (25°C  to  200®C) 

0.09  cm^/sec  (25°C) 

0.146  w/cm/°K  (25°C) 

0.40  j/g/“K 


Taken  from  Wimmer  (1979) . 

Crack  propagation  parallel  to  columnar  grains. 
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TABLE  2.  Knoop  Microhardness  Measurements  (100  g  Load) 


MATERIAL 

KNOOP  HARDNESS 

STANDARD  DEVIATION 

A 

2 

238  kg/mm 

5  kg/mm^ 

B 

244  kg/mm^ 

2 

8  kg /mm 

C 

2 

194  kg/mm 

2 

8  kg/mm 

TABLE  3.  Characteristic  Grain  Dimensions 

for  CVD  ZnS  lormulations  Investigated. 


AVERAGE  LINEAR  INTERCEPT  (ym) 
MATERIAL  A  MATERIAL  B  MATERIAL  C 


Growth  plane  3-4  3-5  5-10 

Cross  section  perpendicular  5-30  5-30  10-100 

to  growth  plane  in  the  direction 
perpendicular  to  growth  plane 


Aspect  ratio  range 


1.25-10 


1-10 


1-20 


General  grain  structure  and  growth  patterns  were  investigated  by 
chemically  etching  polished  surfaces.  The  etchant  was  equal  parts  con¬ 
centrated  HCl  and  H2O  at  45“C.  Etching  time  was  1  minute.  Surfaces 
both  parallel  and  perpendicular  to  the  growth  plane  were  etched.  Re¬ 
presentative  grain  structure  micrographs  are  shown  in  Figure  1  for 
Material  A.  Materials  B  and  C  exhibited  the  same  qualitative  features 
shown  in  Figure  1. 

Due  to  the  highly  anisotropic  nature  and  local  variations  of  the 
grain  size  and  structure,  an  average  grain  size  dimension  is  somewhat 
meaningless.  Average  intercept  distances  were  measured  both  in  the 
growth  plane  and  perpendicular  to  it  for  quantitative  comparison  of 
the  grain  sizes  in  these  materials.  The  grains  appear  equiaxed  in 
the  growth  plane  but  columnar  in  the  perpendicular  plane  as  can  be 
seen  in  Figure  1.  The  intercept  distances  are  given  in  Table  3. 

No  difference  in  grain  structure  was  noted  between  Section  B-2  and  G-2 
for  the  baseline  material. 

In  addition  to  having  columnar  grain  morphologies,  the  three  ZnS 
formulations  exhibited  local  grain  structures  which  will  be  referred  to 
as  haystacks  which  describes  their  general  appearance.  Examples  of  the 
haystack  grain  growth  patterns  are  shown  in  Figure  2.  A  haystack  con¬ 
sists  of  a  nucleation  site  and  localized  distortion  of  the  characteristic 
columnar  grain  growth  pattern.  At  the  nucleation  site,  the  major  axes  of 
the  grains  radiate  out  into  the  half  space  above  the  growth  plane.  As  more 
grains  are  formed,  their  major  axes  tend  to  approach  the  growth  direction 
leaving  the  haystack  structure.  Partially  developed  haystack  structures 
appear  to  have  been  noticed  by  Graves,  et  al.  (1977),  in  thin  zinc  sulfide 
layers  deposited  on  zinc  selenide  substrates.  They  were  referred  to  as 
rosettes  by  these  investigators.  The  haystack  structures  nucleate  at  clinker 
and  rod  inclusions.  These  inclusions,  which  have  a  maximum  size  on  the  order 
of  100  ym,  tend  to  exist  in  growth  bands  parallel  to  the  growth  plane.  The 
inclusions  are  of  a  larger  average  size  (c.  80  ym)  in  material  C  and  approx¬ 
imately  ten  times  more  numerous  than  in  material  A  or  B  (average  size 
c.  50  ym) . 
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Growth  Plane  Section  at  the 
Nucleation  Site  of  a  Haystack 
Structure  in  Material  A. 


Transmitted  Polarized  Light 
View  of  Material  A  Showing 
the  Haystack  Structure  and 
Deposition  Bands  of  Differ¬ 
ing  Grain  Texture, 


Figure  2.  Examples  of  Microstructural  Features  in  CVD' ZnS. 
(Arrows  indicate  growth  direction.) 
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In  addition  to  bands  of  inclusions  all  materials  exhibited  bands  of 
altered  grain  structure.  These  bands  can  be  seen  in  Figure  2d  photo¬ 
graphed  in  transmitted  polarized  illumination. 

2.  SAMPLE  PREPARATION  PROCEDURES 

Specimens  for  water  drop  impact  were  fabricated  in  the  form  of  right 
circular  cylinders  17  mm  diameter  by  12  mm  thick.  The  specimens  of  mate¬ 
rial  C  were  only  4  mm  thick  due  to  limited  plate  thickness. 

The  circular  impact  faces  were  oriented  parallel  to  the  growth  plane, 
and  then  mechanically  polished  to  optical  quality  by  A.M.F,  Optics,  Inc,, 
Woburn,  MA.  Additional  specimens  with  their  impact  faces  oriented  perpen¬ 
dicular  to  the  growth  plane  were  prepared  at  ETI. 

In  the  initial  phases  of  the  program  all  of  the  impact  specimens 
were  polished  at  ETI.  Tliese  surfaces  were  found  to  contain  excessive 
polishing  induced  flaws  which  interacted  strongly  with  the  impact  loading. 
The  quality  of  the  surface  finish  achieved  at  ETI  with  further  polishing 
experience  has  approached  the  level  supplied  by  A.M.F,  Optics,  Inc. 

3.  POST-TEST  EXAMINATION  PROCEDURES 

After  impacting  a  water  drop  in  the  ETI  Liquid  Drop  Impact  Facility, 
the  specimen  is  carefully  removed  from  the  sabot  (the  experimental  procedure 
is  described  in  Section  III.l).  Recovery  debris  are  removed  with  methanol 
alcohol  in  a  low  energy  ultrasonic  cleaner  followed  by  a  methanol  and  then 
ether  rinse.  The  major  fractures  and  outflow  damage  are  photographed  in 
both  reflected  and  transmitted  illumination  for  future  reference.  The 
specimen  is  typically  etched  (50%  HCl  @45°C)  to  reveal  the  fine  surface 
fractures  present.  The  enlianced  fracture  pattern  is  then  photographed. 

The  subsurface  morphologies  of  the  fractures  are  revealed  by  sectioning 
the  specimen  on  a  plane  that  is  perpendicular  to  the  impact  plane  and  passes 
through  the  center  of  the  impact.  The  sectioning  is  accomplished  by  first 
encapsulating  the  impact  face  with  Buehler  mounting  epoxy.  The  sectioning 
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process  involves  rough  material  removal  by  wet  sanding  with  progressively 
finer  grits  of  SiC  paper,  terminating  with  5  ym  A.I2O2  in  water  on  nylon. 
This  is  followed  by  removal  of  approximately  100  ym  of  material  with  0.3  ym 
AI2O2  in  water  on  nylon  and  finally  polishing  the  desired  section  with 
.05  ym  silk  wet  with  50%  HCl.  This  sectioning  method  exposes  the 

internal  fractures  essentially  undisturbed  for  optical  microscope  examina¬ 
tions  using  various  illuminations. 

The  finer  details  of  the  fracture/grain  structure  interactions  on 
the  impacted  face  of  the  specimen  are  observed  from  etched  surfaces  using 
scanning  electron  microscopy.  Surface  replication  procedures  were  used 
to  examine  the  details  of  the  damage  within  surface  grains  using  trans¬ 
mission  electron  microscopy. 


SECTION  III 
WATER  DROP  IMPACTS 
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This  program  concentrates  on  understanding  water  drop  impact  damage 
in  CVD  zinc  sulfide.  The  surface  loading  conditions  associated  with  a 
water  drop  collision  are  still  somewhat  ill-defined,  although  advances  are 
being  made  in  analyzing  water  drop  collisions  at  high  subsonic  to  low  super¬ 
sonic  impact  velocities  (Adler,  1979a).  Since  the  water  drop  impact  condi¬ 
tions  are  not  well  defined  quantitatively,  the  nature  of  the  loading  condi¬ 
tion  as  well  as  the  response  of  the  target  becomes  an  integral  part  of  the 
general  problem.  The  ETI  Liquid  Drop  Impact  Facility  currently  provides 
ideal  (spherical)  water  drop  impacts  at  velocities  from  222  to  684  ms  ^ 

(730  to  2240  fps)  along  with  well-documented  data  on  individual  drops  just 
prior  to  collision  with  the  target. 

The  operation  of  the  ETI  Liquid  Drop  Impact  Facility  is  described  in 
Section  III.l,  The  system  components  and  instrumentation  used  in  the  facility 
were  for  the  most  part  drawn  from  what  was  readily  available.  Quite  accept¬ 
able  results  are  presently  obtained  in  this  facility  as  shown  in  Section  III, 2. 
The  potential  exists,  however,  for  enhanced  and  extended  capability  through 
additional  modifications  and  component  upgrading.  These  modifications  are 
being  made  as  justified.  The  mechanics  of  water  drop  collisions  required 
for  better  understanding  of  the  internal  fracture  observations  reported  in 
Section  III. 2  is  described  in  Section  III. 3. 

1.  DESCRIPTION  OF  LIQUID  DROP  IMPACT  FACILITY  AND  PROCEDURES 

A  schematic  of  the  ETI  Liquid  Drop  Impact  Facilic>  is  shown  in 
Figure  3,  The  specimen  is  mounted  at  the  front  end  of  a  saoot  which  is 
propelled  down  the  range  by  a  small  charge  of  gunpowiler.  Water  drops  are 
formed  at  the  tip  of  a  hypodermic  needle  mounted  above  thu  test  section 
and  fall  in  the  path  of  the  moving  specimen.  A  slotted  rail  system  is  used 
to  guide  the  sabot  through  the  test  section  and  into  the  rc  -.overy  tube.  As 
a  single  drop  falls  from  the  tip  of  the  hypodermic  needle  it  interrupts 
dual  fiber  optics  links  above  the  test  chamber  thereby  initiating  the 
firing  sequence.  The  system  is  calibrated  so  the  sabot  reaches  the  drop 
when  it  has  fallen  to  a  point  within  the  limits  of  the  exposed  surface  of 
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Figure  3.  Liquid  Drop  Impact  Facility 


the  specimen.  After  impacting  the  water  drop  the  sabot  enters  the  recovery 
tube  where  a  graded  distribution  of  material  controls  the  deceleration  of 
the  sabot  and  protects  the  face  of  the  specimen  from  extraneous  damage. 

The  barrel  of  the  powder  gun  has  an  inside  diameter  of  30  mm  and  is 
capable  of  accelerating  the  sabot  to  a  maximum  speed  of  1370  ms  ^  (c.  4500 
fps).  Circular  specimens  up  to  22  mm  (0.875  in.)  in  diameter  can  be  accom¬ 
modated  with  the  existing  sabot  configuration.  Non-circular  specimens  with 
lateral  dimensions  less  than  22  mm  and  a  wide  range  of  thicknesses  can  also 
be  used.  The  zinc  sulfide  specimens  used  in  this  program  are  17  mm  in  diam¬ 
eter  and  range  from  6  to  10  mm  thick.  These  dimensions  are  adequate  to 
avoid  edge  effects  and  to  eliminate  any  stress  wave  interactions  with  the 
water  drop  impact  damage  due  to  stress  wave  reflections  from  the  back  face 
of  the  specimen.  Specimens  of  zinc  sulfide  can  be  recovered  after  water 
drop  collisions  over  the  velocity  range  of  interest,  222  to  684  ms  ^  (730 
to  2240  fps  corresponding  to  an  aircraft  traveling  at  500  mph  through  mach  2), 
without  cracking  the  specimens  or  producing  additional  damage  which  obscures 
that  due  to  the  drop  impact.  As  shown  in  Figure  3  the  total  length  of  the 
recovery  system  is  only  8  m  (26  ft.).  The  extent  of  the  water  drop  impact 
damage  on  zinc  sulfide  over  this  velocity  range  is  described  in  Section  III. 2. 
Multiple  water  drop  impact  damage  can  be  obtained  by  sequential  firings  of 
the  same  specimen  down  the  range.  Only  a  few  interacting  drop  impacts  are 
required  at  342  ms  (1120  fps)  to  significantly  degrade  the  infrared  trans¬ 
mittance  of  zinc  sulfide.  The  ETI  facility  (described  more  fully  by  Adler 
and  James,  1979)  is  ideally  suited  for  documenting  each  drop  impact  in  con¬ 
junction  with  its  effect  on  the  general  degradation  process. 

The  tests  were  conducted  in  a  helium  atmosphere  at  an  ambient  pressure 
of  300  pm  Hg  in  order  to  minimize  distortion  of  the  drop  caused  by  the  bow 
shock  prior  to  impact  with  the  specimen's  surface.  This  effect  is  distinct 
from  that  of  blow-by  which  results  from  the  expansion  of  the  gases  associated 
with  the  gunpowder  burn.  Escape  of  these  gases  generated  by  the  detonation 
of  the  gunpowder  has  been  eliminated  by  means  of  a  self-sealing  sabot  configu¬ 
ration.  A  considerable  developmental  effort  has  been  completed  to  maintain 
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the  integrity  of  the  water  drop  and  to  minimize  its  aerodynamic  distortion 
^prior  to  impact.  labile  shock  wave  and  boundary  layer  effects  are  an  inte¬ 
gral  part  of  the  erosion  problem  for  supersonic  aircraft,  the  drop  geometry 
prior  to  impact  in  actual  flight  environments  is  an  aerodynamic  problem. 

The  intent  here  is  to  understand  the  material's  response  to  water  drop  colli- 
sions  which  requires  a  reproducible  and  relatively  simple  drop  geometry. 

In  order  to  operate  at  pressures  well  below  the  vapor  pressure  of 
water,  a  solution  of  80  volume  percent  of  distilled  water  and  20  volume 
percent  ethylene  glycol  (C2Hg02)  is  used.  From  a  pragmatic  viewpoint  this 
solution  behaves  essentially  the  same  as  pure  water  during  the  collision 
process. 

The  water-glycol  solution  is  deaerated.  The  temperature  of  the 
solution  is  typically  slightly  below  room  temperature:  approximately 
15°C.  Chilling  the  solution  has  the  advantage  of  lowering  the  vapor 
pressure,  but  this  advantage  has  to  be  balanced  against  an  increase  in 
viscosity.  The  deaerated  water-glycol  solution  is  placed  in  a  reservoir 
above  the  test  chamber  where  it  is  drawn  into  a  vertically  mounted  syringe. 
Manually  applied  pressure  to  the  syringe  initiates  drop  formation  at  the 
tip  of  the  hypodermic  needle.  If  the  drops  are  forming  satisfactorily  at 
the  test  pressure  when  viewed  through  an  optical  microscope  with  a  magni¬ 
fication  of  about  seven  times,  a  single  drop  is  photographed  as  it  falls 
between  the  guide  rails  at  the  elevation  of  the  specimen.  The  integrity 
and  general  appearance  of  the  drop  is  checked  from  this  photograph.  If 
it  passes  this  visual  inspection,  the  pressure  in  the  system  is  checked 
to  confirm  that  it  matches  the  desired  test  pressure,  the  pre-set  firing 
sequence  is  unlocked,  and  detonation  of  the  gunpowder  takes  place  as  a 
subsequent  drop  interrupts  the  dual  fiber  optics  links  passing  above  the 
test  chamber  corresponding  to  station  in  Figure  3 .  When  the  sabot 
is  25  mm  from  the  impact  point,  it  passes  through  a  laser  beam  at  station 
d)  which  sends  a  pulse  to  a  delay  generator  which  subsequently  triggers 
a  flash  unit  to  photograph  the  drop  prior  to  Impact,  The  delay  circuit 
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allows  the  drop  to  be  photographed  as  close  to  the  impact  site  as  the 

jitter  in  the  system  will  permit.  In  all  cases  this  photograph  will 

record  the  water  drop  when  the  specimen  is  less  than  25  mm  away.  Stated 

differently,  the  sabot  will  be  a  maximum  of  113  ys  away  from  impacting 

-1  -1 

the  water  drop  when  traveling  at  222  ms  and  36  ys  away  for  a  700  ms 
impact  condition.  The  velocity  of  the  sabot  is  measured  by  the  length 
of  time  required  for  the  sabot  to  pass  through  the  laser  at  station  @ 
in  Figure  3.  Knowing  the  length  of  the  sabot  and  this  time,  the  veloc¬ 
ity  can  be  calculated.  A  second  determination  of  the  impact  velocity  is 
made  by  measuring  the  time  required  for  the  sabot  to  pass  through  laser 
beams  (^'  and  @  in  Figure  3.  The  distance  between  these  stations 
divided  by  the  recorded  time  increment  yields  a  second  value  for  the 
impact  velocity. 

During  the  initial  phases  of  this  investigation  the  water  drops 
were  tending  toward  ellipsoidal  shapes  prior  to  colliding  with  the  speci¬ 
men  due  to  the  interaction  with  the  gas  remaining  in  the  gun  barrel.  The 
pumping  capability  was  subsequently  improved.  At  a  system  pressure  of 
100  to  200  ym  Hg  in  a  helium  atmosphere,  the  drop  distortions  appear  to 
be  negligible  prior  to  impacting  the  specimen.  A  series  of  tests  were 
carried  out  for  confirmed  spherical  water  drop  impacts  as  shown  in 
Figure  4.  The  significance  of  ellipsoidal  drop  collisions  compared 
with  spherical  drop  impacts  will  be  considered  along  with  the  nature  of 
the  loading  conditions  for  spherical  water  drop  collisions  on  ZnS. 

2.  OVERVIEW  OF  WATER  DROP  IMPACT  DAMAGE 

Water  drop  impacts  have  been  obtained  on  the  three  ZnS  materials 

described  in  Section  TI.l,  and  additional  impacts  were  obtained  on  specimens 

from  an  unidentified  ZnS  FLIR  window  over  the  velocity  range  from  222  to 
-1 

684  ms  .  The  complete  range  of  test  conditions  and  the  general  features 
of  the  resulting  impact  fractures  will  be  described. 
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Figure  4.  Typical  Appearance  of  Water  Drop  Prior  to  Impact 
for  Impact  Velocities  from  200  to  700  ms“l  at  a 
Pressure  of  0.1  torr  or  less. 


a.  Test  Conditions 

Tables  4  and  5  list  the  impact  conditions  for  each  test  carried 
out  within  this  program. 

The  impacts  obtained  from  the  initial  sequence  (Table  4)  were  for 
surfaces  prepared  at  ETI.  At  the  time  the  initial  test  sequence  was  under¬ 
taken  the  drops  were  deforming  to  ellipsoidal  shapes  prior  to  impact,  how¬ 
ever  the  water  drop  impacts  listed  in  Table  5  are  for  confirmed  spherical 
drop  collisions.  An  ellipsoidal  drop  will  be  considered  equivalent  to  a 
spherical  drop  when  they  have  the  same  mass.  Let  bj^,b2  be  the  semi-major 
and  semi-minor  axes  of  the  ellipsoid,  then  the  equivalent  spherical  drop 
radius  r  is  determined  from  the  relation 


The  equivalent  drop  diameter  is  calculated  using  the  dimensions  measured 
from  the  photographs  of  the  drop  prior  to  impact  as  shown  in  Figure  4, 

The  ratio  of  the  major  to  minor  axis  of  the  elliptical  cross  section  is 
listed  in  Table  4  to  provide  a  measure  of  the  distortion  which  occurred. 

The  influence  ellipsoidal  drop  impacts  have  on  the  magnitude  of  the  ob¬ 
served  damage  will  be  examined  in  Section  III (a).  The  missing  entries  under 
the  heading  b^/b2  are  for  shots  in  which  the  drop  was  already  obscured 
by  the  sabot  when  the  photograph  of  the  drop  prior  to  impact  was  taken; 
for  these  cases  which  arise  due  to  the  inherent  jitter  in  the  system  the 
diameter  of  the  pretest  drop  is  used  for  the  equivalent  spherical  drop 
diameter. 


Improvements  in  the  ETI  Liquid  Drop  Impact  Facility  were  made  during 
the  early  stages  of  this  program  which  eliminated  drop  distortions  prior 
to  impact:  spherical  drop  impacts  are  now  obtained  at  all  impact  velocities. 
A  second  series  of  water  drop  impact  tests  were  carried  out  for  this  condi¬ 
tion  for  optically  polished  surfaces  as  listed  in  Table  5.  The  nature  of 
the  impact  damage  during  the  initial  sequence  did  not  show  any  change  in 
its  general  form  over  the  velocity  range  from  237  to  640  ms  ,  so  the  second 
test  series  concentrated  on  two  impact  velocities:  342  and  512  ms  ^ 
corresponding  to  mach  1  and  imich  1.5. 
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Specimens  sectioned  for  microscopic  examination 


To  our  knowledge  single  water  drop  impacts  on  zinc  sulfide,  or  zinc 
selenide,  above  342  ms  ^  (1120  fps)  have  not  been  investigated,  Hackworth 
and  Kocher  (1978)  have  reported  single  drop  impacts  on  ZnS  and  ZnSe  at 
velocities  to  340  ms  Their  experimental  facility  does  not  allow  veri¬ 
fication  of  the  drop  sizes  and  shapes.  A  number  of  multiple  drop  impact 
sequences  have  been  reported  in  their  work  (Hackworth  and  Kocher,  1977, 
1978).  A  series  of  interacting  spherical  drop  impacts  at  512  ms  ^  was 
initiated  in  order  to  assess  the  erosion  pit  nucleation  process  at  this 
higher  velocity  (shot  numbers  1049,  1051,  1052),  However  shortly  after 
obtaining  these  three  shots,  the  quality  of  the  impacts  was  declining  due 
to  gun  barrel  erosion  which  did  not  allow  the  gases  generated  by  the  gun¬ 
powder  burn  to  be  completely  sealed  by  the  sabot.  The  gun  barrel  is  being 
replaced,  so  the  multiple  impact  sequence  has  been  terminated  until  the 
previous  high  level  of  perfection  in  the  water  drop  collisions  is  once 
again  achieved. 

Qualitatively  the  nature  of  single  drop  impact  damage  in  polycrystal¬ 
line  ZnS  over  the  velocity  range  investigated  (237  to  640  ms  ^)  shows  a 
characteristic  fracture  pattern  whose  extent  depends  strongly  on  the  im¬ 
pact  velocity.  The  characteristic  fracture  pattern  as  viewed  in  the  direc¬ 
tion  of  impact  with  transmitted  illumination  is  apparent  in  Figure  5. 

Both  the  velocity  dependence  of  damage  and  the  characteristic  cross  sec¬ 
tional  fractures  due  to  single  drop  impacts  can  be  seen  in  Figure  6. 

Figure  6  shows  the  dramatic  difference  in  the  extent  of  the  damage  for 
a  water  drop  impact  at  299  ms  ^  compared  with  a  drop  impact  at  517  ms 
The  same  magnification  is  used  for  both  micrographs. 

When  viewed  with  reflected  light  single  drop  impacts  characteris¬ 
tically  appear  as  in  Figure  7(a).  Figure  7(b)  is  the  same  impact  (as  in 
Figure  7(a))  but  photographed  with  transmitted  illumination.  In  both 
transmitted  and  reflected  illumination,  it  is  apparent  that  the  fracture 
pattern  in  Figure  7  has  been  altered  by  interaction  with  linear  surface 
scratches.  The  surface  of  the  specimen  in  Figure  7  (shot  no.  836)  is 
from  the  initial  sequence  which  had  a  high  occurrence  of  surface  polishing 
scratches.  The  impact  damage  for  the  multiple  impact  sequence  is  shown  in 
transmitted  Illumination  in  Figure  8. 


Figure  5.  Shot  No.  985  (Transmitted  Light  Illumination) 
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Photographed  with  Reflected  Light  Showing  the 
Lateral  Outflow  Chipping  Damage, 


7.  Impact  Face  of  Shot  No.  836  Impacted  at  640  ms 
Note  Strong  Interaction  of  Fracture  Pattern  with 
Polishing  Marks. 


b)  Photographed  with  Transmitted  Light  Showing 
Radial  Extent  of  Subsurface  Fractures. 


Figure  7 


.  (continued).  Impact  Face  of  Shot  No.  836  Impacted 
at  640  ms“l.  Note  Strong  Interaction  of  Fracture 
Pattern  with  Polishing  Marks. 
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Crack  Morphologies  and  Dimensions 


In  order  to  obtain  data  on  the  nature  of  the  fractures  produced  in 
the  interior  of  a  specimen,  cross  sectional  views  on  a  plane  normal  to  the 
impact  surface  passing  near  and  through  the  center  of  the  damage  zone  are 
employed.  Three  distinct  fracture  zones  can  be  identified  by  the  crack 
trajectories  observed  in  these  cross  sections.  The  three  primary  fracture 
zones  are  believed  to  correspond  to  the  nature  of  the  surface  loading  con¬ 
ditions  and  the  resulting  transient  stress  states  within  the  target.  The 
fracture  zones  will  be  labeled  numerically  beginning  at  the  point  of  impact 
and  moving  outward  in  a  radial  direction. 

A  typical  impact  fracture  pattern  at  Mach  1  is  shown  in  Figure  9* 

A  portion  of  the  impact  face  and  the  principal  cross  section  is  shown. 
Readily  visible  in  the  center  of  this  figure  as  a  region  free  of  fracture 
is  zone  1,  alternatively  called  the  central  undamaged  zone.  The  central 
undamaged  zone  exists  since  the  initial  stresses  within  the  material  are 
predominantly  compressive  (Adler,  1977), 

Radially  outward  from  the  undamaged  zone  is  an  annulus  of  densely 
spaced  (both  circumferentially  as  well  as  radially)  fractures  which  in 
cross  section  are  monotonically  increasing  in  depth  with  Increasing  radial 
distance.  The  fractures  adjacent  to  zone  1  are  numerous  and  shallow  ranging 
from  3  to  50  ym  in  depth,  which  is  on  the  order  of  the  columnar  grain  dimen¬ 
sions  (Table  3 ) .  Progressing  radially  outward  the  fractures  become  less 
densely  distributed,  deeper,  and  approach  the  surface  at  a  shallow  angle 
(that  is,  more  nearly  tangent  to  the  surface).  The  cross  sectional  views 
of  the  midplane  fractures  show  a  very  distinct  transition  to  fractures 
which  are  nearly  normal  to  the  impact  surface  at  the  face  of  the  specimen 
and  then  penetrate  the  interior  of  the  specimen  in  an  initially  concave 
outward  direction;  however,  if  the  penetration  is  deep  enough,  the  frac¬ 
ture  path  will  exhibit  a  reversed  curvature.  Fracture  paths  with  double 
curvature  (concave  outward  near  the  surface  and  concave  in  at  greater 
depths)  are  observed  for  all  the  velocities  and  appear  to  form  a  transi¬ 
tion  from  the  fractures  which  are  essentially  tangent  to  the  surface  to 
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Figure  9 . 


Impact  Face  and  Cross  Section  of  Shot  No.  1016 
(331  ms”l)  Showing  Characteristic  Fracture  Zones 


those  which  tend  to  be  more  normal  to  the  surface  at  larger  radial 
distances.  This  fairly  well-defined  transition  in  the  orientation  of 
the  fracture  trajectories  will  define  the  outer  boundary  of  zone  2. 

The  general  character  of  the  fractures  within  zone  2  is  shown  in 
Figure  10.  The  deepest  fractures  occur  predominantly  at  the  outer 
radius  of  zone  2,  or  occasionally  at  the  extreme  inner  radius  of  the 
next  zone,  zone  3.  Zone  3  encompasses  all  of  the  fractures  outside 
of  zone  2. 

The  zone  3  fractures  are  characterized  by  a  nearly  vertical  approach 
to  the  impact  surface  for  the  first  10  to  50  ym  of  depth  depending  on  the 
impact  velocity.  The  innermost  fractures  which  are  the  deepest  in  zone  3 
exhibit  a  double  curvature  as  mentioned  earlier  which  can  be  seen  fairly 
clearly  in  Figure  6  .  The  depth  of  fracture  in  zone  3  decreases  initially 
quite  rapidly  and  the  fractures  lose  their  double  curvature  becoming  rela¬ 
tively  shallow  and  both  radially  and  circumferentially  sparse. 

In  addition  to  bulk  fracturing  of  the  target  material  water  drop 
impacts  also  produce  lateral  outflow  chipping  at  the  surface  crack  open¬ 
ings  (Adler,  1979a).  The  chipping  results  from  the  impingement  of  the 
lateral  outflow  jets  upon  the  outer  upraised  edges  of  the  bulk  fractures 
which  are  a  characteristic  feature  of  water  drop  impact  fractures.  The 
outflow  chipping  damage  is  most  readily  seen  with  reflected  illumination 
as  in  Figure  7(a). 

Characteristic  dimensions  of  the  radial  extent  and  maximum  fracture 
depths  for  the  fracture  zones  were  measured  on  the  impacts  which  were  sec¬ 
tioned  as  noted  in  Table  4  and  5  .  The  dimensions  of  the  various  fea¬ 
tures  of  the  impact  damage  described  above  are  summarized  in  Tables  6 
and  7  .  Due  to  the  level  of  variability  inherent  in  the  fracture  forma¬ 
tion  process,  a  considerably  greater  number  of  measurements  would  have  to 
be  made  to  provide  statistical  data  which  could  then  be  used  to  interpret 
the  trends  in  the  observed  fracture  zones.  The  measurements  in  Tables  6 
and  7  are  simply  intended  to  illustrate  the  velocity  dependence  of  the 
characteristic  features  of  the  damage  region  which  have  been  identified. 
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Some  of  these  dimensions  are  better  defined  than  others,  such  as  the 
radius  of  the  central  undamaged  zone,  zone  1,  the  onset  of  lateral  outflow 
damage,  and  the  maximum  radius  of  zone  2.  No  attempt  has  been  made  to 
develop  a  normalization  procedure  to  account  for  differences  in  the  drop 
diameters  or  shapes. 

3.  MECHANICS  OF  WATER  DROP  IMPACTS 

The  experimentally  observed  damage  was  produced  by  both  spherical 
and  ellipsoidal  shaped  drops.  The  difference  between  the  surface  loading 
characteristics  of  ellipsoidal  drops  will  be  investigated  for  comparison 
with  spherical  drop  impact  conditions.  In  addition,  the  general  features 
of  the  transient  stress  conditions  within  the  ZnS  specimens  will  be  con¬ 
sidered  in  relation  to  the  Impact  conditions  represented  by  the  experimen¬ 
tal  program  in  order  to  interpret  and  correlate  the  experimental  results. 
The  need  for  further  refinements  in  the  analytical  approach  is  pointed  out 
and  an  initial  effort  to  broaden  the  scope  of  the  stress  wave  analysis  is 
undertaken  in  Section  V. 


a.  Ellipsoidal  Drop  Collisions 

The  analysis  of  ellipsoidal  drop 
formulated  by  Adler  and  James  (1979). 
since  it  is  relevant  to  the  assessment 
observed  in  the  ZnS  specimens. 


impacts  on  rigid  surfaces  has  been 
This  analysis  will  be  extended  here, 
of  the  water  drop  impact  damage 


The  geometi'y  associated  with  a  perfectly  compressible  ellipsoid 
colliding  with  a  rigid  plane  is  indicated  in  Figure  11  .  The  equation  for 
the  elliptical  section  is 


2.2  22 
r  fiin  rt  r  rn<5  n 


(3.2) 


In  terms  of  the  relevant  impact  parameters. 


(3.3) 


(3.4) 


where 


c  is  the  radius  of  the  circular  contact  area 

V  is  the  impact  velocity 
o 

t  is  the  time  elapsed  from  initial  contact 

<l>  is  the  contact  angle  between  the  surface  of  the  ellipsoid 
and  the  rigid  plane. 


The  time-dependent  radius  of  the  contact  area  is  found  from  Eq.(3.3), 


,(t)  -  Uj  -  (v_,t)' 


(3.5) 


Differentiating  Eq.(3.5)  yields  the  rate  of  expansion  of  the  contact  radius, 
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According  to  Figure  10,  Eq.(3,6)  can  also  be  written 
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Using  Eq.(3.4),  Eq.(3.7)  becomes 
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tan  ((i(t) 


(3.8) 
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The  expression  in  Eq.(3.8)  is  identical  to  that  for  a  spherical  body. 
This  observation  implies  that  the  velocity  of  the  boundary  of  the  ex¬ 
panding  contact  area  will  be  the  same  for  these  two  cases  when  the  con¬ 
tact  angles  are  equal.  The  corresponding  contact  radii  at  which  this 
occurs  can  be  evaluated.  Using  the  equation  for  an  elliptical  section 
in  Eq.(3.2),  it  is  easy  to  show 


(3.9) 


Denoting  the  radius  of  the  contact  zone  for  a  spherical  drop  collision  by 
a,  then 


a  =  r  sind) 
o 


(3.10) 


where  r^  is  the  radius  of  the  spherical  drop. 

As  stated  in  Section  III, 2, a  an  ellipsoidal  drop  will  be  considered 
equivalent  to  a  spherical  drop  when  they  have  the  same  mass.  The  equiva¬ 
lent  drop  radius  is  determined  from  Eq,(3,l),  The  major  and  minor  semi¬ 
axes  of  the  equivalent  ellipse  can  be  determined  by  reorganizing  the  param¬ 
eters  in  Eq.(3.1)  such  that 


(3.11) 


When  the  radius  and  distortion  of  the  spherical  drop  are  prescribed, 
and  bj^  can  be  evaluated.  As  an  example,  consider  when  r^  =  1  mm  and 
b^/b2  =  1.5  ,  then  b2  =  0.765  mm  and  b^^  =  l,5b2  =  1.145  mm.  The 
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dimensions  of  the  equivalent  ellipsoid  can  be  found  in  general  from  the 
graphical  representation  of  the  relations  for  the  equivalent  ellipsoid 
provided  in  Figure  12  .  The  cross  sectional  shapes  of  ellipsoids  covering 
the  general  range  of  possible  interest  are  also  shown  in  Figure  12,  so 
some  idea  of  the  extent  of  the  drop  deformation  can  be  realized. 

The  previous  analysis  assumes  the  drop  is  perfectly  compressible 
which  is  a  reasonable  assumption  for  the  initial  stages  of  the  impact 
event.  However  incompressibility  effects  will  become  significant  in  con¬ 
junction  with  the  onset  of  lateral  outflow  (Adler,  1979a),  and  it  would 
seem  that  differences  in  the  nature  of  the  lateral  outflow  velocities 
would  occur  between  a  spherical  and  ellipsoidal  drop  impact.  It  is  our 
j:onjecture  that  whatever  differences  may  arise  would  be  relatively  small. 
Furthermore,  it  is  Adler’s  finding  (Adler,  1979a)  that  an  adequate  cri¬ 
terion  for  the  onset  of  lateral  outflow  for  a  spherical  drop  does  not 
exist,  although  a  number  of  analyses  can  be  found  in  the  literature.  A 
fairly  wide  dicotomy  prevails  between  the  elementary  analyses,  numerical 
calculations,  and  experimental  measurements  for  disc-shaped  drops.  To 
date,  there  are  no  direct  measurements  of  the  critical  lateral  outflow 
parameters  for  spherical  drop  impacts.  In  view  of  this  situation  the 
widely  accepted  lateral  outflow  condition  stipulated  by  Bowden  and  Field 
(1964)  is  generally  employed.  Their  criterion  is  simply 

V 

sin<^c  "  (3-^2) 
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where  is  the  compressional  wave  speed  for  water  which  is  approximately 
1.5  mm/ys.  For  the  range  of  Impact  velocities  under  investigation  the  wave 
speed  should  be  replaced  in  Eq.(3.12)  by  the  shock  wave  speed  for  water, 
U^.  According  to  Heymann  (1968), 

U  =  C  (1  +  2M  ) 
w  w  o 


(3.13) 


The  relation  in  Eq,(3,13)  Is  valid  for  Impacts  on 


where  M  =  v  /C 
o  o  w 

a  rigid  surface  when  M  <  1,2  .  Eq.(3.12)  now  takes  the  form 
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(3.14) 


This  is  a  lateral  outflow  criterion  for  spherical  drops,  however  the 
same  condition  pertains  to  elliptical  drops,  since  this  is  a  condition 
on  velocity  components  at  the  periphery  of  the  contact  zone^  namely 


*  2  2  2 
c  +  V  =  U 
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(3.15) 


It  was  already  shown  that  the  expansion  velocities  of  the  contact  zone 
were  identical  for  a  sphere  and  an  ellipsoid  for  the  same  value  of  the 
contact  angle  <p  (Eq.(3.8)).  The  critical  contact  radii  for  which  this 
condition  pertains  are  determined  from  Eq,(3.9)  and  (3.10)  upon  substi¬ 
tution  of  Eq.(3.14),  then 
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(3.16) 


Alternatively  the  condition  that 


c  =  U 
c  w 


(3.17) 


has  also  been  used  as  a  criterion  for  estimating  the  onset  of  lateral 
outflow.  However  the  relation  was  conjectured  to  be  related 

to  the  location  of  the  peak  pressures  in  the  pressure  distribution  found 
from  numerical  studies  of  water  drop  collisions  on  rigid  surfaces  (Rosen¬ 
blatt,  et  al.,  1979).  The  occurrence  of  the  peak  pressures  at  the  criti¬ 
cal  radii  computed  from  Eq.(3.17)  appears  to  conform  to  the  experimental 
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results  obtained  for  spherical  water  drop  impacts  on  polymethylmethacrylate 
specimens  (Adler,  1979b).  The  critical  contact  radii  for  ellipsoidal  water 
drops  are  found  from 
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(3.18) 


The  numerical  evaluation  of  this  expression  is  quite  close  to  the  value 
obtained  from  Eq.(3.16):  the  values  from  Eq.(3.16)  being  slightly  higher. 
The  critical  radii  are  plotted  in  Figure  13  as  a  function  of  the  im¬ 
pact  velocity  v^  . 

The  influence  of  water  drop  distortions  prior  to  impact  on  the 
dimensions  of  the  contact  zone  between  the  drop  and  the  target  (which  is 
assumed  to  be  directly  related  to  the  central  undamaged  zone)  can  now  be 
evaluated  explicitly.  Consider  a  2  mm  spherical  water  drop  impinging  at 
500  ms  The  radius  of  the  contact  zone  is  0.2  mm  as  determined  from 
Eq.(3.18).  However,  based  on  the  plots  in  Figures  12  and  13,  if  the 
initially  spherical  water  drop  deforms  into  an  ellipsoid  prior  to  colli¬ 
sion  with  the  target,  the  critical  contact  radius  increases  to 
0.328  mm  (b^/b2=1.5),  0.466  mm  (t)^/b2  =  2.0),  and  0.734  mm  (b^/b^  =  3.0). 

It  is  seen  that  the  size  of  the  contact  zone  can  be  increased  significantly 
for  even  relatively  minor  drop  distortions  assuming  the  stated  criteria 
for  lateral  outflow  and  the  location  of  the  peak  pressure  are  applicable 
in  general. 


The  time  after  impact  for  which  the  critical  contact  radii  are 
reached  corresponding  to  Eq.(3,16)  and  (3.18)  can  also  be  specified. 
The  critical  time 


t^  corresponding  to  Eq.(3.16)  is 
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(3.19) 


and  the  critical  time  corresponding  to  Eq.(3.18)  is 


The  duration  of  the  loading  time  is  also  greater  for  ellipsoidal  drops 
compared  with  spherical  drops. 

The  aforementioned  effects  should  have  relatively  little  influence 
on  the  pressure  distribution  applied  to  the  surface  of  the  target,  since 
lateral  outflow  for  the  ellipsoidal  drops  considered  will  occur  before 
the  shock  wave  in  the  liquid  reaches  the  back  surface  of  the  drop.  In 
contrast,  note  the  different  conditions  which  would  prevail  for  a  pancake¬ 
shaped  drop.  On  the  basis  of  the  previous  analysis  the  major  difference 
is  that  the  direct  pressure  would  be  applied  over  a  larger  contact  zone 
for  a  longer  duration  for  the  ellipsoidal  drop  as  compared  with  a  spheri¬ 
cal  drop  of  the  same  mass,  but  the  magnitude  of  the  applied  pressure  should 
be  roughly  the  same  in  both  cases.  The  pressure  distribution  on  the  inter¬ 
face  would  have  to  be  determined  from  a  detailed  analysis  of  the  flow  con¬ 
ditions  at  the  expanding  interface  as  already  carried  out  by  Huang  (1971), 
Hwang  (1975),  and  Rosenblatt,' et  al.  (1979). 

The  estimates  for  the  dimensions  of  the  contact  zone  are  derived 
for  a  rigid  surface,  however  the  difference  for  water  impacts  on  a  zinc 
sulfide  surface  is  not  too  significant;  the  difference  in  the  shock  wave 
velocity  for  the  two  cases  is  less  than  6%.  The  magnitude  of  the  uniform 
pressures  which  develop  over  the  contact  zone  according  to  a  one-dimensional 
wave  analysis  were  determined  and  are  shown  in  Figure  14  as  a  function  of 
the  impact  velocity. 
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Aerodynamic  effects  should  have  a  perturbing  influence  on  water  drop 
geometry  prior  to  impact,  in  all  water  drop  impact  facilities  employing  a 
specimen  moving  at  high  subsonic  to  supersonic  velocities.  If  severe  per¬ 
turbations  are  not  produced  (small  eccentricities),  it  is  extremely  diffi¬ 
cult  to  identify  the  ellipsoidal  drop  collisions  from  post-test  examinations 
of  the  impact  specimens.  Therefore,  the  photographic  records  of  the  actual 
collision  in  the  ETI  Liquid  Drop  Impact  Facility  are  very  useful  because 
they  afford  the  opportunity  to  quantify  and  correct  these  subtleties  in  the 
water  drop  impact  conditions  and  the  extent  of  the  resulting  damage  which 
have  gone  unnoticed  in  most  liquid  drop  erosion  studies. 

b.  Transient  Conditions  Within  the  Target  Materials 

The  generation  of  stress  waves  due  to  water'  drop  collisions  on  zinc 
sulfide  specimens  will  be  considered  in  order  to  better  realize  the  manner 
by  which  the  material  fails  and  for  the  development  of  correlations  with 
the  loading  cycle.  Blowers*  analysis  (1969)  has  been  employed  by  Adler 
(1977),  Hackworth  and  Kocher  (1977,  1978),  and  Hackworth,  et  al.  (1979) 
to  describe  the  state  of  stress  in  an  elastic  half-space  due  to  a  water 
drop  collision.  Rosenblatt,  et  al.  (1977,  1979)  have  employed  finite 
difference  procedures  to  address  this  problem.  Representative  results  from 
the  analytical  approach  have  been  described  (Adler,  1977,  1979a)  and  only  a 
few  observations  will  be  made  which  will  be  used  to  show  the  relation  be¬ 
tween  the  fracture  patterns  for  the  zinc  sulfide  specimens  and  the  form  of 
the  stresses  imposed.  A  new  and  more  general  analytical  approach  for  the 
determination  of  the  transient  stress  states  due  to  an  arbitrary,  local 
pressure  distribution  has  been  formulated  in  Section  V.  The  analysis  of 
the  water  drop/surface  crack  interactions  is  derived  in  Section  VI. 

The  characteristic  fracture  pattern  for  water  drop  impacts  on  ZnS 
over  a  velocity  range  from  222  to  684  ms  ^  is  shown  in  Figure  6.  The 
fractures  are  arrayed  along  circumferential  paths  around  a  central  undam¬ 
aged  zone,  zone  1  defined  in  Section  III.2.b.  The  fracture  concentration 
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decreases  with  radial  distance.  The  subsurface  character  of  the  fractures 
(depth  and  general  trajectory;  also  varies  with  radial  distance  as  shown 
on  a  central  cross  section  in  Figure  9  .  The  characteristic  fractures 
are  formed  at  or  near  the  surface;  all  of  the  damage  produced  over  the 
specified  velocity  range  lies  within  an  annulus  around  the  central  undam¬ 
aged  zone  and  extends  to  relatively  shallow  depths  within  the  specimen. 

The  analyses  which  follows  will  accordingly  focus  on  the  conditions  at 
and  near  the  surface  of  the  specimen  during  the  water  drop  loading  cycle. 

When  a  spherical  water  drop  strikes  the  surface  the  boundary  of  the 

contact  zone  will  expand  in  accordance  with  Eq,(3,5)  when  b^=b2=r^  .  The 

initial  velocities  of  the  boundary  of  the  contact  zone,  Eq,(3.6),  are  quite 

-1/2 

high,  but  decrease  essentially  as  t  ,  As  described  by  Blowers  (1969) 
a  condition  is  reached  whereby  the  dilatational  wave  for  the  target  material 
moves  ahead  of  the  contact  zone  followed  by  the  shear  wave  at  a  later  time. 
The  time  and  radial  distance  at  which  these  separations  occur  can  be  found 
by  equating  the  velocity  of  the  contact  zone,  a(t)  ,  in  turn  to  the  dila¬ 
tational  wave  speed  for  zinc  su3fide,  Cj^=4,96  mm/ys  ,  and  to  the  shear 
wave  speed,  €2=2. 60  mm/ys  .  The  derivation  is  analogous  to  that  carried 
out  for  a  similar  condition  in  Eq.(3.17),  The  final  results  are  plotted 
in  Figure  14,  which  clearly  shows  the  relation  between  the  propagation  of 
stress  waves  and  the  hypothesized  attainment  of  the  peak  pressure  and  the 
ensuing  precipitous  pressure  rele£,se  phase  of  the  collision  process.  It 
is  interesting  to  note  that  for  impacts  above  550  ms  ^  this  analysis  indi¬ 
cates  that  the  shear  wave  will  move  ahead  of  the  contact  zone  boundary  after 
the  peak  pressure  has  been  reached.  Also  due  to  the  increasing  shock  wave 
velocity  in  water  as  the  impact  velocity  increases,  the  time  required  for 
the  attainment  of  the  peak  pressure  remains  fairly  constant  over  an  extended 
range  of  impact  velocities. 

The  relations  provided  in  Figure  15  provide  some  perspective  on  the 
size  of  the  loaded  region  during  the  pressure  buildup  phase  and  the  time 
duration  required  for  the  indicated  sequence  of  events  to  occur.  The  tem¬ 
poral  development  of  the  radial  stresses  for  a  2  mm  water  drop  impacting  a 
ZnS  target  at  200  ms  ^  is  provided  in  Figure  16  .  The  stresses  are  computed 
at  a  depth  of  5  ym  below  the  surface  from  Blowers*  solution  to  the  water 
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drop  impact  problem  which  is  for  a  uniformly  loaded  circular  region  whose 
boundary  is  expanding  as  the  square  root  of  the  time  after  impact.  Neg¬ 
lecting  the  variation  in  the  pressure  distribution  during  the  loading  phase 
the  computed  stress  levels  up  to  35  ns  should  be  reasonable.  The  subsequent 
decline  in  pressure  predicted  from  Figure  15  to  occur  at  27  ns  is  not  re¬ 
flected  in  these  calculations.  The  succeeding  time  increments  are  provided 
to  show  the  evolution  of  the  radial  stress  for  this  special  case  and  to  il¬ 
lustrate  the  general  character  of  the  radial  stress  as  it  propagates  through 
fracture  zones  2  and  3.  The  valley  between  the  two  peaks  in  the  tensile 
stress  distribution  seen  when  the  contact  time  exceeds  35  ns  corresponds  to 
the  location  of  the  Rayleigh  surface  wave:  the  shear  wave  is  propagating 
just  slightly  ahead  of  the  Rayleigh  wave.  The  maximum  extent  of  the  radial 
stress  values  corresponds  to  the  location  of  the  dilatational  wave  front. 

The  magnitude  of  the  stress  is  normalized  with  respect  to  the  magnitude  of 
the  impact  pressure.  For  a  water  drop  impact  on  zinc  sulfide.  Figure  14 
indicates  the  magnitude  of  the  impact  pressure  is  340  MPa.  For  example  the 
magnitude  of  the  peak  radial  tensile  stress  when  t  =  45  ns  is  279  MPa. 

A  significant  state  of  triaxial  tensile  stress  develops  in  conjunction 
with  the  Rayleigh  wave.  At  a  depth  of  5  ym  and  when  t  =  45  ns,  a^^=179  MPa, 
*^00"89  MPa,  and  a^^=169  MPa.  The  spatial  extent  of  this  stress  condition  is 
on  the  order  of  a  few  micrometers  and  is  quite  localized  in  the  vicinity  of 
the  surface  of  the  target. 

Referring  to  Figure  16  ,  it  is  seen  that  significant  radial  tensile 
stresses  are  present  at  a  radial  distance  of  0.1  mm  when  t  =  25  ns  .  The 
potential  exists  for  the  initiation  of  circumferential  fractures.  Accord¬ 
ing  to  the  plot  in  Figure  15,  the  peak  pressure  on  the  surface  will  be 

achieved  when  a  =0.105  mm.  At  this  time  the  distortional  wave  would  have 
c 

advanced  beyond  the  contact  zone,  but  the  occurrence  of  this  condition  is 
dependent  on  the  magnitude  of  the  stress  wave  velocities  in  the  target 
material. 


49 


? 


Similar  plots  can  be  given  for  the  higher  impact  velocities.  The 
magnitude  of  the  impact  pressure  would  be  increased  as  shown  in  Figure  lA, 
so  the  magnitude  of  the  stresses  would  be  correspondingly  increased;  how¬ 
ever  the  general  evolution  of  the  radial  stress  component  follows  a  pat¬ 
tern  analogous  to  that  in  Figure  16  .  The  temporal  evolution  of  the 
radial  stress  component  for  a  2  mm  drop  impacting  at  600  ms  ^  is  shown  in 
Figure  17  .  The  larger  radial  distance  compared  with  the  collision  at 
200  ms  ^  prior  to  the  development  of  significant  tensile  stresses  is  ap¬ 
parent.  This  effect  would  be  reflected  in  the  larger  dimension  of  the 
central  undamaged  zone  as  the  impact  velocity  is  increased.  When  t = 45  ns 
the  plot  in  Figure  17  indicates  the  peak  radial  tensile  stress  is  0.75  p^ 
Referring  to  Figure  14  for  the  value  of  p^  ,  the  magnitude  of  the  ra¬ 
dial  stress  is  about  1  GPa.  This  stress  level  is  more  than  sufficient  to 
initiate  fracture.  The  time  t = 45  ns  corresponds  to  the  peak  pressure 
loading  as  seen  from  Figure  15  when  v^=600  ms  .  The  continued  in¬ 
crease  in  the  radial  stress  components  after  t  =  45  ns  due  to  the  continu¬ 
ous  loading  condition  inherent  in  Blowers*  analysis  has  to  be  viewed  with 
caution.  However  for  times  greater  than  45  ns  the  characteristic  form  of 
the  radial  stress  component  is  indicated  after  the  Rayleigh  wave  has 
separated  from  the  contact  zone. 


The  rapid  change  in  the  nature  of  the  radial  stress  with  depth  into 

-1 


the  target  is  shown  in  Figure  18  when  v^=600  ms  and  t  =  65  ns  .  The 
need  for  an  accurate  analysis  of  the  transient  stresses  at  shallow  depths 
below  the  surface  of  the  target  is  amply  demonstrated  for  establishing  the 
fracture  response  of  the  target  material. 


The  development  of  significant  radial  tensile  stresses  should 
correspond  to  the  initiation  of  fractures  at  the  periphery  of  the  central 
undamaged  zone  (zone  1).  This  condition  appears  to  occur  in  the  vicinity 
of  the  location  of  the  peak  pressure  on  the  surface  of  the  target.  However 
this  criterion  cannot  be  taken  as  a  general  relationship,  since  the  radial 
location  of  the  peak  pressure  is  determined  solely  from  the  shock  wave  re¬ 
lations  for  water  and  the  stress  state  in  the  target  is  dependent  on  the 
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Temporal  Development  of  the  Normalized  Radial  Stress  Component 
for  a  2  mm  Water  Drop  Impacting  a  Zinc  Sulfide  Target  at  600  m 


Variation  of  Normalized  Radial  Stress  Component  with  Depth  into  a  Zinc  Sulfide 
Target  for  a  2  mm  Water  Drop  Impacting  at  600  ms“l  when  t==65  ns. 
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stress  wave  velocities  inherent  in  the  target  material.  With  this  proviso 
the  relationship  between  the  location  of  the  peak  pressure  and  the  onset  of 
fracture  (boundary  of  zone  1)  will  be  evaluated  from  the  available  data. 

The  relevant  data  from  Tables  5  and  7  is  consolidated  in  Table  8  along 
with  the  calculated  values  of  the  critical  contact  radii  determined  from 
Figure  15  .  The  comparison  between  the  measured  values  of  the  outer  bound¬ 
ary  of  fracture  zone  1  and  the  calculated  values  of  the  location  of  the  peak 
pressure  indicates  a  fairly  close  correspondence  within  the  level  of  cumula¬ 
tive  experimental  error  in  the  measurements  required.  The  variation  in  the 
measured  dimensions  is  due  to  the  circumferential  nonuniformity  in  the  ini¬ 
tiation  of  fracture  and  the  differences  which  arise  depending  on  whether  or 
not  the  specimens  have  been  etched.  This  measure  of  the  fracture  response 
is  independent  of  the  zinc  sulfide  formulation  investigated  as  demonstrated 
by  the  comparison  in  Table  8. 


53 


TABLE  8.  Comparison  of  the  Calculated  and  Measured  Values 
of  the  Critical  Radii  for  Spherical  Drop  Impacts 
on  Zinc  Sulfide. 


Shot  No. 

Material 

Impact 

Velocity 

(ms“l) 

Drop 

Radius 

(mm) 

a  /r 
c  o 

ac 

(calculated) 

(mm) 

ac 

(measured) 

(mm) 

1036 

B 

299 

0.85 

0.141 

0.120 

0.115 

1016 

A(B-2) 

331 

.85 

.151 

.128 

.11-. 125 

1012 

A(G-2) 

336 

.81 

.153 

.124 

.125-. 130 

1017 

C 

338 

.78 

.153 

.119 

.105-.  no 

1040 

FLIR 

401 

1.01 

.172 

.174 

.180 

1041 

FLIR 

497 

1.05 

.196 

.206 

.192 

1051 

A(G-2) 

505 

1.05 

.198 

.208 

.200 

978 

C 

512 

0.84 

.199 

.167 

.15-.  16 

1052 

A(G-2) 

514 

1.07 

.200 

.214 

.200 

1048 

FLIR 

516 

1.05 

.200 

.210 

.200-. 220 

1038 

B 

517 

0.86 

.200 

.172 

.160-. 165 

1049 

A(G-2) 

518 

1.05 

.200 

.210 

.200 

1042 

FLIR 

594 

1.05 

.216 

.227 

.215-. 225 

1050 

FLIR 

600 

1.05 

.217 

.230 

.230 
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SECTION  IV 


POST-TEST  EXAIIINATION  OF  ZINC  SULFIDE  SPECIMENS 


The  general  features  of  the  grain  structure  in  CVD  zinc  sulfide  have 
been  described  in  Section  II. 1.  The  influence  of  the  grain  size,  grain 
orientation,  grain  growth  patterns,  inclusions,  surface  imperfections,  and 
hardness  on  the  nature  of  the  damage  produced  by  water  drop  collisions  has 
been  considered  here.  Overviews  of  water  drop  impact  damage  on  CVD  zinc 
sulfide  have  been  given  previously,  but  it  does  not  appear  that  a  system- 
matic  analysis  of  the  above-mentioned  characteristics  has  been  pursued 
(Adler  and  Hooker,  1976,  1978;  Graves,  et  al.,  1977;  Hackworth  and  Kocher, 
1977;  Hackworth,  et  al.,  1979;  Peterson,  1975,  1979).  Detailed  microscopic 
examination  of  both  the  impact  face  and  cross  section  of  impacts  on  the 
three  materials  described  in  Section  II. 1  over  the  velocity  range  from  237 
to  640  ms  have  been  carried  out  to  examine  the  most  damaging  aspects  of 
the  inherent  microstructure  in  CVD  ZnS. 

Variations  in  surface  finish  have  the  greatest  effect  upon  the 
resultant  extent  and  type  of  fractures  observed  for  a  given  velocity-drop 
size  condition.  The  strong  effect  of  surface  scratches  upon  the  fracture 
pattern,  especially  the  fractures  in  zone  3  (the  fracture  zones  were  de¬ 
fined  in  Section  Ill.Z.b  can  be  readily  seen  in  Figure  7b.  Zone  3  now 
contains  in  addition  to  fractures  decreasing  in  depth  with  radius  some 
large  and  deep  fractures  nucleated  at  surface  scratches. 

The  influence  of  the  surface  grain  structures  and  orientation  was 
investigated  using  scanning  electron  microscopy  (SEM)  and  transmission 
electron  microscopy  (TEM)  of  surface  replicas.  The  results  of  the  SEM 
investigation  will  be  discussed  first.  SEM  micrographs  of  shot  no.  838 
are  provided  in  Figures  19  and  20  .  The  Impact  face  of  shot  no.  838  is 
oriented  to  coincide  with  the  growth  direction;  i.e.,  along  the  major 
axis  of  the  columnar  grains  as  can  be  seen  in  the  figures. 

The  circumferential  fractures  on  the  impact  face  display  significant 
interaction  with  the  surface  grains  in  fracture  zone  2.  Figure  19  shows 
two  views  of  the  f’"actures  in  zone  2  of  the  etched  impact  face  of  shot 
no.  838  for  the  major  axes  of  the  grains  oriented  tangent  and  perpendicular 
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to  the  stress  wave  front.  The  stress  wave/grain  orientation  interaction 
controls  the  character  of  the  circumferential  fractures.  For  the  condi¬ 
tions  represented  in  Figure  19,  major  axes  of  the  grains  tangent  to  the 
stress  wave  propagation  direction,  the  fractures  are  almost  entirely  inter¬ 
granular,  When  the  grain  major  axes  are  perpendicular  to  the  stress  wave 
propagation  direction^  Figure  19b  ,  the  circumferential  fractures  tend  to 
be  shorter  than  in  Figure  19a  ,  less  continuous,  and  there  is  relatively 
little  tendency  for  the  fractures  to  follow  grain  boundaries. 

Figure  20  shows  the  radial  variation  in  the  circumferential  fracture 
concentration  in  a  region  originating  at  the  boundary  of  the  central  un¬ 
damaged  zone  (right  hand  side  of  figure)  to  the  region  at  wnich  lateral 
outflow  damage  is  initiating  (left  hand  side  of  figure).  The  composite 
micrograph  is  along  the  stress  wave  propagation  direction  which  is  from 
right  to  left.  The  major  axes  of  the  grains  are  oriented  at  approximately 
50®  to  the  stress  wave  propagation  directions.  A  higher  concentration  of 
short  circumferential  fractures  is  evident  at  the  periphery  of  the  undam¬ 
aged  zone  than  at  larger  radial  distances.  The  circumferential  fracture 
lengths  tend  to  increase  with  radial  distance.  The  onset  of  lateral  out¬ 
flow  damage  would  indicate  that  the  loading  phase  ef  the  water  drop 
collision  has  terminated. 

Two  preferential  crack  orientations  can  be  distinguished  in  the  high 
crack  concentration  annulus  (zone  2),  One  system  of  fractures  represents 
Intergranular  fractures  oriented  predominantly  along  the  grain  boundaries 
parallel  to  the  major  axes  of  the  grains.  These  fractures  are  accordingly 
directed  at  angles  from  rough].y  AO  to  60®  to  the  wave  propagation  direction 
extending  from  lower  left  to  upper  right  in  Figure  20  ,  Tl'  >.  second  system 
of  fractures  is  due  to  a  distribution  of  fine  polishing  scratches  and  are 
oriented  at  an  angle  in  the  vicinity  of  75®  to  the  wave  propagation  direc¬ 
tion  extending  from  the  lower  right  to  upper  left  in  Figure  20,  The  orien¬ 
tations  of  the  major  axes  of  the  grains  and  the  wave  propagation  direction 
represent  an  intermediate  condition  to  those  shown  in  Figure  19a  and  b. 
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As  the  radial  distance  from  the  point  of  impact  increases  the 
concentration  and  orientation  of  the  circumferential  fractures  change. 

The  orientation  of  the  circumferential  fractures  becomes  normal  to  the 
wave  propagation  direction  thereby  essentially  ignoring  the  grain  struc¬ 
ture.  The  fracture  mode  is  now  transgranular  with  little  preference  for 
grain  boundaries,  unless  the  grain  orientation  is  highly  favorable.  This 
mode  of  fracture  prevails  in  the  intermediate  fracture  zone  as  the  concen¬ 
tration  of  circumferential  fractures  drastically  decreases  and  the  fracture 
lengths  increase  significantly.  The  occasional  cracks  observed  far  from  the 
point  of  impact  (at  the  outer  boundaries  of  zone  3)  generally  are  associated 
with  some  form  of  micros true tural  defect.  Some  form  of  stress  concentration 
is  necessary  to  nucleate  and  propagate  a  crack  at  the  low  stress  levels 
which  would  prevail  in  this  region. 

The  surface  of  shot  no.  838  contained  some  polishing  scratches  which 
interacted  with  the  fracture  pattern  on  a  very  small  scale,  although  the 
gross  fracture  behavior  did  not  show  interactions  with  the  polishing  flaws. 

In  contrast  shot  no.  1012  which  was  polished  by  a  commercial  polisher  to 
above  average  optical  window  tolerances  showed  no  scratch  interaction  with 
the  fracture  pattern  on  either  a  micro  or  macroscopic  scale.  Shot  no.  1012 
was  on  a  surface  parallel  to  the  growth  plane  (i.e.,  the  columnar  grains 
were  perpendicular  to  the  impact  plane)  at  an  impact  velocity  of  336  ms 
The  impact  surface  was  examined  by  TEM  microscopy  of  a  shadowed  carbon 
replica.  Figure  21  is  a  region  within  zone  2  adjacent  to  the  central  un¬ 
damaged  zone  which  appears  in  the  lower  left  of  the  figure.  The  fractures 
tend  to  be  of  two  types:  those  that  are  small  (<ly)  and  wholly  contained 
within  the  grains  and  those  which  are  somewhat  longer  ('\'3vim)  and  are  asso¬ 
ciated  with  grain  boundaries.  This  behavior  is  very  similar  to  that  of 
Figure  19b  without  the  perturbing  influence  of  the  surface  scratches. 

Figure  22  shcws  the  fractures  at  the  inner  radius  of  zone  3  (i.e.,  adja¬ 
cent  to  zone  2)  in  which  the  fractures  are  transgranular  with  no  influence 
from  the  grain  structure  or  orientation. 

The  critical  stress  states  due  to  a  V7ater  drop  impact  on  an  elastic 
half-space  have  been  described  in  Section  III.3.b.  For  brittle  materials  the 
magnitude  and  duration  of  the  radial  stress  component  when  it  is  tensile  at 
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Figure  21. 


TEM  Micrograph  of  a  Surface  Replica  of  Impact  Face 
of  Shot  No.  1012  in  the  Region  of  the  Inner  Radius 
of  Zone  2,  Note  the  Very  Short  Circumferential 
Extent  of  the  Fractures. 


and  near  the  specimen's  surface  are  the  dominant  contributions  to  fracture 
initiation  and  propagation.  The  magnitude  and  duration  of  these  tensile 
stresses  are  related  to  the  size  of  the  drop  and  its  velocity  at  impact. 

Once  the  tensile  stresses  develop  in  advance  of  the  expanding  water  drop/ 
target  interface,  they  are  of  a  magnitude  which  more  than  exceeds  the  cri¬ 
tical  stress  intensity  factor  of  the  zinc  sulfide  target.  The  surface 
cracks  surrounding  the  central  undamaged  zone  consisting  of  the  major  por¬ 
tions  of  fracture  zones  2  and  3  show  relatively  little  interaction  with  the 
surface  grain  structure  and  near  surface  inclusions.  Figure  19a  shows  that 
in  this  region  the  grain  boundaries  provide  paths  for  the  crack  trajectories 
when  the  major  axes  of  the  grains  are  approximately  noirmal  to  the  stress 
wave  propagation  direction.  As  the  stress  waves  propagate  away  from  the 
impact  site,  they  are  attenuated  due  to  the  inherent  attenuation  of  the 
material  as  well  as  due  to  the  energy  lost  in  creating  and  extending  cracks. 
Thus  fractures  at  the  more  distant  points  in  zone  3  will  tend  to  show  the 
greatest  dependence  on  the  material’s  micros true tural  features.  The  inter¬ 
actions  between  the  fractures  found  in  this  region  and  those  features  which 
are  able  to  provide  sites  for  stress  amplification  of  the  low  amplitude 
stress  waves  of  moderately  long  duration  propagating  through  zone  3  were 
examined  in  moderate  detail.  Interactions  did  occur  with  some  of  haystack 
structures  described  in  Section  II.  1,  however  this  contribution  to  the  general 
degradation  of  the  specimen  is  essentially  Inconsequential. 

The  fracture  interaction  with  the  haystack  grain  formations  may  be  a 
result  of  the  presence  of  the  same  inclusion  that  nucleated  the  haystack  or 
could  be  a  result  of  the  grain  orientation  interface.  Wiiraner  and  Graves 
(1977)  note  than  the  columnar  grain  structure  results  in  some  anisotropy  in 
strength,  hardness,  and  Poisson’s  ratio.  Their  strength  and  elastic  property 
tests  were  on  a  macroscopic  scale,  so  it  is  conceivable  that  the  effects  of 
local  anisotropy  could  be  much  greater  than  their  mechanical  property  eval¬ 
uations  indicate.  The  local  anisotropy  of  the  CVD  zinc  sulfide  could  in 
itself  produce  critical  stress  conditions  for  fracture  initiation. 


The  major  crack.;  for  specimens  v;ith  a  moderate  distribution  of  fine 
polishing  scratches  are  relitively  large  semi-elliptical  fractures  nu¬ 
cleated  at  surface  scratches  as  can  be  seen  in  Figure  7b,  These  frac¬ 
ture  surfaces  exhibit  a  large  orientational  dependence  between  the  stress 
wave  direction  and  the  surface  scratch.  Surface  scratch  nucleated  cracks 
are  typically  bounded  in  extent  by  '"he  points  at  which  the  included  angle 
reaches  45°  between  the  line  of  the  scratch  and  the  tangent  to  an  expanding 
circle  centered  at  the  impact  site  intersecting  the  scratch.  Even  if  the 
scratch  is  composed  of  a  series  of  small  disjointed  segments  oriented  dif¬ 
ferently  than  the  composite  scratch  directions  the  resultant  fractures  will 
only  interact  with  these  segments  very  weakly,  if  at  all,  until  they  satisfy 
the  included  angle  requirement  described  above.  The  scratch  -lucleated  frac¬ 
tures  tend  to  be  more  pronounced  in  fracture  zone  3  than  in  fracture  zone  2, 
This  may  be  due  to  the  longer  duration  but  decreased  amplitude  of  the  stress 
pulse.  The  decreased  stress  amplitude  is  compensated  by  the  ease  of  fracture 
nucleation  at  a  fine  surface  scratch. 

In  addition  to  the  influence  of  grain  boundaries  and  surface  polishing 
scratches  on  impact  fracture  behavior,  the  influence  of  grain  size,  inclusion 
content,  and  hardness  were  also  investigated. 

Three  different  formulations  of  ZnS  (materials  A,  B,  and  C)  as 
described  in  Section  II. 1, whose  properties  are  tabulated  in  Tables  1  to  3, 
provided  variations  in  grain  size  and  hardness.  Material  C  has  both 
a  lower  microhardness  and  a  larger  grain  size  than  both  material  A  of  B 
which  are  similar  in  grain  size  and  hardness.  Comparison  of  the  impact  dam¬ 
age  due  to  single  drop  impacts  on  these  materials  reveals  no  significant 
difference  in  the  extent  or  nature  of  the  fractures  produced.  Cross  sec¬ 
tional  views  of  the  fracture  patterns  are  shown  in  Figures  6  and  23  for 
impact  velocities  of  517,  515,  and  512  ms  on  materials  B,  A,  and  C,  re¬ 
spectively.  It  is  apparent  from  these  figures  that  the  variations  in  hard¬ 
ness,  grain  size,  grain  shape,  and  inclusion  content  represented  by  these 
three  materials  does  not  significantly  affect  the  impact  damage.  Comparisons 
of  the  measurements  of  the  fracture  zone  boundaries  and  fracture  depth  dimen¬ 
sions  tabulated  in  Tables  6  and  7  corroborate  this  conclusion. 
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In  order  to  extend  the  range  of  available  grain  structures,  specimens 
were  also  cut  from  material  A  parallel  to  the  deposition  direction  to  yield 
specimens  in  which  the  exposed  face  has  the  elongated  grain  structure  illus¬ 
trated  in  Figut':'  la  .  Comparable  water  drop  Impact  conditions  (shot  no.  813 
and  838),  representing  the  difference  in  grain  structure  shown  in  Figure  1  , 
produced  no  obvious  difference  in  the  fracture  pattern  observable  on  both  a 
microscopic  and  macroscopic  scale. 
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SECTION  V 

ANALYSIS  OF  WATER  DROP  IMPACTS  ON  ELASTIC  SOLIDS 


Knowledge  of  the  stress  states  in  a  target  material  due  to  a  water 
drop  collision  is  an  essential  prerequisite  for  analyzing  the  extent  of 
the  resulting  crack  growth  which  occurs  when  the  impact  velocity  exceeds 
the  fracture  threshold.  In  general,  the  motion  of  the  drop  and  the  mate¬ 
rial  is  a  coupled  problem  and  should  be  treated  as  such,  since  the  contact 
pressure  is  a  function  of  the  deceleration  of  the  drop  which  depends  in 
turn  on  the  deformation  of  the  surface.  This  leads  to  a  complex  non-linear 
interdependence  between  the  pressure  and  the  surface  deformation.  We  will 
side-step  this  issue,  however,  by  assuming  that  for  the  purpose  of  calcu¬ 
lating  the  pressure,  the  material  surface  can  be  assumed  to  be  rigid. 

This  will  result  in  an  c  erestimation  of  the  pressure  which  for  materials 
with  the  hardness  of  zinc  sulfide  is  not  too  significant.  This  can  be 
corrected  if  desired  by  using  an  iterative  procedure  once  the  general 
solution  for  the  surface  deformation  has  been  obtained. 

The  analytic  approach  developed  here  is  based  on  the  use  of  Green's 
functions  and  a  Green's  theorem  relation  derived  from  the  divergence  theorem. 
For  these  techniques  to  be  useful,  it  is  necessary  that  the  problem  be  linear 
which  imposes  some  restrictions  on  the  drop  speed  and  material  properties. 

In  pa-ticular,  deformations  must  be  sufficiently  small  that  a  linear  stress- 
strain  relation  is  valid.  This  eliminates  from  consideration,  for  example, 
those  events  which  result  in  substantial  cratering.  The  linearity  assumption 
is  bound  to  fail  at  crack  tips  however,  since  a  strictly  linear  material 
would  never  crack  and  consequently,  some  special  treatment  of  these  regions 
is  required.  From  the  ETI  experimental  work  on  ZnS,  it  appears  that  in  the 
velocity  range  of  interest,  namely  222  to  684  ms  the  response  of  the 
material  is  reasonably  linear  except  for  the  creation  of  the  cracks.  This 
is  deduced  from  the  fact  that  any  permanent  deformations  of  the  surface  as 
a  result  of  the  impact  are  negligible  compared  to  the  dimensions  of  the  impact 
zone.  For  simplification,  the  material  will  also  be  assumed  isotropic. 

This  assumption  is  supported  by  experiments  on  ZnS  conducted  by  ETI  which 

do  not  show  any  strong  dependence  of  impact  damage  on  grain  orientation 

for  the  typical  grain  dimensions  in  the  current  production  of  CVD  ZnS  windows. 
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PREVIOUS  WORK 


A  few  connnents  concerning  some  previous  work  on  this  problem  will  be 
made.  Rosenblatt,  et  al»  (1976,1977),  have  developed  a  numerical  solution 
to  the  crack  growth  problem  in  ZnSe  using  explicit  numerical  methods.  In 
comparing  their  results  with  the  experimental  results  obtained  in  this  pro¬ 
gram  for  ZnS,  it  is  seen  that  the  numerical  solution  produces  cracks  whose 
special  extent  essentially  encompasses  the  annulus  corresponding  to  zone  2 
cracks  (ETI  notation)  observed  in  the  cross  sections  in  Section  IV.  The  cracks 
in  the  numerical  solution  are  not  as  clearly  defined  as  in  the  experiments: 
they  tend  to  appear  more  in  the  form  of  a  completely  cracked  zone  than  as 
a  distribution  of  distinct  cracks.  Also,  the  experimentally  observed  zone  3 
cracks  appear  to  be  almost  completely  absent  in  the  numerical  solution. 

Since  ZnS  has  a  greater  strength  than  ZnSe,  this  discrepancy  is  not  likely 
to  be  explained  by  the  difference  in  the  materials.  A  more  likely  explana¬ 
tion  is  that  the  crack  model  used  by  Rosenblatt  and  co-workers  is  incomplete, 
since  according  to  their  fracture  criterion  only  tensile  cracks  are  allowed. 
Cracks  can  occur  in  any  cell  in  the  finite  difference  grid  once  the  tensile 
stress  criterion  is  satisfied.  This  is  not  an  appropriate  representation 
for  the  preferential  crack  formations  controlled  by  the  imposed  transient 
stresses  and  a  surface  flaw  distribution.  Considerable  improvements  in  the 
basic  fracture  model  appear  to  be  required  before  the  finite  difference  com¬ 
putations  can  begin  to  reproduce  the  physics  of  the  actual  impact  induced 
fracture  patterns.  This  points  to  an  area  of  concern  which  requires  further 
study  and  points  out  the  danger  of  relying  too  heavily  on  numerical  results 
alone. 


A  second  point  examined  by  Rosenblatt,  et  al.  (1977)  was  the  pressure- 
surface  deformation  interdependence.  It  was  found  that  assuming  a  rigid 
surface  for  the  purposes  of  calculating  the  pressure  does  not  introduce 
intolerable  errors  over  the  velocity  range  of  interest. 
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The  analytic  approach  used  by  Hackworth  and  Kocher  (1978)  will  also 
be  mentioned,  since  their  analysis  was  to  provide  the  local  stress  condi¬ 
tions  at  various  locations  within  the  half-space  for  the  purpose  of  esti¬ 
mating  the  crack  growth  which  would  occur.  However  detailed  examination 
of  their  work  revealed  an  apparent  error  in  their  contribution  to  an  ex¬ 
tension  of  the  initial  analysis  of  the  stress  wave  problem  (Blowers,  1969; 
Adler,  1977).  Their  starting  point  is  the  TURBAN  computer  program  which 
evaluates  the  stresses  in  an  elastic  half-space  due  to  a  pressure  distri¬ 
bution  given  by  P(r,t)  =  P^6(k>/t^r)  where  0(x)  is  the  Heaviside 
function 


0(x)  = 


1  x>0 


0  x<0  . 


(5.1) 


The  analytic  solution  was  determined  by  Blowers  (1969)  for  this  pressure 
distribution  iu  terms  of  several  integrals  which  are  evaluated  by  TURBAN. 
Hackworth  and  Kocher ’s  extension  of  Blowers'  idealization  of  the  water 
drop  impact  problem  to  time-dependent  pressure  loadings  is  based  on 
Duhamel's  superposition  relations  which  utilized  Blowers  solution,  0^  . 
They  proposed  that  0(r,z;t)  ,  the  solution  for  the  general  time-dependent 
distribution  P(r,t)  ,  is  given  in  terms  of  the  solution,  a^(r,z;t)  ,  for 
the  restricted  distribution  expression 


P  (t 

a^(r,z,t-T)  dx 


(5.2) 


where  the  final  pressure  distribution  is  given  by 


P(r,t)  =  P  (r)  P(t)  . 
o 


(5.3a) 


^  f 


(5.3b) 
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This  is  correct  as  far  as  it  goes,  but  Eq.(5.2)  is  valid  only  for  the 
limited  class  of  distributions  f. _  which  the  restricted  distribution  is 
Independent  of  time,  a  point  which  is  not  clearly  made  in  the  report. 

This  however  excludes  Blowers'  solution  upon  which  TURBAN  is  based.  A 
more  general  relation  similar  to  Eq.(5.2)  can  easily  be  found,  however, 
which  is  valid  for  Blowers'  distribution  although  the  final  pressure  dis¬ 
tribution  does  not  have  the  simple  product  form  sought  by  Hackworth  and 
Kocher.  The  general  solution  as  will  be  shown  later  (see  Eq.(5.50))  for 
the  physical  stress  matrix  at  location  (y,t)  in  a  half-space  due  to  a 
prescribed  pressure  distribution  on  the  surface  is  given  by 

ay(y.t)  P(x,t-T)  E^^Cx.y;!)  (5.4) 


where^da^  denotes  the  integral  over  the  surface  S  ,  specified  by  X2=0, 
of  the\alf-space.  On  this  surface  the  coordinate  vector  is  x  ,  the  pres¬ 
sure  distribution  is  P(x,t)  and  the  appropriate  Green's  tensor  is 
E3ij(x,y;T)  .  Taking  the  Laplace  transform  of  this  expression,  we  find 


5^j(y,s)  -i  da^P(x,s) 

s 


23ij(x,y;s) 


(5.5) 


Since  the  problem  is  axlsymmetric,  Eq.(5.5)  can  be  written  in  cylindrical 
coordinates.  Thus  (dropping  the  tensor  indices  and  incorporating  numerical 
constants  into  the  definition  of  P(r',s)) 

00 

a(r,z;s)  =J  dr'  P(r',s)  Z(r',r,z;s)  (5.6) 

o 

In  order  to  reproduce  the  result  given  by  Hackworth  and  Kocher  (1978),  we 
will  first  assume  a  pressure  distribution  which  is  independent  of  time. 
Later,  we  will  extend  these  results  for  a  general  pressure  distribution. 
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Suppose  we  know  the  solution  for  a  constant  pressure  •  Then 


P^(r  .=)  ./  P^(r)dt  - 


(5.7) 


and 


a(r,z,s)  -  a  (r,z,s)  =  f  dr'  E(r',r,z;s)  1  - 

O  •r  ^ 


rsP(r'.s)-P  (r') 


■]<3. 


8) 


Since 


P(r,s)  =  =  /  e  l^t  =  sP(r,s)  -  P(r,o) 


-St  dP 


(5.9) 


and  provided  that  Pq(’^)  =  P(r,o) 


a(r,z,s)  -  o^(r,z,s)  =J  dr'  E(r',r,z;s) 


P  (r’.s) 
s 


r“  -  I  (r'  s) 

=  /  dr'E(r',r,z;s) 


P  (r’) 

o 


(5.10) 


The  Laplace  transform  of  Eq.(5.2)  gives 


-  -  -  P(s) 

o(r,z,s)  -  a^(r,z,s)  =  o^(r,z,s) 


(5.11) 


Comparing  the  expressions  in  Eq.(5.10)  and  (5.11),  it  is  seen  that  Eq.  (5.2) 
can  only  be  valid  for  those  pressure  distributions  for  which 


P(r'.s) 

Po(r') 


is  independent  of  r'  . 


This  requires 


P(r,s)  =  P^(r)H(s)  .  (5.12) 

where  H(s)  is  an  arbitrary  function  of  s  . 

In  view  of  Eq.(5.9),  we  find 
P(r,.) 

and  it  follows  that 

P(r,t)  =  P^(r)  F(t) 

where 

F(t)  =  1+  j  H(T)dT  . 

0 

Thus,  Eq. (5.2)  is  only  valid  for  restricted  pressure  distributions  P  (r) 

o 

which  are  functions  only  of  r  .  This  excludes,  for  example.  Blowers’ 
distribution  since  it  depends  on  both  r  and  t  .  It  is  natural  to 
ask  then  whether  a  product  form  might  in  fact  be  more  useful  than  Blowers' 
solution  since  then  one  would  be  able  to  include  the  time  dependence  of 
the  pressure  at  any  given  point.  This  is  not  likely  to  be  the  case,  how¬ 
ever,  since  such  a  distribution  does  not  model  the  expansion  of  the  con¬ 
tact  zone  of  the  drop.  Since  the  waves  are  expanding  outward  from  the 
contact  zone,  the  fact  that  the  source  is  expanding  is  certainly  as  impor¬ 
tant  as  the  fact  that  its  magnitude  changes  with  time  at  a  fixed  point. 


(5.13) 


(5.14) 


(5.15) 
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These  results  are  easily  extended,  however,  to  general  distributions 
P^(r,t)  which  do  include  Blowers’  distribution  as  a  special  case,  although 
the  resulting  final  pressure  distributions  are  not  simply  the  products  of 
the  restricted  distributions  and  arbitrary  functions  of  time. 

If  we  know  the  stress  state  due  to  a  pressure  distribution  P^(r,t)  , 

then 


a(r,z;s)  -  a^(r,z;s)  =J’  dr'  E(r ' ,r,z;s)  [  P(r ' ,s)  -  P^(r',s)] 

°  (5.16) 


CO 

f  dr'  I(r',r,z;s)  P  (i',s)  x 

'p(r'.s) 

J  0 

0 

P  (r’,s) 

0 

As  before,  the  right  hand  side  of  this  relation  can  be  factored  provided 


-  1  =  H(s)  (5.17) 

P  (r',s) 
o 


in  which  case 


a(r,z;s)  -  a^(r,z;s)  =  a^(r,z;s)  H(s) 


or 


From  Eq.(5.17) 


(5.18) 


a(r,z;t)  =  0^(r,z;t)  +  J  ^  dT  a^(r,z;t-T)  H(t)  (5.19) 


P(r,s)  =  PQ(r,s)  |l+H(s)j 


(5.20) 


which  implies 


P(r,t)  =  P„(r,t)  +f^  P„(r,t-T)  H(t)  dr 
o  »  o 


(5.21) 
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It  is  easy  to  show  that  this  reduces  to  the  previous  result  when 
PgCrjt)  is  independent  of  time.  In  the  general  case,  H(t)  is  not 
equal  to  the  normalized  time  derivative  of  the  pressure  function  nor 
is  the  final  pressure  a  simple  product  of  the  restricted  pressure, 
r^(r,t)  ,  and  an  arbitrary  function  of  time.  For  example,  using 
Blowers'  distribution 


P^(r,t)  =  P^e(ki^-r) 


(5.22) 


where  P^  is  a  constant.  We  find  after  some  manipulation. 


P(r,t)  =  Pj,(r,t)  [1+I(r,t)] 


(5.23) 


where 


.t- 


I(r,t)  *  9(k/t-r)  J  ^  dt  H(t)  . 
o 


(5.24) 


This  does  not  have  the  simple  product  form  sought  by  Hackworth  and  Kocher 
because  of  the  radial  dependence  of  I(r,t).  Using  the  fact  that 


/  2  . 

0(k/t-r)  = 


(5.25) 


we  notice  that  Eq.  (5.23)  can  be  rewritten  in  the  form 


P(r,T)  =  P  (r,T)  [l+I(T)] 


(5.26) 


where 


T  =  t  -  — 
k' 


(5.27) 


P  (r,T)  =  P  0(T) 
o  ’  o  ' 


(5.28) 
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Thus  it  is  apparent  that  within  the  confines  of  Eq.(5.21)  generalizations 
of  Blowers*  problem  can  be  generated  through  the  use  of  Eq.(5.23).  If  our 
only  interest  was  the  stress  state  for  the  axisymmetric  problem  and  this 
rather  restricted  class  of  pressure  distributions,  we  would  in  fact  use 
these  results  without  introducing  the  additional  complexity  of  the  Green’s 
tensor.  The  latter  is  central  also,  however,  in  developing  a  systematic 
treatment  of  wave  diffraction  from  cracks  and  for  treating  asymmetric  pres¬ 
sure  distributions.  An  example  of  the  latter  is  the  distribution  due  to  a 
drop  impact  at  an  angle  of  incidence  other  than  90®. 

2.  A  FOKMAL  SOLUTION  TO  THE  PROBLEM  OF  CALCULATING  THE  STRESS  DISTRIBUTION 
IN  AN  ELASTIC  HALF-SPACE  DUE  TO  AN  ARBITRARY  PRESSURE  DISTRIBUTION  ON 
ITS  SURFACE 

The  formal  solution  for  the  determination  of  the  stresses  in  a  half¬ 
space  due  to  a  localized  pressure  distribution  on  its  surface  is  formulated 
in  terms  of  a  Green's  theorem  derived  from  the  divergence  theorem.  The  de¬ 
rivation  of  this  relation  will  be  outlined  because  a  clear  understanding  of 
its  origin  will  be  necessary  later.  Consider  two  solutions  to  the  wave  equa¬ 
tion  in  the  half-space  denoted  by  (Ujj^,tj^^)  and  (u^*,t^j^)  where  u^  and 
denote  the  displacement  vector  and  stress  tensor,  respectively*  (The 
physical  stress  matrix  is  not  a  tensor  because  it  does  not  transform 

as  such  for  curvilinear  coordinates.  If  needed,  the  relation  between  tj^^ 
and  is  easily  found  using  the  metric  tensor  for  the  coordinates 

involved.  For  cartesian  coordinates,  this  tensor  is  unity.) 


Form  the  quantity 


(5.29) 


where  *  denotes  the  convolution  defined  by 

t 

a*b  =  dr  a(t-T)  b(T)  .  (5.30) 

0 
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The  upper  limit  t  denotes  the  limit  Lim  t+€  whenever  it  is  necessary 

€-H3 

to  avoid  ambiguity  in  the  definition  of  the  integral.  In  particular. 


f  6(t-T)  b(T)  dt  =  b(t) 


(5.31) 


where  6(t)  is  the  Dirac  delta  function.  A  repeated  subscript  such  as 
in  Eq,(5.29)  indicates  a  summation  over  the  range  1  to  3  unless  otherwise 
indicated. 

Taking  the  divergence  of  Eq.(5.29)  and  using  the  equations  of  motion. 


*^k£,k 


(5.32) 


*^k£,k 


(5.33) 


(The  notation  a  denotes  the  partial  derivative  .)  we  find, 

5  It 


w  =  t  *e’  -  t'  *6  +  oti  *u'  -  dU'  *u 


-  ‘’'tS  + 


(5.34) 


where  the  strain  components  are  defined  by 


'k»  =  "k.»^ 


(5.35) 


^k£*®k£  “  ®k£mn®mn*®k£  "  ®k£*®k£mn®m  "  ®k£*^k£  “  ^k£*®k£ 
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where  we  have  used  the  general  stress-strain  relation  t,  „=  E,  .  e 

®  k)l  kJtmn  mn 

and  the  symmetry  property  ~  ^mnkil  elastic  constants. 

Also,  from  the  definition  of  the  convolution,  Eq.(5.30),  it  is  easily 


seen  that 


u*u'  =  (u*u')  -  v(o) ‘u’  -  u(o)*u* 


(5.37) 


-y- 

where  v(o)  =  u(o)  .  A  similar  expression  holds  with  the  primed  and 
unprimed  variables  reversed.  Thus 


Wk,k=  V-w  = 


-pf*u'  +  pf'*u  -  pv(o)*u’  +  pv'(o)*u 


4- 

pu(o) ’U*  +  pu* (o) ‘U 


(5.38) 


The  divergence  theorem  states  that  for  any  closed 
surface  S  surrounding  a  volume  V  , 


J^V*w  dv  =  w*n  da 


(5.39) 


where  n  denotes  the  unit  outward  normal  vector  to  the  surface.  Thus 


J* da  t^^j*u'  +  J* dv  p[f*u'  +  v(o)*u'  +  u(o)*u'] 


da  t^^j*u  +  J  dv  pff '*u  +  v’(q)'  u  +  u*  (o)  'u  j 


(5.40) 


where 


(5.41) 


77 


In  a  later  discussion  it  will  be  important  that  this  relation  is  true 
for  any  closed  surface  in  the  material  and  not  just  for  the  actual 
physical  surface  of  the  material. 

Now  for  u'  ,  we  chose  the  solution  of  the  wave  equation  with  the 
singular  source  term 


pf^  =  6(t)  6(x-1)  , 


(5.42) 


where 


is  the  Kronecker  delta,  and  the  initial  conditions 


u'  (o)  =  V  ’  (o)  =■  0  ^ 


(5.43) 


We  denote  this  solution  by 


Gij(x,t;t) 


and  require 


=  0  for  t<0  . 


(5.44) 


The  term  involving  f  is  easily  evaluated 

t 


dv^  p(f'*u)j  =  ^^dv^6(x4)  /  dT6(t-T)  u^(x,t) 


=  /  dv  (S(x-t;  u  (x,t) 


u  (';,t) 

j 


5  e  V 
C  £  V 


(5.45) 


where  V  denotes  the  complement  of  V  ,  i.e..  the  region  outside  the 
boundary  of  V  .  We  will  define  u^(5>t)  on  the  surface  as  the  limit  as 
the  surface  is  approached  from  the  inside.  Thus  using  Eq. (5.40)  (dropping 
the  vector  notation  for  x  and  ?  ) 
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/ 


Uj(?,t)  =  (x,t)*G^^(x,?;t)  -  Uj  (x,t)*T^^j^j(x,C;t)] 

+  f  dv^  p|^f^(x,t)*G^^(x,  5,t)  +  v^(x,o)  G^j(x,C,t) 

+  u^(x,o)  •  G^j(x,C;t)J  (5.46) 

where  the  stress  matrix  T,  x..  =  T...  n.  is  obtained  from  G,.  ,  by 

(n)ij  ijk  k  xT 

using  the  stress-strain  relation. 


"^ijk  '^ik  ^kj,i^ 


(5.47) 


The  partial  derivatives  are  with  respect  to  the  x  coordinates. 


Since  both  t,  ..  and  u.  cannot  be  independently  specified  on 
(n) 1  i  ^ 

a  closed  surface,  we  must  impose  an  additional  constraint  on  the  Green  s 
function  to  eliminate  the  dependence  on  one  or  the  other  of  these  quanti¬ 
ties.  Since  we  wish  to  specify  the  pressure  on  the  surface,  the  natural 
choice  is  to  require 


^(n)ij 


=  0  on  S 


(5.48) 


where  S  denotes  the  surface  of  the  half-space.  If  we  further  assume 
the  absence  of  sources  and  homogeneous  initial  conditions,  then 


(5.49) 


This  is  the  fundamental  solution  to  the  half-space  problem  with  a 
prescribed  stress  distribution  on  the  surface.  The  stresses  in  the 
material  are  found  by  applying  the  stress-strain  relation  to  both 
sides  of  Eq.(5.49).  Thus 

tk^(5,t)  =  ^  da^  ‘^(n)i^^’^^  * 

s 
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where  is  determined  by  Eq.(5.47)  with  the  exception  that 

the  partial  deriratives  are  now  with  respect  to  the  5  coordinates. 


These  expressions »  while  providing  a  formal  solution  to  the 
problem  are,  of  course,  of  no  practical  value  until  an  explicit  ex¬ 
pression  for  the  Green's  tensor  is  found. 

3.  DERIVATION  OF  THE  GREEN'S  TENSOR 


The  Green's  tensor  we  require  is  the  solution  to  the  wave  equation 


~^t  ^  Sj,£i  “  ^ij]  “  ~ 


Ij 


(5.51) 


with  the  source  term 


pf^j  =  6(t)  6(x-t) 


(5.52) 


which  has  homogeneous  initial  values  and  satisfies  the  boundary 
condition  on  the  surface  of  the  half-space  given  by 


’^(n)ij  “  °  • 


(5.53) 


The  surface  stress  tensor  i®  given  in  terms  of  the  Green's 

tensor  by 


•(n)ij  “  ^%j,k‘^kj,i^]  \ 


(5.54) 


where  X  and  y  are  the  Lame  coefficients  which  are  related  to  the 
longitudinal  and  transverse  wave  speeds  by 


2.JL 

:  P 

2  X4-2y 


(5.55) 
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Since  the  governing  equations  are  linear,  we  can  always  express 
the  Green's  tensor  as  a  sum,  each  term  of  which  satisfies  the  wave  equa¬ 
tion  but  not  the  complete  set  of  boundary  conditions.  By  making  a  judi¬ 
cious  choice  in  this  separation,  the  problem  of  determining  reduces 

to  several  simplier  problems  whose  solutions  when  added  determine  G^^  . 

A  familiar  example  from  electrostatics  is  the  problem  of  determining  the 
electric  field  due  to  a  point  charge  near  an  infinite  conducting  plane. 

The  solution  is  most  easily  found  by  writing  it  as  the  sum  of  the  fields 
due  to  the  charge  alone,  which  will  not  satisfy  the  boundary  conditions 
at  the  plane,  and  that  due  to  a  virtual  image  charge  located  behind  the 
plane  with  properties  chosen  so  that  the  sum  of  the  two  fields  does  satisfy 
the  boundary  conditions  on  the  plane.  The  first  of  these  fields  will  be 
singular  at  the  location  of  the  charge,  whereas  that  of  the  second  will  be 
everywhere  regular  since  the  image  charge  does  not  lie  inside  the  physical 
region. 


We  will,  in  fact,  use  this  analogy  to  write  the  Green's  tensor  as  a 
sum  of  three  fields  with  the  first  being  the  singular  field  due  to  the 
source  Eq.(5.52)  with  no  boundaries.  The  second  term  will  be  the  field 
due  to  the  image  of  Eq.(5.52)  and  the  third  that  which  is  required  to 
satisfy  the  boundary  conditions  on  the  p]ane.  T’-ese  iatlec  two  terns  are 
regular  everywhere  inside  the  material.  In  view  of  Eq.(5.53),  the  boundary 
condition  for  the  third  term  is  just  the  negative  of  the  sum  of  the 
contributions  from  the  first  two  terms. 

Let  us  now  choose  a  definite  coordinate  system.  The  material  will 
be  assigned  to  the  half-space  x  <0  and  the  origin  will  be  placed  directly 
above  the  source  point  |  ,  Figure  24.  The  solution  for  an  arbitrarily 
located  source  point  can  be  found  by  a  simple  translation  of  coordinates. 
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where 


The  field  due  to  the  source  given  in  Eq.(5.73)  is  well  known, 

1 


g(s)  ^ 
ij  Airp 


3r.r. 
1  1 


X(S(t-Xr)dX 


r  .r. 
+  -^ 


-4  6(t-^)  -  -4  «(t-— ) 


— il  — )  1 


6.  . 


+  ^  6(t  ^ 

rct  t 


(5.56) 


^  ^  ^ 

r  =  X  -  ^  .  The  integral  is  easily  evaluated 
1 


X6(t-Xr)dX  =  -^ 

9(t-4-)  -  9(t-^) 
c 

1 

1  t 

(5.57) 
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«e  see  that  is  singular  at  fO  where  it  varies  as 

tensor  corresponding  to  this  field  is  found  by  substituting 


1/r. 

into 


The  stress 
Eq.(5.54). 
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The  second  contribution  to  the  Green's  tensor  will  be  chosen  as  the 

field  due  to  the  image  of  the  source  Eq.(5.52)  found  by  reflecting  in  the 

X2=0  plane.  The  reason  for  including  this  field  explicitly  rather  than 

combining  the  second  and  third  terms  is  that  when  the  combined  fields  of 

the  source  and  its  image  are  substituted  into  the  boundary  condition 

Eq.(5.53),  a  number  of  terms  cancel.  The  fields  due  to  the  image  are  given 

by  Eqs.(5.56)  and  (5.58)  when  |  is  replaced  by  i.e., 

(0,0,-O  and  t- 
’  ’  image 


source 

=  (0,0j5)  with  .  Thus,  since  r=x-5 


r  =  (x, ,x_,Xo  +  O 

source  1’  2’  3 


(5.60) 


'image  * 


On  the  boundary,  ,  we  find 


r  =  r.  i=l,2 

source^  image^  ’ 


source. 


-  r 


image. 


(5.61) 


source' 


|r.  I 

‘  image ' 


The  combined  contribution  to  the  boundary  condition  is  found  by  adding 
the  individual  contributions.  From  Eqs.(5.59)  and  (5.60),  it  is  apparent 
that  all  terms  will  cancel  which  are  proportional  to  an  odd  power  of  r^  . 
By  contracting  Eq.(5.58)  with  fi=(0,0,l)  ,  we  easily  find  that 


(s)  .  M 

(n)ij  (n)ij 


=  0 


for  i,j=l,2  or 


i=j=3  . 


(5.62) 
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When  i  =  3,  j  f  3  or  j  =  3,  i  3,  the  sum  is  not  zero.  For  example^ 
for  i  =  1,2, 


t(8)  +  1  f  *i] 

^(n)i3^^(n)i3-2^f^\  ^  ^5 - sl  J 


X6(t-Xr)dX 


[•^41(»-->4 

L  r  r  J  \  t  c. 


-  (^)  - -i  L(t_^)  6(t-^)|  J 

‘'i  ‘'i  r^  ‘^t  ‘'t  ‘'t  J  ) 


(5.63) 


Each  term  is  of  the  form  x^/r  multiplied  by  either  a  delta  function  or  its 
derivative.  The  problem  of  finding  the  Green's  function  thus  reduces  to 
finding  solutions  to  the  homogeneous  wave  equation  which  satisfy  these  two 
types  of  boundary  conditions.  When  the  individual  solutions  are  found, 
it  only  remains  to  chose  the  linear  combination  which  satisfies  Eq. (5.53) 
to  complete  the  derivation  of  the  Green's  tensor. 

We  seek  solutions  then  to  the  following  problems 


(5.64) 


T,  -  0  i,i  =  1,2  or  i  =  1  =  3 


[•1«1 


@ir 

p?.  .' 


m 


Mm 


;«a 


r 

where  F(t - )  =  <  . 

(6(t-f) 


(5.66) 


Notice  that  this  problem  can  be  treated  by  isolating  each  of  the  inhomo¬ 
geneous  boundary  conditions.  For  example,  one  such  problem  would  be  set 
by  requiring  =  0  for  all  i  and  j  except  for  ~  • 

While  these  problems  do  not  appear  to  have  been  solved,  the  solutions  to 
similar  problems  do  exist. 


First,  we  notice  that  since  for  i,j=l,2  and  i=j=3  and 

there  are  no  sources,  we  can  set 


G. .=0  i,j=l,2  and  i=j=3 

ij 


(5.67) 


The  nonzero  components  of  are  then  and  G^^,  j=l,2. 

We  will  separate  the  problem  into  two  parts  corresponding  to  these  two  cases. 
The  first  will  be  treated  explicitly;  the  second  can  be  treated  in  an  anal¬ 
ogous  manner.  In  order  to  simplify  the  notation,  we  define 


u.  =  G . „ 
1  i3 


*^ik  '^ik3 


9\\ 

‘^(n)i  "  ^n)i3  =  “^ik  +  +  ^k 


(5.68) 


(5.69) 


(5.70) 
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We  define  cylindrical  coordinates  (p,(p,z)  such  that  (see  Figure  24) 


X  =  (pcoscj),  psin(j),z)  (5.71) 

t  =  (0.0,-5)  (5.72) 

It  follows  that 

5  9  2  2  2  2 

r  =  x+C  +2zg  =  p  +z  +C  +2zC  (5.73) 

The  tmit  vector  normal  to  the  surface  is 

ft  =  (0,0,1)  .  (5.74) 

The  problem  now  becomes 

•• 

2  9m.  9  2  -f 

c/v  U  +  (c/-Cj.  )  V  (Vu)  -  u  =  0  (5.75) 


^  ^  /  3u  3u  \ 

=  ^13 

with  the  boundary  conditions 


*1  p  F(t-f-) 

1  1  cos<p  i=l 

Y  1 

(5.77) 

*"(11)  i 

A  / 

1  sinj)  i=2 

(5.78) 

*'(n)3 

=  0 

(5.79) 

z=0 
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/ 


Using  Eq.(5.76),  Eqs.(5.77)  to  (5,79)  become  in  rectangular  Cartesian 
coordinates 


/  3u  3u  \  A, 

(  —  +  —  =  -^  -P.9.qs<^,  n  F  ft- 

z=0  \p  +C  / 


(5.80) 


( 


3u  3u  \ 
_i  +  _JL 
3z  3y  / 


_  _1  psinij) 
z=0  ^  (p^+5^ ) 


n  F 
7 


(t-f) 


(5.81) 


/ 3u  3u  \ 

\  8x  3y  / 


+  (X+2y) 


3u 
_ z 

dz 


z=0 


=  0 


z=0 


(5.82) 


We  now  express  u  in  terms  of  a  scalar  and  a  vector  potential 


u  •••  Vq  +  V  X 


(5.83) 


where 


V.i?=0 


(5.84) 


The  wave  equation  then  separates  into  equations  for  q  and  ij^ 


2  1  •* 

V  q  -  q  =  0 

c  ^ 


(5.85) 


2-v  1  -> 

V  il>  -  -^2  -^  =  0 


(5.86) 


The  problem  thus  becomes  that  of  solving  Eq,(5.85)  and  (5,86)  with  the 
subsidiary  condition  Eq,(5,84)  and  the  boundary  conditions  Eqs.(5.80) 
to  (5.82). 
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By  using  the  Laplace  transform,  we  can  eliminate  the  time  dependence 
in  these  equations.  Define 


Ti(s)  »  /  e“®^  n(t)  dt 


■/ 


(5.87) 


=J^ e“®*^  ij)^(s)dt 


(5.88) 


/ 


F(r,s)  =  /'*"  e“®*^  F(t-  ^)  dt 


(5.89) 


The  latter  transform  can  be  given  explicitly  using  the  definition  of 

F(t-  f)  , 

rs 


e  ^  for  F(t-  ^)  =  6(t- 


F(r,s)= 


rs 


se  -  c6(r) 


for  F(t-  |)  =  6(t- 


(5.90) 


Further, 


n(s)  =  s^Ti(s)  -  f|(0)  -  sn(0) 


^^(s)  =  s^l^(s)  -  v^(0)  -  S1|;^(0) 


(5.91) 

(5.92) 


Using  the  initial  conditions 


Tl(0)  =  ri(0)  =  i^'^(O)  =  i{»^(0)  =  0 


(5.93) 
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''7  '■'"■i'.J"'- 


/ 


Tli6  wav6  GcjustionSj  E(is*(5*85)  and  (5*86)j  bacoinc  (omitting  the  bair 
notation  for  transformed  quantities) 


(V^-k^^)Ti  =  0 

(V^-k^^)l  =  0 


(5.94) 

(5.95) 


where 


=  — 


=  — 

2  c. 


(5.96) 

(5.97) 


These  equations  written  in  cylindrical  coordinates  become 


(V^-k^^)  n  =  0 


where 


0  1  o  2 

V  ij>  -  4-  -rr  -  k/ijj  =  0 

2  2  d4>  2  ^p 

H  p  K  p 


2  ^"^p  2 

k/  =  0 


<p  2  ^6  2  3  <j>  2  ^<l> 

i'  p  y  p  y 


(V^-k2^)  =  0 


„2  _  1  9  /  9  \  ^  1 
p  9p  (*^9p  )  p2 


2  2 
9  _92. 

2  ^  2 
9(|)  9z 


The  divergence  condition  Eq.(5.84)  becomes 


(5.98) 

(5.99) 

(5.100) 

(5.101) 

(5.102) 


9i|)  i|;  .  9iJ;  9iJj 

^+^  +  i^  +  _^=  0 

9p  p  p  9(fi  9z 


(5.103) 
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The  boundary  conditions  can  also  be  expressed  in  cylindrical  coordinates. 
The  results  are 


where 


2  2  2 
r  =  p  +C  and 
o 


AAA 

U  =up+u,  |J)  +  U  z 
P  <{>  Z 


u 

p 


In  +  i  z  _  (j) 

9p  p  9<|>  9  z 


iin 

P  9<j» 


9  z  ”  9p 


(5.104) 


(5.105) 


(5.106) 


(5.107) 

(5.108) 

(5.109) 


u 


z 


9'l»,  .  9i|;^ 

«%  •  _  fV  ■ 


(5.110) 


We  notice  that  the  boundary  conditions  do  not  depend  on  <j)  when  written 
in  cylindrical  coordinates.  This  being  the  case,  it  is  natural  to  seek 
a  solution  which  is  also  independent  of  (j)  .  Such  a  solution  does  in  fact 
exist. 


The  solution  in  terms  of  Hankel  transforms  which  satisfies  the  wave 
equations,  the  divergence  condition,  and  the  boundary  condition  Eq.(5.105) 
is 


«  Ve^+k z 

n(s)  =  /  B  (6)  e  J^(pB)  dB 

y  n  0 


(5.111) 


B^(e)  e  Jj^CpB)  dB 


(5.112) 


4)  (s)  =  4^  (s)  =  0 

P  Z 


(5.113) 


The  coefficients  and  B^  are  found  by  substituting  into  Eq. (5.104) 

and  (5.106).  After  some  algebra,  the  result  is  found  to  be 


\  "  “d(s,B)  ^ 


(5.114) 


-B  (  B^+k2/2) 
"  D(s,B) 


(5.115) 


where  I  da  aJ^(Ba)  ^-^2  ^'2'jn/2  ^  (Va^+5^)s) 


D(s,B)  =  2 


7  9?  2  I  2  2  2  2* 

(B^+k2/2)^  -  B  V  (3  +kp  (3  +k2) 


zs 


D(?) 


(5.116) 

(5.117) 

(5.118) 
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(5.119) 


D(?)  =  (1/2+0^ 


(Y^+0 


(1+0 


(5.120) 

(5.121) 


In  this  last  relation,  V  is  Poisson’s  ratio.  Since  for  all  materials, 

2 

0<v<l/2  ,  it  follows  that  0<Y  <  1/2  .  In  particular,  for  V=l/3  which 

2 

is  a  value  typical  of  many  materials,  y  =1/4  or  Y-1/2  . 


From  Eq. (5.108)  to  (5.110),  we  find  (including  again  the  bar 
notation) , 


Up(s) 


0 

-/ 


■J.2.,  2’ 

n  +k_ 


-I  Z  I  oV 

B  3e  ^  +  B.V3>k.  e 

h  9  2 


J^(Bp)dB  (5.122) 


u^(s)  =  0 


(5.123) 


u^(s) 


■/ 


B  V?+k  ^  e 
n  1 


1 


v;2..  2- 


,  .8  +k-  z 

+  B.3e 

9 


jQ(8p)d3 


(5.124) 


The  inverse  Laplace  transform  gives  the  time  dependent  displacements 


u^(p,z,t) 


a+i" 


a-i“ 


st 


e 


u^(s)  ds 


(5.125) 


This  completes  the  formal  solution  to  the  problem.  The  remaining 
problems  are  to  evaluate  the  required  integrals  and  to  collect  the  various 
contributions  to  form  the  complete  Green’s  tensor. 
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4 


PROCEDURES  FOR  EVALUATING  THE  RESULTANT  INTEGRAL  EXPRESSIONS 


We  will  examine  the  first  Integral  of  Eq. (5.122)  for  the  case 
F(r,t)  =  6(t-r/c)  .  The  others  have  the  same  general  form  and  can  be 
evaluated  in  a  similar  fashion.  We  define 


H(t) 


27ri 


/ 


0C+ia> 


e®*^  H(s)  ds 


a-i“ 


(5.126) 


00  VpH^r  z 

H(s)  =  -  /  3  e  ^  Jj^(ep)d3 

o 


(5.127) 


Our  intention  is  to  reverse  the  order  of  the  s  and  the  3  and  0  inte¬ 
grations,  In  order  to  justify  this  step,  it  is  necessary  to  show  that 
the  integrals  converge.  We  will  assume  for  the  present  that  z<0  and 
leave  until  later  the  limiting  case  z=0  . 

For  s  real  and  positive,  it  is  obvious  that  the  a  integra¬ 
tion  converges  and  when  in  addition  z<0  ,  it  is  equally  obvious  that 
the  6  integration  also  converges  provided  the  denominator  does  not 
vanish  for  B  positive  and  real.  It  has  been  shown,  however,  (Lamb,  1907) 
that  D(?)  has  a  simple  zero  for  a  value  of  C<-1  and  thus  in  view  of  (he 
definition  of  ?  ,  D(s,B)  does  not  vanish  for  3  in  the  range  of  inte¬ 
gration  and  s  real.  Since  both  integrals  converge  for  any  s  positive 
and  real,  the  integrated  result  defines  an  analytic  continuation  for  s 
not  real;  in  particular  along  the  line  a-i<»  to  a+i°°  where  Oi=Re(s) 
lies  to  the  right  of  any  singularities  in  the  s  plane. 


Because  the  integrals  do  converge,  we  can  reverse  the  order  of 
integration  and  study  the  properties  of  the  s  integration  first.  We 
define  a  function  G(t) 


a+i<» 


G(t) 


2Tii  J 


st 


-Ve^+k.^'izi 


+K 


r» 


ds  e  B(s,3)  e 


(5.128) 


a-ioo 
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then 


00 


(5.129) 


OC  00  2 

H(t)  =  -  A  (■  dB  3^  J  (Bp)  f  J  (3a)  G(t) 


where 


B(s,B)  = 


D(s,B) 


(5.130) 


For  notatlonal  convenience,  we  will  change  the  variable  of  integration 
from  s  to  the  real  variable  y  where  s=a+iy  .  As  we  will  show,  all 
the  singularities  lie  along  the  imaginary  s  axis  or  equivalently,  the 
real  y  axis  and  consequently  we  can  set  a=0  with  the  understanding 
that  the  integration  contour  passes  below  all  the  singularities  in  the 
y  plane.  In  terns  of  y  ,  the  function  G(t)  becomes 


eU)  -  -i/  dy 
••  00 


(5.131) 


The  sji.gularities  in  y  re' ult  from  the  branch  cuts  of  the  two  radials, 

the  simple  zero  of  D(y,B)  at  5<-i  and  the  asnmptotic  behavior  of 

D(C)  for  large  ?  .  In  the  latter  case,  it  happens  that  D(y,P)  is 
2 

proportional  to  y^  and  thus  the  integrand  has  a  second  order  pole  at 
y“0.  To  show  this,  wa  expand  about  y“0  (which  in  view  of  Eq. (5.120) 
is  large  X,  ) 
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D  =  2 


3  - 


q2  2 

iJL_ 


_ 


2c. 


0  1/2  „  ^ 


1/2 


(5.132) 


=  2 


2  2 

-  3^ 

"t 


L  -4^  -  -4t)^  O(y') 


_  -JL 
23^c; 


(5.133) 


,2  r 


1-Y 


+  O(y^) 


(5.134) 


As  mentioned  earlier,  D(5)  has  a  simple  pole  at  a  value  of  real  !^<-l  . 

/  Sc  \2 

Since  C,  =  >  this  will  result  in  a  pair  of  simple  poles  in  the  y 

plane  at 


2 


=  -1-c 


6>0  ,  real 


(5.135) 


3c. 


or 


y  =  ± 


t/I+5 


(5.136) 


The  equation  which  determines  6  is 

(1/2+5)^-  (1+6)  V6(1+6-Y^)  =  0  (5.137) 

and  thus  6  is  a  function  only  of  Y  (where  Y=Cj./c^)  which  in  turn  depends 
only  on  Poisson’s  ratio. 


The  remaining  singularities  are  the  branch  cuts  at  y=+3c^  and 
+BCj^  ,  We  chose  the  phase  of  the  two  radicals  by  requiring  these  to  be 
real  when  B  and  s  are  real.  Then 
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(5.138) 


(5.139) 


where  the  0^  are  defined  in  Figure  25.  Also  Indicated  are  the  branch 
cuts,  the  poles  and  the  integration  contour. 


Let  us  first  consider  the  limit  of  large  y  in  the  lower  half  plane. 
By  continuing  in  y  between  the  cuts ,  (Figure  26) ,  we  find 


where 


Lim 

lyh* 


6^  =  02  2Tr-  l0l 
03  =  0^^-161 
0  «  -|0| 


Thus 


Lim 

lyk” 


Lim 

iy(-)-co 


In  this  same  limit 


(5.140) 

(5.141) 

(5.142) 


(5.143) 


(5.144) 


(5.145) 
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so  that  the  Integrand  becomes 


We  see  than  that  if  T  <  0  where 


(5.14( 


T  =  t-  -  ^  (5.14 

then  the  Integrand  will  vanish  everywhere  on  the  semicircle  as  • 

We  can  thus  close  the  Integration  contour  by  including  the  semicircle 
since  its  contribution  is  zero.  Since  there  are  no  singularities  in  the 
lower  half  plane,  by  Cauchy’s  theorem  we  find  that  G(t)  vanishes  for 
T  <  0  .  Thus 


where 


(5.141 


(5.14! 


In  the  upper  half  plane  in  the  limit  of  large  |y|  ,  we  find 


and  thus 


Lira 

lyh- 


jL  =:1JL 


Llm 

|y|->eo 


r* 

2 


=  ::i2L 


(5.151) 


(5.152) 


The  integrand  becomes 


-21c/ 


(5.153) 


and  since 


-  ~  Vo^+E^'  +  =  T  + 


(5.154) 


the  integrand  vanishes  on  the  upper  semicircle  for  171*^“  whenever 
T  >  0. 


We  can  now  deform  the  integration  contour  into  the  upper  half¬ 
plane  as  indicated  in  Figure  27.  From  this  it  follows  that 

00  -c6  c6 

Jdy  ■‘‘Z 


(5.155) 


«  00  —  00 


-CtB  -c^3  «>  “ 

+ z  +  z  z 

Ct^  -C^.3  O 


c^B  Cj.B 


As  we  have  just  seen,  the  integral  over  the  semicircle  at  lyl-»-“  vanishes. 


f 


(5.156) 
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the  two 


p 


o  » 


Integrals  around  the  bran'''.!  poir:^-'^  also  vanish. 


•f  •<f  -0 


(5.157) 


The  remaining  contributions  are  nonzero.  The  contribution  from  the 
second  order  pole  at  y=0  is 


_  ♦  liini  d  I  /  \ 


(5.158) 


where  I(y)  is  the  integrand  of  G(t)  , 


i(y)  ■=  2?  ^ 


iy  [ t-1/2  I 

1-1  -  c/ 


D(y,B) 


(5.159) 


It  is  easily  determined  that 


f‘-Xr-  (t-i  VA?) 

{  3(1-7“)  '  ' 


-3h 


(5.160) 


Next,  consider  the  simple  poles  at  y^j  =  +  - -  •  Since  D(y,3)  has 

simple  zeros  at  these  points,  we  can  write 


D(y,6)  =  D(yg,3)  +  (y-Yg)  D'(yg,B)  + 


(5.161) 
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where  D(yg,3)  -  0.  The  residue  theorem  gives 


r- 

^  +  f  =  i  e  VB^  - 


2 


C2B  -c^B 
/1+6  v^+6 


iyj  ( t-  i  V7I?’)  -iy .  ( t-  i  V?^) 


D'Cyj) 


D'C-yj) 


I - ^ 

, _ ,  r  I _ • 

°  'fSs  “ _  sin  [yj  (t-  fVaV) 

B'Cyj) 


Here 


r-v 


(yg)  = 


2B' 


(1+26)  ‘ 


3±6£zWl  -4^2(1.62, 
1+6 


To  simplify  notation,  we  will  define 


e(y2)  _ _ 


^«V(1-62) 


(5.162) 


(5.163) 


(5.164) 


(5.165) 
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then 


^6  "^6 


1+6 


(5.166) 


The  integrals  along  the  right  hand  cut  can  be  combined  to  give 


c  3 


00  CO  P  •y 

f  f  i(y-io)  -  i(y+io) 


c^3  c^B 


0  2 

2 


=  real 


(5.167) 


This  must  be  separated  into  the  regions  c^S<y<c^3  and  y^c^B  . 
In  the  region  c^6<y<Cj^B  >  the  radicals  have  the  values 


(5.168) 


2 

Z_ 

2 


-  B 


2  X 


I  “i  above  cut 
I  i  below  cut 


(5.169) 


and  thus 


D(y,B)  =  2 


^  2c.  c/  c^" 


(5.170) 


with  the  +  sign  corresponding  to  y  above  or  below  the  cut.  The  resulting 
integral  is  obtained  by  substituting  into  Eq.(5.167).  If  we  include  also 
the  contribution  from  the  left  hand  cut,  we  obtain 


_ _  ^ 

C,e  -c,B  CjB  V  -4  -  b'  '(b'-  fs) 

r  r  ,  /•  C.  '  2c  / 


/„■"  /  ^  T~ 

'tB  -'tB  CjB  /  2_ 


)  -b^(b^-4)(4-b^ 


X  e  sin 


[y(t-  i 


Next,  for  Y^c^B  >  we  have  in  addition  to  Eq. (5.169) 


’  /  2 


-i  above  cut 


2*  2  ^^ji  below  cut 


D(y,B)  =  2  B 


(b2.^)  +bV(4-6')(4-b') 


Note  that  D(y,B)  is  cut  only  between  c^B<y<c^B  .  The  combined 
integrals  give 


.  ^  Ulj  sin[y(t-^^^)_ 


CnB  “C„B  C„B  (  ^2  v^  \  .  «2  HV^  „2  \/v^  .2  \ 


The  function  G(t)  is  obtained  by  adding  these  contributions. 
Defining  a  new  integration  variable 

w  «  y/c^B  (5.175) 


we  obtain 


G(t)  =  e(T)  -4 


T  e“^i^l  2 

- ECy'") 

1-Y 


(5.176) 

-Vi-yV  Blzl 

e  sin(Tw) 


-  —  J  dw  ^ Vy^w^-IB  I  z  I  _ 

(l-^)+  V(w\^-1)  (w^-1) 

where  , 

T  =  c^B  (t-  i  V?+?  )  >  0 


(5.177) 


The  inequality  follows  froiu  T  >  0.  The  integrands  of  the  above  integrals 
are  well  defined  in  the  range  of  integration.  H(t)  is  now  found  by  sub¬ 
stituting  into  Eq. (5.129).  Because  of  the  limits  imposed  by  0(T)  ,  the 
range  of  integration  for  the  a  integral  becomes 
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(5.178) 


and  consequently  this  Integral  always  converges. 

The  8  integration  of  the  first  three  terms  of  Eq. (5.176) 
also  converges  because  of  the  decaying  exponentials.  To  show 
that  the  last  term  results  in  a  convergent  integral,  we  need  an  esti¬ 
mate  of  the  asymptotic  behvaior  of  the  last  term  of  Eq.  (5.176).  The 
denominator  is  a  monotonically  increasing  function  of  w  and  can  be 
written 


(l-  +  V(w^^-l)(w^-l)  =  ^ 

w^  2 

=  ^  F(w  )  (5.179) 

2  2 
where  F(w  )  is  also  a  monotonic  function  of  w  (this  is  not  ai> 

obvious  conclusion).  To  obtain  a  bound  on  the  integral,  we  can 
2 

replace  F(w  )  by  its  largest  value,  namely  unity.  Also  in  the  numera¬ 
tor,  we  can  replace  the  radical  by  its  large  w  limit  and  the  cosine 
by  unit,  thus 


00 


1/Y 


'\/w^-l  cos  Vw^Y^-1  gizl  si 


sin  Tw 


Ii-^)^iV-i)  (w^-i) 


(5. ISO) 


continued. . . 
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Tr  'i''  'G^(;T-‘T55<^£^^i>' 


CO 

f 

.  w  sin  TW 
dvj  ^  = 

cc 

‘■f 

1/Y 

w 

4 

1/Y 

dw 


siix  Tw 


w 


(5.180) 


Integration  fay  parts  shows  that  this  last  integral  is 


/ 


dw 


1/Y 


sin  Tw 


w 


4y^  .3y^ 

cos  t/y  +  ~2  ^/Y  ■*■ 


(5.181) 


Th'iis  the  least  convergent  of  the  3  integrations  is  bous'^ed  by 


i  - 


/ 


d3  3  Jj^(3p)  J3^(3a)  X 


(5.182) 


which  is  a  convergent  integral  [see  W.N.  Bailey  (1936)]. 

We  have  now  shown  that  the  required  integrals  are  convergent.  They 
can  be  evaluated  by  a  combination  of  analytic  and  numerical  techniques. 

The  limiting  cases  of  |z|  =  0  and  ^=0  can  be  studied  by  consideration 
of  the  3  and  a  integrals  in  Eq. (5.176)  as  these  limits  are  approached. 
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SECTION  VI 

FRACTURE  ASSOCIATED  WITH  WATER  DROP  COLLISIONS 


The  critical  fracture  problems  identified  as  the  major  contributors 
to  the  degradation  of  zinc  sulfide  infrared-transparent  windows  are  frac¬ 
ture  initiation  due  to  the  transient  stresses  generated  by  the  water  drop 
impact  and  crack  propagation  due  to  hydraulic  penetration  of  the  cracks 
which  exist  after  the  initial  water  drop  collision.  The  initial  fractures 
in  polycrystalline  zinc  sulfide  display  a  high  level  of  insensitivity  to 
subsequent  stress  wave  interactions  as  observed  microscopically  for  inter¬ 
secting  ring  fracture  formations.  At  a  limited  number  of  sites,  crack  bi¬ 
furcation  is  noted  where  one  fracture  surface  propagated  in  the  direction 
of  the  critical  stress  associated  with  the  first  drop  impact  and  the  second 
fracture  surface  propagated  in  conformity  with  1;he  stress  states  generated 
by  the  second  drop  in  the  absence  of  any  initia..  fractures.  While  limited 
crack  growth  can  be  assigned  to  the  stress  wave  Interaction  mechanism,  it 
is  conjectured  that  crack  growth  by  hydraulic  penetration  is  the  dominant 
damage  mechanism  up  to  the  time  at  which  a  significant  population  of  erosion 
pits  exists  on  the  specimen's  surface.  At  this  point,  material  removal  can 
be  envisioned  to  be  due  to  both  lateral  outflow  jetting  within  erosion  pits 
and  hydraulic  penetration. 

The  need  for  a  more  precise  analysis  of  these  loading  conditions  and 
the  associated  fracture  problems  was  pointed  out  in  Section  Ill.S.b 
to  remove  the  uncertainties  in  the  response  of  the  target  to  water  drop  im¬ 
pact  collisions  have  been  initiated  and  are  described.  Additional  work  is 
required  before  the  results  from  these  analyses  will  become  available,  but 
the  general  approach  and  the  type  of  problem  which  can  be  realistically 
addressed  is  clearly  Indicated. 


Ill 
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CRACK/ STRESS  WAVE  INTERACTION 


A  Green's  tensor  formulation  will  be  developed  for  the  propagation 
of  a  surface  crack  ijito  an  elastic  half-space  due  to  the  stress  waves  gen¬ 
erated  in  the  half-space  by  a  water  drop  collision.  Eventually  the  propaga¬ 
tion  of  stress  waves  in  a  region  containing  many  cracks  will  be  considered. 


Our  starting  point  is  Eq.  (5.46).  We  evaluate  the  surface  integral 
over  the  combined  surfaces  of  the  half-space  and  the  crack.  Again  assuming 
no  sources  and  homogeneous  initial  conditions,  we  find, 


where  is  the  same  Green's  function  used  previously.  With  the 

exception  of  the  crack  opening,  the  first  integral  is  identical  to 
Eq.(5.49),  Clearly,  since  the  crack  opening  is  small,  Eq,(5.49)  yields 


da  t.  ..  *G. . 
X  (n)i  ij 


(6.2) 


except  for  a  small  region  close  to  the  crack  opening.  Since  does 

not  necessarily  vanish  on  the  crack  boundary,  we  obtain  an  inhomogeneous 

integral  equation  for  u^(5,t)  when  5  is  evaluated  on  the  crack.  We 

could,  of  course,  simply  require  T,  ,  . ,  to  vanish  on  all  surfaces  thereby 

i  J 

reducing  the  problem  formally  to  Eq.(5.49)  with  a  new  Green's  tensor.  This 
is  not  done,  however,  because  the  resulting  problem  for  the  new  Green's 
tensor  would  be  impractical  to  solve.  The  idea  in  this  approach  is  to  make 
the  Green's  tensor  problem  as  complete  as  possible  without  rendering  it  in¬ 
soluble.  The  stress  wave  propagation  problem  in  Section  5  is  complicated 
enough  without  introducing  the  crack.  The  remaining  part  of  the  problem  is 
then  treated,  for  this  case,  by  solving  an  integral  equation.  While  inte¬ 
gral  equations  are  not  in  general  easy  to  solve  in  closed  form,  they  are 
often  amendable  to  approximation  methods. 


The  surface  stress  on  the  crack  boundary  will  be  determined 

by  the  details  of  a  particular  situation.  The  simplest  case  occurs  when 
the  crack  is  air-filled  since  then  the  required  boundary  conditions  are 
(In  Cartesian  coordinates,  the  physical  stress  matrix 


ij 


and  the  stress  tensor  t^^  are  identical.)  More  generally,  if  the  crack 
is  filled  with  water,  then  ^(n)i  equal  the  hydrodynamic  stress. 


In  writing  Eq.(6.1),  no  mention  was  made  of  the  fact  that  the  crack 
is  not  stationary.  We  must  therefore  determine  the  exact  nature  of  the 
integral  expressions  in  this  equation.  Recall  that  during  the  derivation 
of  the  Green's  theorem,  it  was  noted  that  the  divergence  theorem  is  valid 
for  any  closed  surface.  It  follows  chat  we  can  always  choose  the  surface 
at  any  instant  of  time  to  coincide  with  combined  surfaces  of  the  half¬ 
space,  S'  ,  and  the  instantaneous  surface  of  the  crack. 


The  time  dependence  of  the  crack  morphology  only  enters  into  the  first 
integral  of  Eq.(6.1)  through  its  dependence  on  S*  and  thus  this  term  will 
only  be  affected  when  evaluated  close  to  the  crack  opening.  Close  to  the 
crack  means  within  a  distance  on  the  order  of  the  width  of  the  crack.  The 
second  term  is  more  complicated,  however,  due  to  the  time  dependence  of  the 
integrand  introduced  by  the  convolution. 


Consider  the  first  of  these  integrals. 


=  S  da  t ,  . .  *  G. . 

X  (n)i  n 


(n)i  ij 
t 


C  0  ^ 


(6.3) 


At  time  t=0  ,  t^^^(x,t-T)  =  t^^^(x,t)  which  is  the  stress  on  the 

crack  surface  at  time  t  .  For  t>0  however,  t,  is  evaluated 

(n)i 

in  the  interior  of  the  material  since  at  earlier  times  t-x  ,  the 
crack  is  smaller  and  consequently  the  points  specified  by  x(t)  ,  the 
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integration  surface  at  time  t  ,  lie  not  on  the  crack  surface  of  time 
t-T  ,  but  in  the  interior  of  the  material.  Thus  t^^j^(x,t-T)  is  not 
given  by  the  boundary  values  on  the  crack  surface,  but  must  be  determined 
from  the  dynamic  solution  at  each  instant  of  time  as  the  solution 
progresses.  It  is  clear  that  the  required  information  is  available  at 
any  time  t  to  proceed  with  the  calculation,  although  the  need  to  re¬ 
compute  the  integrand  at  each  instant  would  be  tedious.  Fortunately, 
we  can  simplify  this  result  considerably  when  the  crack  width  is  small. 
The  integration  surface  c(t)  is  independent  of  x  and  we  can 
interchange  the  order  of  integration  (dropping  the  tensor  indices). 


=  ^  da  /  dx  t,  ^(x,t-x)  G(x) 
C  o 

t 

=  /  dx  ^  da  t.  , (x,t-x)  G(x) 


=  /  dx[^da'^ 


t^^j(x  ,t-x)  G(x  ,x) 


^ da  t^^jCx  ,t-x)  G(x  ,x) 


(6.4) 


where  (+)  and  (-)  denote  the  sides  of  the  crack  as  indicated  in 
Figure  28.  Since  the  crack  width  is  assumed  small,  at  each  x 


n_|_  =  -n_ 


(6.5) 


G(x  ,x)  =  G(x  ,x)  . 


(6.6) 


The  latter  follows  from  the  fact  that  the  Green’s  tensor  is  the 
solution  to  the  problem  without  the  crack  and  is  thus  continuous 
everywhere  in  the  half-space. 
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Figure  28.  Schematic  of  the  Variation  of  the  Crack 
Boundary  with  Time. 


115 


Thus 


■  ^ik  ("""’t-T)]  G^j(t)  (6.7) 


Using  the  continuity  of  t^^^  and  ,  we  now  suppose  that  the 
surface  of  integration  can  be  collapsed  to  the  surface  of  the  crack  at 
time  t-T  and  a  line  extending  from  the  location  of  the  crack  tip 
at  time  t-x  to  the  location  of  the  crack  tip  at  time  t  .  The  latter 
however  vanishes  by  the  continuity  of  t^^^  and  .  These  concepts 

are  pictorially  represented  in  Figure  29. 

The  integral  ^da^  inEq.(6.7)  thus  reduces  to  ^da^(t-T) 
which  denotes  the  integral  over  the  crack  surface  at  time  t-x 

and  consequently  =  ^'ik^^^^  ~  '"ik^^  ^  need  only  be  evaluated  on 

the  crack  surface.  This  is  a  considerable  simplification  since  only 
the  boundary  values  of  t^j^  are  then  involved  in  the  calculation.  A 
similar  expression  holds  for  the  second  integral  in  Eq.(6.1). 

^2j“  /  ^da'*'(t-x)  |^u^(x'^,t-x)  -  u^(x“,t-x)j  n^^  (6.8) 


Along  the  crack  surface,  u^(x  )  and  u^(x  )  generally  have  opposite 
signs  so  the  integrand  does  not  vanish.  Thus  if  we  know  the  Green's 
function,  we  obtain  an  integral  equation  for  u^(5,t)  which  depends 
only  on  the  values  of  ^(n)i  crack  faces. 


The  approximation  of  small  crack  widths  as  used  above  can  be  examined 
by  using  the  two-dimensional,  plane  strain  solution  for  an  edge  crack  in  a 
semi-infinite  plate  subjected  to  mode  I  deformations.  If  r  and  G  are 
measured  from  the  instantaneous  crack  tip,  the  displacements  of  the  material 
point  at  the  location  (r,9)  are 


"x  -  2(l+v)  cos  I  [l-2vfslti^  I] 

“y  '  'i'Ri,  f  I] 


(6.9a) 


(6.9b) 
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collapses  to 


Figure  29.  Illustration  of  the  Collapse  of  the  Integration 

Contour  to  Include  just  the  Boundary  of  the  Crack 


The  location  of  the  crack  surface  at  time  t  ,  is  determined  by  setting 
r=Jl  and  0=H-Tr  in  Eq.(6,9).  (Refer  to  Figure  30). 


=  0 


X 

.  t 


K, 


(6.10) 


With  regard  to  the  crack  at  time  t-T  ,  the  displacements  of  the  points 
located  at  the  position  of  the  integration  surface  at  time  t  are  found 
by  substituting  into  Eq.(6.9)  values  of  r'  and  0’  given  by 


r’  =  Vu^^  +  (n-!l)^ 

e-  - 


(6.11) 


It  will  be  useful  at  this  point  to  introduce  some  numbers  into  these 
expressions.  The  stress  intensity  factor  for  this  mode  I  stress  condition 
is  =  1.12  o/ira  where  a  is  the  length  of  the  crack  (Paris  and  Sih, 
1965).  We  chose  as  a  representative  value  of  stress  the  flexural  strength 
of  ZnS  ,  0^=100.8  MPa  .  Young’s  modulus  for  this  material  is  E=74.1  GPa 
and  Poisson’s  ratio  is  0.3.  Thus  Eq.(6.10)  can  be  evaluated  to  yield 


u’  =  +  3.9x10  ^V'a’ 

y  - 


(6.12) 


where  a’  is  the  length  of  the  crack  at  time  t  .  Then 


=  +  3.9x10' 
n-i  - 


-3  ✓a^ 
n-£ 


(6.13) 


-3 

and  except  for  n-X.  <  0(10  )  , 


u 

JL. 


«  1 


(6.14) 
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Figure  30.  Coordinate  System  Used  in  Argument  Justifying 
the  Collapse  of  the  Integration  Contour  to 
Include  just  the  Crack  Boundary, 
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The  integrand  in  Eq.(6.8)  is  proportional  to  Au=u^(x^,t-T)-u  (x”,t-T). 
For  0'<Tr/2  ,  Eq.(6.9b)  is  evaluated  as  before  with  the  assumption  that  6 
is  small,  so  using  the  previous  material  properties  and  formulae. 


5 . 6x10 


|u’^+(n-A)^j  •  1/2 


(*)  [2(l-v)  -l] 


==  3.9x10 


-6  ^  a*a**£ 
’  n-il 


(6.15) 


except  for  n-il  1  0(10“^)  .  For  |0'|>ir/2  , 


3l»^.+  1-i 

(■^)  +  -  J  =  ±i 

(6.16) 

1/4 

and 

Au  =  5.6x10“^ 

|u^^+(n-£)^  I  •  2(l-v)  ==  7.8x10“^ 

(6.17) 

Since 

^ijk  continuous  at  the  crack,  the  integration  around  the 

crack 

will  be  approximately  given  by 

"k  ^ijk  . 

1  d£  Au 

'o 

(6.18) 

Breaking  the  integral  into  the  two  regions,  il<n(=a’-a")  which 
corresponds  to  e’<Tr/2  and  n<il<a'  (B^tt/2)  ,  we  find 


d£  Au  =  6.2x10"^  n[a'(a'-Ti)J 
=  /  dll  Au  =  5. 3x10"^ (a '-q)^ 


n  1/2 


n 


(6.19a) 

(6.19b) 


and  thus  except  at  a’-n  <  0(10~^)  where  vanishes  because  the 
crack  vanishes. 


<  -3 

—  10  «  1 


(6.20) 
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The  region  very  close  to  the  crack  tip  where  <  0(10  )  must 

be  re-examined  but  the  general  conclusion  will  not  be  changed.  Consequently 
the  surface  of  integration  can  safely  be  collapsed  to  enclose  just  the  crack 
at  time  t-t  without  introducing  appreciable  errors. 

We  now  have  a  procedure  then  for  calculating  the  stress  and  displacement 
distributions  from  the  Green's  function  and  the  boundary  stresses  provided  we 
know  the  morphology  of  the  crack.  The  remaining  steps  are  then  to  determine 
the  Green's  tensor  and  the  relation  between  the  dynamic  field  variables  and 
the  growth  of  the  crack. 

2.  HYDRAULIC  PENETRATION  OF  SURFACE  CRACKS 

Hydraulic  penetration  of  existing  surface  cracks  has  been  identified 
as  a  significant  crack  extension  mechanism  (Adler*  1979a)  with  characteris¬ 
tics  distinct  from  those  of  the  stress  wave  extension  mechanism.  One  prin¬ 
cipal  difference  is  in  the  direction  of  crack  growth  induced  by  these  two 
mechanisms.  It  has  been  noticed  for  example  that  hydraulically  activated 
cracks  tend  to  undermine  sections  of  the  surface  thus  weakening  large  pieces 
of  the  surface  which  are  subsequently  removed.  This  behavior  is  strongly 
dependent  on  the  characteristics  of  the  material  and  thus  the  first  step 
toward  understanding  and  describing  this  process  is  to  determine  the  rela¬ 
tionship  between  the  angle  of  crack  extension  and  the  other  parameters  of 
the  problem  such  as  the  crack  depth,  inclination  angle,  and  material  proper¬ 
ties.  In  general  the  analysis  of  this  problem  is  quite  complicated  because 
of  the  hydrodynamics  of  the  water  flow  in  the  crack  and  the  interaction  of 
the  crack  with  impact  induced  stress  waves.  A  much  simpler,  but  relevant 
problem  is  analyzed  here:  namely,  a  two-dimensional  static,  pressurized 
surface  crack  of  finite  length,  L  ,  inclined  to  the  surface  at  an  arbi¬ 
trary  angle.  The  immediate  goals  of  this  study  are  to  determine  the 
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stress  distribution  surrounding  the  crack  and,  using  one  or  more 
fracture  model,  to  determine  the  initial  angle  of  crack  extension. 

Of  particular  interest  is  the  effect  of  the  proximity  of  the  surface 
on  the  crack  extension  angle. 


The  stress  distribution  in  two  dimensions  is  most  conveniently 
stiidied  by  introducing  complex  variables.  As  indicated  in  Figure  31,  the 
material  occupies  the  upper  half  plane.  The  location  of  a  point  in  the 
material  is  specified  by  a  complex  number  z  given  by 

z  =  X  +  iy  .  (6.21) 

Complex  conjugates  will  be  denoted  by  a  bar,  e.g., 

z  =  X  -  iy  .  (6,22) 

The  equation  of  equilibrium  becomes  in  complex  notation 

9.\(z,z)  ^  Q  ^g^23) 

9  z3  z 


where  x(z,z)  is  the  Airy  stress  function.  The  general  solution  of  this 
equation  is  given  in  terms  of  two  arbitrary  functions  ())(z)  and  i|>(z)  which 
are  analytic  in  the  region  interior  to  the  material  boundary.  In  terms  of 
these,  the  Airy  stress  function  is 


x(z»z) 


z 

i|<(z)dz 


(6.24) 
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The  stresses  are  given  by 


a  +0=4  Re$(z)  (6.25) 

X  y 

Oy  -  o^  +  2iT^  =  2  [s$'(z)  +  'l'(z)]  (6.26) 

where 

$(z)=«i‘(z)  (6.27) 

•l'(z)  5  4»»(2)  (6.28) 

The  boundary  condition  which  determines  these  functions  is 

N(z)-iT(z)  =  $(z)  +  $(z)  -  e^^“[i$’(z)  +  'l'(z)]  (6.29) 

A 

where  N(z)  is  the  projection  of  the  surface  stress  along  n  ,  the 
outward  unit  normal  vector,  and  T(z)  is  the  projection  along  the 
surface  tangent  pointing  to  the  left  of  n  when  viewed  from  inside  the 

A 

material.  The  angle  between  n  and  the  x-  axis  is  denoted  by  a  . 

Since  we  are  considering  hydrostatic  forces, 

N(z)  =  -P(z),  T(z)  =  0  (6.30) 

where  P(z)  is  the  pressure  which  in  general  will  depend  on  position. 

We  are  thus  confronted  with  the  problem  of  finding  two  analytic 
functions  which  satisfy  the  boundary  condition  Eq.(6.29).  Rather  than 
approach  this  problem  directly  however,  we  will  instead  utilize  a  method 
introduced  by  Muskhelishvili  (3,953)  who  showed  that  by  extending  the 
definition  of  <{)(z)  into  the  lower  half-plane,  the  function  ’l'(z)  can 
be  eliminated. 


,,,'.v*  =s-\'  ,  »-5.^y^ 


We  define  a  function  $(z)  analytic  in  the  lower  half-plane  by 


4>(z)  =  -  5(z)  -  z4>*(z)  -  ^(z) 


zeS 


(6.31) 


Here  the  bar  notation  signifies 


$(z)  =  $(z) 


We  now  set  z-»z  and  take  the  complex  conjugate.  Rearranging  gives 


(6.32) 


4<(z)  =  -  'J>(z)  -  $(z)  -  z$'  (z) 


Z€S 


(6.33) 


where  now  z  is  defined  in  the  upper  half-plane.  Thus  by  extending 
the  definition  of  ^(z)  into  the  lower  half-plane,  the  function  '^(z) 
is  determined  in  the  upper  half-plane.  Introducing  this  result  into 
Eq.  (6.29)  we  obtain 

2ia 


N(z)-iT(z)  =  'I>(z)+  $(z)  -  e  £(z-z)  $'(z)  -  $(z)  -  i(z)]  (6.34) 

Consider  now  a  region  of  the  real  axis  for  which  N  and  T  vanish. 

TT 

Since  a=  -  —  along  this  material  surface,  we  have 


0  =  $(z)  +  $(z)  +  (z-z)$'(z)  -  4>(z)  -  $(z) 


(6.35) 


0  =  4>(z)  -  $(z)  +  (z-z)$’(z) 


(6.36) 


Since  z-z=0  on  the  real  axis,  provided  <J>'(z)  is  not  too  singular 
we  obtain 


^(z)  =  $(z)  at  z  =  z 


(6.37) 
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This  means  that  $(z)  as  extended  by  Eq.  (6.31)  is  continuous  across  the 
real  axis  and  it  follows  that  3>(z)  is  analytic  everywhere  in  the  plane, 
except  perhaps  on  portions  of  the  boundary  where  N  and/or  T  are  non-zero. 

The  asymptotic  behavior  of  $(z)  has  been  investigated  by 
Muskhellshvili  (1953)  who  showed  that 

K  =  3  -  4v 

where  T  is  related  to  the  limit  of  the  stresses  as  ,  X+iY  is 

the  resultant  of  the  surface  stresses  over  the  finite  boundaries,  and  v  is 

Poisson’s  ratio.  If  the  stresses  vanish  as  z-»-“  ,  r=0  and  if  the  surface 

stresses  are  self-equilibrating,  then  X+iY=0  .  The  function  (z)  is  analy- 

-2  ^ 

tic  everywhere  except  at  the  origin  and  vanishes  as  z  as  .  In  order 

to  simplify  the  calculations  as  much  as  possible,  we  will  suppose  that  the 
pressure  is  self-equilibrating  and  further  that  the  pressure  vanishes  every¬ 
where  on  the  real  axis.  Clearly  by  superposition,  other  more  general  solu¬ 
tions  can  be  found  by  the  addition  of  solutions  for  the  cases  in  which  this 
pressure  is  not  zero. 


Along  the  crack  (see  Figure  31)  ,  we  have 


2i« 

e 


21 


'2—2/ 


e 


iXir 

-e 


(6.39) 


with  the  plus  or  minus  sign  depending  on  which  side  of  the  crack  we  are 
considering.  Thus  the  boundary  condition  becomes 

$(z)  (l-e^""^)  +  -  e^’^^  5(z)  +  e^’'^(z-z)$’(z)  +  P(z)  =  0  (6.40) 
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The  boundary  condition  is  not  in  a  convenient  form  because 
the  crack  shape  does  not  allow  a  simple  power  series  expansion  for  $(z)  . 
To  remedy  this  problem,  we  Introduce  a  conformal  mapping  which  maps  the 
crack  into  a  unit  semicircle  (Figure  31)  .  The  mapping  which  accomplishes 
this  is  (Bowie,  1973), 

z  =  w(?)  =  A  Y  (5-1)^  (C+1)^'^  (6.41) 


where 


A  = 


(6.42) 


The  angle  0^  defined  by 


0Q  =  tan 


I  2  1/2^ 

>1  I  (2X-A^) 

^  1-X  I 


(6.43) 


determines  the  location  of  the  crack  tip  on  the  unit  circle  in  the 
4  plane.  It  is  easily  shown  that  the  positive  and  negative  real  axes 
in  the  z  plane  map  into  the  corresponding  real  axes  for  | C 1 >1  in 
the  ?  plane.  The  crack  maps  into  the  unit  semicircle  in  the  ?  plane. 
Introducing  the  notation 


«’(w(C))  =  $(?) 


(z) 


w'(?) 


i0 

the  boundary  condition  becomes  (on  the  unit  circle  ?=a=e  ) 


i7rX 


iirX 


$(a)  d-e""'')  +  $(o)  -e  '  '  $(a) 


■2SP^-U — $'(a)  +p(a)  =  0  O<0<Tr 


0  +2(X-l)a+l 


(6.44) 

(6.45) 


(6.46) 


■ . ■■' .  /  '  ■  ■"■" 


fS. 

I? 


Sw. 

* 

e>'4 

>v- 


i 

M 

M' 

I 

Ss 


The  coefficient  of  $'(a)  has  poles  at 

=  1-X+i  VaX-X^ 


(6.47) 


By  counting  singularities,  one  finds  that  v/riting  $(C)  in  the  form 

.4 


$(?)  = 


A(L)-.;. 


iZ-o^)  (^-oj- 


(6.48) 


defines  a  function  G(?)  which  is  analytic  external  to  the  unit  circle  and 
has  no  severe  singularities  on  the  unit  circle.  It  follows  that  G(^)  can 
be  represented  by  a  power  series  of  the  form 


G(0  =52  Aj^  C 


-k 


k*2 


The  inversion  of  this  expression  gives 


J  =  1 


2x 


/*■"  i0,  inG 

de  G(a  )  e 

0 


n 


n  £  1 
n  >  2 


(6.49) 


(6.50) 


Thic  integral  can  be  separated  into  two  regions 


•'n  ■  i,/ 


iG,  inG  1 
dG  G(e  )  e  +  2^ 


'•2tt 


-J  d9  G(e  )  a 


In  the  firit  of  these,  we  introduce  Eq,  (6.49)  to  find 
ir 


JO  d®  '  * 

de  G(a  J  e  =  ^  + 

k=2 


1.  V  A  lS=lfll=l\ 


n-k 


kj4n 


(6.51) 


(6.52) 


In  the  second  of  these  integrals,  G  ,  which  is  evaluated  in  the  lower 
half-plane,  is  given  by  the  boundary  condition  Eq.(6.46).  If  we  define 
a=2'?T-0  ,  we  obtain 


i/ 


o/  1  /*  j  n/  “ina 

dS  G(e  )  e  =  —  /  da  G(e  )  e 


1  r^ja  nt 

-  J  d9  G(e  )  e 


(6.53) 


From  the  boundary  condition,  it  follows  after  some  algebra  (recall 
G(a)  =  G(a)  ,  a=  ^  ,etc.) 


G(a)=G(a)  F^(a)  +  G(0)  F2(a)  +  G'(o)  ?^ia)  +  P(a)  F^(a) 


0  <  0  <  ir 


where 


F^(o)  -  (j^) 


(6.54) 


(6.55) 


O  r  2 

F  (a)  =  (e^""^-!)  - ^  (aM7^)^+  ^  (4a-3a^-a_)  (6.56) 

(.0-0  )  (a-a  ) 


7^(0)  = 


(e^’^^-l)  (l-0^)o_^0-a^) 
(o-o  )^ 


(6.57) 


^iirX  2-  .  j  f  V 


F  (a)  =  e"o_(o-o_^)^'  (o-oj 


(6.58) 


i 

! 
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Substituting  Eq.(6,49)  into  (6.54)  and  the  result  into  (6,53)  yields 


TT  00 

i  J[  de  G(a)  e"^’"®  lA^I^(k+n)  +  I2(k-n)-kl3(k-n) 


+I^(n)  (6.59) 


where 


T  1  -i(k+n)0^  ,  10, 

Ij^(k+n)  =  7  J  d0  e  Fj^(e  ) 


(6.60) 


,  1  /•IT  .  i(k-n)0  _  /  i0, 

l2(k-n)  =  -J-  d0  e  '  F2(e  ) 


(6.61) 


TT 


l3(k-n)  =  i  J  d0 


,i(k-n)0  ^  .10, 

■3' 


(6.62) 


T  /•  ^  1  r"  JO  T>/  i6^  T7  /•  i0x  -in0 


(6.63) 


Combining  Eqs.(6.50),  (6.51),  (6.52)  and  (6.59)  yields  an  infinite  set  of 
nonhonogeneous  linear  simultaneous  equations  for  the  infinite  number  of 
complex  coefficients  •  Uniqueness  of  the  solution  is  assured  by  the 
uniqueness  of  analytic  functions  for  a  given  set  of  boundary  conditions. 


Since  the  complex  coefficients  A^  represent  an  infinite  set 
governed  by  an  infinite  set  of  nonhomogeneovs  linear  simultaneous  equa¬ 
tions,  there  is  a  degree  of  arbitrariness  in  the  selection  of  a  trunca¬ 
tion  procedure  for  determining  the  coefficients.  The  first  attempt  in¬ 
volved  solving  n  complex  equations,  equally  divided  into  the  groups 
n^2  and  n<2  ,  for  n  complex  coefficients  Aj^  ,  k=2,3,..on+l  .  It  was 
found  however  that  the  coefficients  so  determined  increased  rapidly  with 
increasing  n  and  thus  did  not  lead  to  a  convergent  expression  for  the 
function  G(5)  in  Eq.(6.49).  In  order  to  Inhibit  this  tendency. 
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additional  equations  were  Included  so  that  n'  equations  (n'>n)  were 
solved  for  n  coefficients  In  a  least  squares  sense.  For  X«  1  ,  this 
procedure  results  In  a  rapidly  converging  sequence  of  coefficients  when 
n*>2n  ,  the  accuracy  of  which  was  checked  by  comparing  with  the  coeffi¬ 
cients  determined  by  Bowie  (1973)  for  the  case  of  a  crack  normal  to  the 
surface  X=1  .  For  X  much  less  than  unity,  on  the  other  hand,  some 
practical  problems  remain  with  this  procedure  because  the  number  of  coef¬ 
ficients  required  for  convergence  appears  to  grow  as  X  decreases  and,  as 
a  consequence,  the  amount  of  computer  time  required  becomes  excessive. 

We  have  not  yet  exhaustively  studied  the  possible  variations  within  this 
class  of  algorithms  and  thus  further  work  may  yield  a  more  efficient  pro¬ 
cedure  which  will  reduce  the  amount  of  computer  time  required  to  obtain 
a  solution.  An  alternate  method  for  determining  the  stress  near  a  crack 
tip  has  been  discussed  by  Hussian,  et  al.  (1974),  for  a  related  problem. 

This  method,  which  is  based  on  the  solution  of  an  integral  equation  in 
the  limit  as  the  crack  tip  Is  approached,  will  be  discussed  briefly  belov; 
after  some  other  aspects  of  this  same  paper  which  are  relevant  to  the 
present  problem  are  discussed. 

Before  doing  so,  however,  we  will  derive  the  expressions  which  relate 
the  stress  intensity  factors  and  the  function  G(C)  .  For  fracture  Modes  I 
and  II,  the  stress  near  the  crack  tip  is  known  to  have  the  form 
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where  the  parameters  r,0  are  defined  in  Figure  33  ,  If  we  set  0=0 
and  use  Eqs. (6.25)  and  (6.26),  we  find 

K  =  Lim  I  l-Zluv  Re  <J>  (z)  1  (6.66) 
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After  some  algebra,  these  reduce  to 
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which  allows  us  to  calculate  the  stress  neai  the  crack  tip  from  G(0^) 
or  equivalently,  the  . 


The  next  problem  is  to  relate  the  stress  near  the  crack  tip  to  the 
growth  of  the  crack.  Many  models  have  been  proposed  to  describe  the  rela¬ 
tion  between  stress  and  crack  growth  but  most  are  too  limited  in  scope  to 
treat  combined  mode  cracking.  TAro  models  which  deal  specifically  with  com¬ 
bined  mode  cracking  are  those  of  Sih  (1973)  and  Hussain,  et  al.  (1974). 

In  Sih's  model,  the  crack  is  assumed  to  extend  in  the  direction  of  the  maxi¬ 
mum  potential  energy  density  where  the  potential  energy  density  P  is  given 
by 


P  =  - 


r 
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(6.69) 
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—  is  the  elastic  free  energy  in  the  limit  as  the  crack  tip 


approaches  (irK)) 
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Since  S  is  a  quadratic  function  of  the  stress  intensity  factors  at  the 
crack  tip,  the  direction  of  crack  growth  predicted  by  this  model  can 
easily  be  calculated  once  these  factors  are  known.  In  fact,  since 


(3-4v-cos9)  (l+cos6) 

Xop  1 
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the  condition  that  P  be  a  maximum  yields 


3S  .  3^S  ,  „ 
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where 
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03  =  2(1-2V)  (Kj-Kjj) 
a^  =  4(1-2V)  Kj  Kjj 


Thus  given  K^  and  ,  a  crack  extension  angle  can  be  determined.  Crack 
extension  is  supposed  to  occur  when  S  reaches  a  critical  value  charac¬ 
teristic  of  each  material. 

More  recently,  Hussain,  et  al.  (1974),  proposed  a  model  for  mixed  mode 
crack  extension.  In  some  respects,  it  is  similar  to  Sih's  model  although 
the  crack  extension  criterion  is  stated  in  terms  of  the  energy  release  rate 
given  by  the  invariant  J  integrals  rather  than  in  terms  of  the  elastic 
free  energy.  The  important  difference  however  is  that  Hussain,  et  al., 
treat  the  stress  in  the  neighborhood  of  the  crack  tip  as  chat  for  a  deflected 
crack  in  the  limit  as  the  propagation  branch  of  the  crack  shrinks  to  zero 
length.  They  found  that  the  stress  determined  by  this  limiting  procedure  is 
dependent  on  the  angle  between  the  propagation  and  main  branches  of  the 
crack.  Consequently  the  energy  release  rate  also  depends  on  this  angle  and 
they  argue  that  the  crack  will  extend  in  the  direction  of  maximum  energy  re¬ 
lease.  This  extension  criterion  is  more  difficult  to  apply  in  general  be¬ 
cause  a  more  complicated  boundary  value  problem  must  be  solved  to  determine 
the  appropriate  stress  intensity  factors. 

The  procedure  used  by  Hussain,  et  al.,  is  to  convert  the  boundary 
value  problem  for  the  stress  function  into  an  integral  equation.  This  step 
is  useful  because  an  exact  solution  can  be  found  in  the  double  limit  as 
the  propagation  branch  shrinks  to  zero  length  and  the  crack  tip  is  approached. 
This  yields  sufficient  information  to  calculate  the  stress  intensity  factors. 
It  is  felt  that  this  technique  might  also  prove  useful  for  determining  the 
stress  intensity  factors  for  the  present  problem.  If  this  is  indeed  the  case, 
then  the  convergence  problems  discussed  earlier  with  regard  to  the  series  solu¬ 
tion  would  be  circumvented.  Further  work  will  be  required  to  verify  this  sup¬ 
position.  We  emphasize  that  the  use  of  integral  equation  techniques  is  dis¬ 
tinct  from  the  concept  of  calculating  the  stress  in  the  limit  of  zero  propa¬ 
gation  branch  length.  The  deflected  crack  problem  could  also  be  solved  using 
the  series  expansion  techniques  discussed  earlier  if  the  integral  equation 
method  proved  to  be  too  cumbersome.  The  remaining  question  is  to  determine 
which, if  either,  of  these  models  is  correct  and  this  can  only  be  done  by 
comparison  with  experiment. 
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SECTION  VII 

TRANSMISSION  LOSSES  DUE  TO  CRACKS 


It  is  intuitively  obvious  that  optical  degradation  will  occur  when 
the  surface  of  an  infrared-transmitting  window  becomes  pitted  on  the  scale 
of  the  diameter  of  a  typical  rain  drop  since  then  the  transmitted  infrared 
radiation  with  a  wavelength  X  on  the  order  of  1  to  10  ym  will  be  scat¬ 
tered  in  various  directions  vis  a  vis  Snell’s  law.  This  is  particularly 
true  for  materials  such  as  zinc  sulfide  and  zinc  selenide  which  have  rather 
large  indices  of  refraction.  At  X=3,3  pm,  n  =2,3;  and  for  X=8-13  pm, 
^ZnSe”^*^*  ^  important  observation  with  regard  to  the  subsurface 

fractures  described  in  Section  4,  however,  is  that  significant  degradation 
occurs  before  there  is  a  measurable  mass  loss  from  the  sxirface.  This  in¬ 
dicates  the  importance  of  light  scattering  by  the  internal  fracture  sur¬ 
faces  which  develop  due  to  single  and  multiple  water  drop  collisions. 

It  is  easy  to  understand  the  origin  of  the  transmission  degradation 
by  considering  the  optical  properties  of  an  idealized  crack  as  indicated 
in  Figure  34  ,  The  crack  faces  are  taken  to  be  parallel  in  order  to  avoid 
unnecessary  complication  in  the  elementary  analysis  which  follows.  When¬ 
ever  electromagnetic  radiation  is  incident  on  an  interface  between  two 
regions  with  differing  indices  of  refraction,  generally  both  transmitted 
and  reflected  waves  result,  Figure  35a.  If  the  index  of  refraction  of 
the  second  medium,  n2  ,  is  less  than  that  of  the  first,  n^  ,  then  a 
critical  angle  of  incidence  exists  such  that  the  reflection  is  total 

for  angles  of  incidence  »  Figure  35b. 

This  angle  is  given  by 

sina  =  —  (/.I) 

c  n^^ 

For  the  interface  between  zinc  sulfide  and  air  and  water,  respectively, 
this  angle  has  the  values 
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Figure  ?5.  Iransinission  of  Electromagnetic  Radiation 
at  the  Interface  between  Two  '''edia. 


25.8“ 


ZnS  -  air 
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34.4“ 


ZnS  -  water 


(7.2) 


For  the  case  of  a  crack  which  has  two  interfaces,  the  situ<ition  is 
more  complicated  because  transmission  can  occur  for  angles  greater  than 
by  a  process  analogous  to  quantum  mechanical  tunneling,  provided  the 
crack  width  d  is  not  too  large.  For  the  electromagnetic  field  does  not 
vanish  instantaneously  at  the  interface  in  Figure  35b  under  conditions  of 
total  internal  reflection,  but  instead  attenuates  exponentially  according 
to 


a>a 

c 


(7.3) 


where  Z  is  the  distance  normal  to  the  interface  in  region  2  and 

2i[ 

^2=  ~  is  the  wave  number  in  region  2.  It  is  easily  shown  that  even 

though  the  field  is  nonzero  in  region  2,  no  energy  is  propagated  in  the 
Z-direction  by  this  field.  If,  however,  a  second  interface  is  placed 
close  to  the  first  in  the  region  where  the  field  is  nonzero,  then  energy 
can  flow  across  the  gap  and  some  transmission  will  occur.  The  transmis¬ 
sion  coefficient,  T  ,  defined  to  be  the  ratio  of  I he  transmitted  energy 
to  the  incident  energy,  for  the  crack  is  determined  by  solving  the  bound¬ 
ary  value  problem  for  the  geometry  shown  in  Figure  34  .  For  unpolarized 
light,  the  transmission  coefficient  is 


T  =  1/2 


(7.4) 
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where  for  a<a  , 
c 


y  =  k2ci  cosa2 


and  for  a>a 

c 


The  eubscripts  J_  and  |  j  refer  to  polarization  r-f,  the  electric  field 
perpendicular  and  parallel  to  the  plane  of  incidence.  These  expressions 
include  the  effects  of  multiple  reflections  in  the  crack.  Notice  that 
since  1?|  |  1  »  transmitted  light  will  in  general  be  partially 

polarized. 

Some  representative  evaluations  of  Eq.(7.4)  are  provided  in 

Figures  36  to  39  .  Figure  36  shows  the  dependence  of  T  on  a  for 

fixed  values  of  the  ratio  d/X  .  X  refers  to  the  wavelength  in  air. 

For  values  of  »  essentially  no  transmission  occurs  for  as 

expected.  The  rapid  fluctuations  in  T  for  a<a  are  due  to  the  con- 

c 

structive  and  destructive  interference  of  the  primary  and  reflected  waves 
emerging  from  the  crack.  For  small  values  of  d/X  ,  on  the  other  hand. 


significant  transmission  can  occur  for  a>a 


Figures  37  and  38  show 


the  dependence  of  T  on  d/X  for  fixed  angle.  Figure  38  shows  the 
results  from  Figure  37  on  an  expanded  scale.  It  is  seen  that  for  a<a^  , 
T  oscillates  indefinitely  as  a  function  of  d/X  in  response  to  the  in¬ 
terference  effects  mentioned  above  never  dropping  below  about  .55;  bat  for 
a>a^  ,  T  decays  rapidly  to  zero  as  d/X  is  increased.  Notice  also  that 
T  approaches  unity  for  all  angles  if  d/X  is  sufficiently  small  which 
would  be  the  case  if  the  crack  were  to  close  after  its  formation. 


The  results  in  Figures  36  to  38  were  calculated  assuming  the 
crack  to  be  filled  with  air.  For  comparison,  the  analogous  resu.lts  for 
water  filled  cracks  are  shown  in  Figure  39  .  For  ease  of  comparison,  the 
dependence  of  T  on  d/X  for  air  when  aj^=50'’  is  also  plotted  in 
Figure  39,  It  is  apparent  that  the  transmission  coefficient  is  enhanced 
by  the  water  which  is  expected  since  increasing  decreases  the  atten¬ 

uation  coefficient  F  ,  Eq.(7.6d).  This  shows  incidently  that  improved 
transmission  will  be  obtained  for  a  given  crack  width  by  decreasing  the 
index  of  refraction  of  the  material. 


142 


Function  of  Angle  of  Incidence 


Transmission  Coefficient 


0^=20° 


a^=50\  \  ai=40° 


0.4  0.6 


1.0  1.2 


Dimensionless  Width  Parameter,  d/> 


Figure  38.  Transmission  Across  an  Air-filled  Inclined  Gap 
in  Zinc  Sulfide  as  a  Function  of  Gap  Width 
(expanded  scale) 


To  relate  these  results  to  the  net  transmission  of  a  window,  it  is 
necessary  to  consider  the  length  of  the  crack.  This  will  eliminate  the 
apparent  anomaly  of  zero  transmission  at  aj^=90°  ,  since  then  the  projected 
c'rea  of  the  crack  is  also  zero  and  the  crack  does  not  cause  a  loss  of 
transmission. 

The  calculations  presented  here  do  not  yet  include  the  effects  of 
attenuation  in  the  material  and  reflections  at  the  window  surfaces.  Since 
these  later  effects  are  multiplicative,  the  influence  on  radiation  trans¬ 
mission  of  the  cracks  developed  during  water  drop  impacts  can  be  adequately 
addressed  without  treating  these  other  effects  in  detail. 

If  the  incident  radiation  is  normal  to  the  window  surface,  a  simple 
expression  is  obtained  for  the  time  dependent  loss  of  transmission  due  to 
crack  formation.  Denoting  the  incident  and  transmitted  energy  per  unit 
area  per  unit  time  by  and  ,  respectively,  a  transmission  coeffi¬ 

cient  T  is  defined  for  an  uncracked  window  by 
o 

'  Vi 


T^  includes  the  effects  of  reflection  at  the  window  surfaces  and  energy 
absorption  by  the  window  material.  When  cracks  are  present,  the  total 
energy  transmitted  will  be  the  sum  of  that  transmitted  by  the  uncracked 
and  cracked  regions.  Thus, 


s^A  =e  T  /a  -T  A  \  +  E  e„  A 
T  0  X  o^  o  n  n/  n  Tn  n 

where  A  is  the  total  surface  area  of  the  window,  A^  is  the  area 
o  th  ” 

projected  onto  the  window  surface  of  the  n  crack  at  time  t  ,  and 

is  the  energy  transmitted  through  this  area,  will  include 

the  attenuation  effects  of  T  in  addition  to  that  of  the  crack  and 

o 

it  follows  that 


'Tn 


e^T  T 
ion 


(7.9) 
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where  is  the  transmission  coefficient  for  the  n  crack.  Introducing 

the  definition  in  Eq.(7.9)  into  Eq.(7.8), 


=  e.T  Fa -I  A  /l-T  ] 
To  i  o  L  o  n  n\  n/J 


(7.10) 


We  have  assumed  that  cracks  do  not  overlap.  If  they  do,  the  total  crack 
attenuation  is  found  by  substituting  for  the  appropriate  products  of 

transmission  coefficients  for  the  particular  cracks  involved.  The  time- 
dependent  transmission  coefficient  due  to  crack  generation  is  determined 
from  Eq.(7.10), 


P  A 

T(t)  =  — ^  =  T  [l- 1  I  A  (l-T  )1 


(7.11) 


In  order  to  proceed  further,  we  need  an  estimate  of  the  average  crack 
inclination  angle.  Experimental  studies  show  that  the  majority  of  the 
crack  surfaces  are  inclined  to  the  surface  at  angles  greater  than  45° 
except  for  those  in  fracture  zone  2.  A  representative  angle  of  a^=55°  , 
an  angle  which  is  much  greater  than  for  either  air-  or  water-filled 

cracks  will  be  selected  for  the  purpose  of  considering  an  explicit  example. 
Referring  to  Figure  38,  we  see  that  such  cracks  will  be  almost  totally 
opaque  for  d/A>0.3.  Suppose  for  simplicity  that  all  the  cracks  are  formed 
at  the  same  angle  with  the  same  projected  area  and  are  created  at  a  con¬ 
stant  rate  of  N  cracks  per  unit  time,  then  the  transmission  coefficient 
for  the  window,  Eq.(7.11),  reduces  to 


T(t)  =  T„[l-  I  (1-T^)] 


(7.12) 


where  t  =  is  the  time  required  to  completely  obscure  the 

surface. 
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is  the  projected  area  of  a  single  crack, 

is  the  transmission  coefficient  for  a  single  crack. 

The  implications  of  this  result  are  that  no  incubation  period  is  required 

before  degradation  begins  and  that  the  transmission  loss  will  increase 

linearly  with  exposure  time.  The  expression  in  Eq.(7.12)  is  valid  only 

for  times  small  relative  to  t  . 

o 

Figure  40  shows  the  degradation  of  transmission  for  aj^=55‘’  for 

several  values  of  d/X  .  It  is  apparent  from  this  figure  that  the  value 

of  d/X  has  a  significant  effect  on  the  rate  of  degradation.  The  line 

labeled  (T  =0)  is  the  limiting  line  for  the  case  of  completely  opaque 
c 

cracks  which  is  reached  for  d/X>*3  (see  Figure  38)  ,  Thus,  for  example, 
for  a  crack  width  of  d=.2  pm  ,  the  crack  would  be  essentially  opaque  for 
wavelengths  of  X=lpm  and  nearly  transparent  for  X=20pm.  Behavior  of 
this  sort  has  been  observed  by  Hackworth  and  co-workers  (Hackrjorth  and 
Kocher,  1977;  Hackworth,  et  al.,  1979).  According  to  their  work,  the 
infrared  transmittance  of  zinc  sulfide  decreased  almost  linearly  with  in¬ 
cremental  exposure  to  the  standard  rainfield  in  the  AFML/Bell  rotating 
arm  erosion  facility  for  wavelengths  in  the  range  0.5  to  2.1  pm.  However 
they  indicate  that  there  is  a  definite  incubation  period  associated  with 
the  development  of  surface  pitting  for  the  degradation  of  transmittance  in 
the  range  of  2.5  to  25  pm.  Once  erosion  pits  are  nucleated  the  rate  of 
transmittance  loss  was  approximately  proportional  to  the  growth  rate  of 
the  total  cross-sectional  area  of  the  pits.  Hackworth  and  co-workers 
speculate  that  the  same  relations  apply  to  zinc  selenide  and  gallium 
arsenide,  however  the  rate  of  material  removal  was  too  great  for  ex¬ 
posures  to  the  standard  rainfield  to  experimentally  confirm  this 
conjecture. 

Obviously  this  model  is  oversimplified  since  many  effects  such  as 
nonparallel  crack  surfaces  and  overlapping  cracks  were  not  considered.  It 
is  expected  however  that  a  more  complete  model  which  includes  these  and  other 
effects  and  the  statistical  nature  of  the  particle  field  would  yield  a 
similar  result  for  the  transmission  coefficient. 
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SECTION  VIII 
DISCUSSION 


The  objective  of  this  effort  was  to  detemine  the  relationships 
between  water  drop  impact  damage  and  the  material  properties  of  infra- 
red~transparent  materials  in  order  to  evaluate  and  guide  the  development 
of  erosion  resistant  window  materials.  A  combined  experimental  and  analy¬ 
tical  approach  was  undertaken  to  determine  and  describe  primarily  the 
water  drop  impact  response  of  CVD  zinc  sulfide,  however  the  results  of 
this  program  are  of  more  general  utility. 

The  impact  damage  mechanisms  were  investigated  for  water  drop  impacts 
at  velocities  from  222  to  684  ms  ^  (730  to  2240  fps)  with  the  major  empha¬ 
sis  on  the  velocity  range  from  342  to  512  ms  Both  experimental  and  com¬ 
putational  procedures  were  utilized  to  relate  the  dominant  damage  modes  to 
the  material’s  surface  condition,  microstructure,  and. bulk  properties.  The 
most  productive  analysis  for  evaluating  the  bulk  response  of  the  material 
to  the  water  drop  impact  loading  involves  continuum  representations  of  the 
target  material. 

The  complete  description  of  the  impact  process  entails  the  solution 
of  a  coupled  problem  including  the  deformations  of  both  the  incident  par¬ 
ticle  and  the  target  material.  At  the  relatively  low  velocities  of  interest 
and  for  highly  deformable  particles  the  problem  can  be  uncoupled  as  a  first 
approximation.  The  effect  of  the  incident  particle  can  then  be  treated  as 
a  known  temporally  and  spatially  dependent  pressure  distribution  on  the  un¬ 
deformed  surface  of  the  target  material.  This  pressure  distribution  has 
been  determined  using  finite  difference  procedures  and  experimentally  for 
water  drop  impact  simulations  as  reviewed  by  Adler  (1979a).  The  analytical 
representation  of  the  surface  pressure  distribution  developed  by  Rosenblatt, 
et  al.  (1977)  with  minor  modifications  will  be  adequate  for  the  initial  im¬ 
plementation  of  the  computational  program  which  will  be  an  outgrowth  of  the 
completed  analysis  described  In  Section  V,  The  reasons  for  undertaking  this 
analysis  are  outlined  in  Section  V.l.  While  the  scope  of  this  analysis  is 
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beyond  the  original  intent  of  this  program,  the  approach  adopted  affords 
a  broader  range  of  generalization  than  would  otherwise  be  possible  and  is 
the  basis  for  describing  crack  growth  due  to  surface  flaw  interactions 
with  the  dynamic  stress  field  (Section  VI. 1). 

The  solution  to  the  general  problem  of  an  elastic  half-space 
subjected  to  an  arbitrary  pressure  distribution  on  its  surface  will  pro¬ 
vide  the  transient  stress  states  within  the  half-space.  The  only  mater! -.I 
properties  included  in  this  analysis  are  the  density,  Young's  modulus,  and 
Poisson's  ratio  for  the  target  material  which  enter  the  solution  through 
the  characteristic  wave  velocities  for  the  material.  The  approach  taken 
is  to  determine  the  Green's  tensor  for  Che  elastic  half-space.  The  tran¬ 
sient  stress  states  are  then  found  in  terms  of  a  surface  integral  of  the 
convolution  of  the  pressure  distribution  and  the  Green’s  tensor.  Trie  for¬ 
mal  problem  of  determining  the  Green’s  tensor  has  been  solved  and  has  been 
reduced  to  evaluating  several  definite  integrals  as  described  in  Section 
V.4.  In  some  cases  singularities  occur  in  these  integrals  which  have  been 
treated  analytically  to  preserve  the  accuracy  of  the  results.  The  singular¬ 
ities  are  understood  and  isolated  so  the  resulting  integrals  caii  be  evalu¬ 
ated  numerically  with  a  high  degree  of  accuracy.  The  solution  is  essen¬ 
tially  complete;  it  now  has  to  be  programed  for  digital  computer  evaluation. 

The  alternative  method  for  determining  the  stress  waves  is  the 
completely  numerical  finite  element  approach.  We  have  avoided  using  these 
methods  because  of  the  large  expense  and  the  fact  that  these  methods  do 
not  give  a  very  detailed  picture  of  the  stress  waves  unless  a  very  small 
grid  size  is  used.  One  is  then  beset  with  problems  resulting  from  the  dis¬ 
parity  in  the  size  of  the  regions  important  for  stress  wave  propagation 
and  the  size  of  the  crack  tip. 

Within  a  continuum  context  the  fracture  problems  to  be  addressed 
Involve  introduction  of  the  fracture  toughness  and  the  distribution  of 
critical  flaw  sizes  for  the  target  material.  For  purely  elastic  material 
response  conditions  the  fracture  toughness,  ,  is  considered  to  be  a 
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material  constant,  and  as  such  is  used  to  define  the  critical  stress 
intensity  factor  for  a  wide  variety  of  crack  configurations  and  loading 
conditions.  The  fracture  problems  identified  as  being  most  significant  in 
relation  to  the  water  drop  impact  damage  produced  in  CVD  zinc  sulfide  are 
the  stress  wave/surface  flaw  interactions  and  hydraulic  penetration  of  pre¬ 
existing  surface  cracks.  It  was  originally  thought  that  the  stress  wave/ 
microstructure  interactions  would  be  much  more  significant  than  the  experi¬ 
mental  results  indicate.  The  general  observation  from  the  experimental 
data  is  that  the  fractures  produced  during  water  drop  collisions  on  CVD 
ZnS  are  essentially  insensitive  to  the  material's  general  microstructure 
and  inclusion  content  over  a  representative  range  of  grain  size  variations. 

Stress  wave/surface  flaw  interactions  will  be  the  first  step  in 
utilizing  the  available  fracture  data  in  a  model.  This  model  will  also  be 
considered  in  terms  of  the  experimental  data  obtained  for  water  drop  impacts 
on  CVD  zinc  sulfide.  Some  siiaple  examples  of  hov»  this  data  can  provide  use¬ 
ful  estimates  of  the  fracture  response  of  CVD  ZnS  will  be  described. 

The  critical  stress  intensity  factor  (fracture  toughness)  for  CVD  ZnS 

has  been  determined  by  Evans  (data  referred  in  Evans  and  Wilshaw,  1976), 

Wiromer  and  Graves  (1977),  and  Shockey,  et  al.  (1977).  Evans  found 
]  /2 

=  1  MPa  m  for  a  material  with  an  average  grain  size  of  10  ym  with 

essentially  equiaxed  grains  (Evans,  1979).  According  to  Wimraer  and  Graves, 

1/2 

=  0.67  MPa  m  .  This  value  was  determined  for  the  baseline  material 

for  this  investigation,  material  A,  using  a  four-point  bend  test.  Wimmer 

and  Graves'  grain  size  measurements  showed  that  the  material  had  columnar 

grains  with  a  cross  sectional  diameter  of  8  to  10  ym  and  a  length  up  to 

60  ym  in  the  growth  direction.  They  also  determined  that  the  bend  strength 

1/2 

was  100.8  MPa.  Shockey,  et  al,,  found  =  0.75  MPa  m  based  on  four 

tests  using  a  modified  expanded  ring  test  technique.  The  material  used  by 
Shockey,  et  al.,  was  composed  of  columnar  grains  with  average  cross  sec¬ 
tional  diameters  of  7  ym  and  lengths  of  approximately  56  ym. 

The  grain  size  was  reported  to  be  30  ym  in  Evans  and  Wilshaw  (1976),  however 
in  a  recent  communication  with  Dr.  EvaiiS  (Evans,  1979)  the  grain  size  was 
restated  as  being  10  ym  and  equiaxed  grains.  Representatives  from  Raytheon 
Corp.  uo  not  recall  producing  equiaxed  grain  material  and  have  only  been 
able  to  achieve  a  minimum  average  aspect  ratio  of  about  3  with  the  CVD  process. 
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There  are  two  crack  geofetrtes  which  are  relevant  to  the  water  drop 
impact  conditions:  an  edge  crack  in  a  plate  and  a  teml circular  crack  in  a 
half-space  as  illustrated  in  Figure  41  .  The  stress  intensity  factors  for 
each  of  these  crses  when  the  cracks  are  subjected  to  an  uniaxial  stress, 
a  ,  are 


Kj  =  1.12av^ 

(3.1) 

[Paris  and  Sih, 

1965]  and 

K  =  C  ^ 

I  TT 

(8.2) 

[Smith,  et  al., 

6  =  90°  (refer 

1967]  where  C  =  1.03  when  6=0° 

to  Figure  41) , 

and  C  =  1.21  when 

We  note  in  passing  that  Shockey  and  co-workers  (1977)  incorrectly 
evaluated  Eq.(8.2)  using  their  value  of  .  They  equated  a  to  the 

flexual  strength  of  110  MPa  and  computed  the  magnitude  of  the  critical 

ic 

flaw  size  which  they  define  as  a  =  2  a  where  a  is  the  value  corres- 

c  c  c 

ponding  to  Eq.(8.2).  Shockey,  et  al.,  find  that  a*  =  0.73  mm  .  However 

* 

t.iey  made  an  arithmetic  error;  the  value  of  a^  should  be  0.073  mm.  Fur¬ 
thermore  if  the  critical  flaw  size  is  properly  interpreted,  then  a^  (the 
depth  of  the  flaw  as  shown  in  Figure  41  b  )  is  36  pm.  This  result  indicates 
the  critical  flaw  size  is  approximately  five  times  the  cross  sectional  diam¬ 
eter  of  the  average  grain  in  contrast  to  the  value  determined  by  Shockey 
and  co-workers  which  is  100  times  this  dimension. 

The  relations  in  Eq.(8.1)  and  (8,2)  can  also  be  used  in  conjunction 
with  the  transient  stress  s  '>tes  due  to  particulate  collisions.  Since  our 
evaluation  of  the  transient  stresses  is  not  yet  operational  for  digital 
computations,  some  of  the  published  results  will  be  used  as  approximations 
to  the  stresses  imposed  on  the  crack. 
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Adler  and  Hooker  (1976)  showed  that  Blowers'  model  and  analysis  of 
a  water  drop  impact  could  be  used  to  estimate  the  stresses  in  an  elastic 
half-space  (Blowers,  1969).  For  brittle  materials  the  critical  stress 
conditions  occur  at  and  near  the  surface  of  the  half-space  outside  of  the 
expanding  water  drop  contact  zone.  Hackworth  and  Kocher  (1977)  used  the 
computer  program  for  Blowers'  analytical  solution  to  evaluate  many  impact 
conditions  including  the  spatial  and  temporal  development  of  the  radial 
stress  component  for  low  velocity  water  drop  impacts  on  zinc  sulfide. 

Hackworth  and  Kocher  recorded  the  values  of  the  peak  radial  tensile  stresses, 
however  they  Ignored  any  criterion  for  the  onset  of  pressure  release.  Their 
subsequent  work  (Hackworth  and  Kocher,  1978)  attempted  to  incorporate  the 
pressure  profile  developed  by  Rosenblatt,  et  al .  (1977)  into  Blowers’  formu¬ 
lation  of  the  water  drop  Impact  problem.  The  deficiencies  in  their  modifi¬ 
cation  of  Blowers'  original  derivation  arc  discussed  in  Section  V.l.  The 
radial  stress  levels  they  report  are  for  relatively  long  times  after  the 
initial  contact  with  the  target.  More  representative  conditions  for  the 
near-field  stress  levels  were  provided  in  Section  III. 3b. 

Assuming  as  a  first  approximation  that  the  radial  tensile  stress  is 

responsible  for  crack  propagation,  the  magnitude  of  this  stress  will  be 

used  in  conjunction  w^th  Eq.(8.1)  and  (8.2)  to  estimate  the  potential  for 

significant  stress  wave/surface  flaw  interactions.  Representative  stress 

levels  will  be  considered  in  relation  to  the  critical  flaw  size,  a  ,  which 

c 

may  be  responsible  for  une  observed  cracks  in  water  drop  impacted  zinc  sul¬ 
fide.  Since  the  critical  stre»is  intensity  fai'.tor  is  known  and  estimates  of 
the  applied  tensile  stress  can  be  made,  a^  can  be  evaluated  using  Eqs.(8.1) 
and  (8.2).  The  parameters  associated  with  the  two-dimensional  case,  F.q.(8.1), 
will  be  denoted  by  primes  and  those  associated  with  the  semicircular  flaw  by 
double  primes,  Eq.(8.2).  When  C  =  i.03  the  fracture  stress  for  the  semi¬ 
circular  flaw  is  considerably  larger  than  that  required  for  an  edge  crack 
when  a^  =  a^"  :  the  depth  of  the  crack  is  the  same  for  both  cases. 

^IC 

ai  =  0.505  —  (8.3) 

c 
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0”  =  0.856 
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According  to  Haclc/orth  and  Kocher  (1977),  the:  radial  tensile:  stresses  for 

a  2  cm  drop  inpacting  at  222  ms  acting  long  enough  to  propagate  cracks 

reach  a  level  of  170  IfPa  with  a  peak  value  of  196  MPa  at  a  deoth  of  3 2. 7  pm 

below  the  surface.  Setting,  a^  =  12..7  ym  in  Eq.(8.3)  and  (8.4),  we  find 

0*  =  95  ?C?a  and  a”  =  .'’■,62  MPa.  The  fracture  stresses  when  a  =  5  ym  are 
F  P  c 

0p  =  151  MPa  and  0^  “  256  MPa^  The  transient  radial  tensile  stresses  for 
this  case  attain  an  effective  level  of  255  MPa  with  a  peak  value  of  285  MPa. 
The  magnitude  of  the  tensile  stresses  generated  depends  on  the  impact  veloc¬ 
ity  acid  drop  siae,  however  these  calculations  indicate  that  the  potential 
for  crack  growth  does  prevail  for  the  applied  stress  conditions  and  flaw 
sizes  which  appear  to  be  consistent  with  our  experimental  observations. 

The  comparable  teet  condition  would  be  shot  no.  815  which  is  a  1.60  mm 
water  drop  impact  at  237  m$  .  The  observed  maximum  fracture  depth  is 
40  ym  and  several  cracks  on  the  order  of  5  ym  were  seen. 


The  plots  of  the  radial  stress  component  in  Section  III.3.b  indicate 

that  appreciably  larger  stresses  are  generated  for  higher  impact  veloc- 
-1 

ities.  At  700  ms  radial  stresses  on  the  order  of  1  GPa  occur  in  the 
target  at  a  depth  of  5  ym.  These  stress  levels  have  the  potential  for 
propagating  cracks  from  critical  flaws  on  the  order  of  0,1  ym.  The  magni¬ 
tude  of  the  radial  stress  will  also  increase  as  the  surface  of  the  target 
is  approached.  Thus  the  surface  finish  has  a  strong  influence  on  the  propa¬ 
gation  of  cracks  into  the  substrate.  Concentrating  on  cracks  which  do  not 
appear  to  originate  from  surface  scratches,  it  would  seem  reasonable  to 
speculate  that  a  large  proportion  of  the  cracks  originate  from  nucleating 
flaws  less  tuan  the  grain  size.  The  correlation  of  the  radial  distance  at 
which  fractures  are  first  observed  and  the  development  of  significant  radial 
tensile  stresses  as  described  in  Section  Ill»3,b  appears  to  be  consistent  with 
the  experimental  data  over  the  j.ir.pact  velocity  range  from  200  to  700  ms 
as  summarized  in  Table  8. 
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Rossnblatt  and  co-workers  (1976,  1977)  have  also  considered  the 
problems  outlined  here  for  a  zinc  selenide  target  material  with  a  mean 
grain  size  of  50  pm.  The  smallest  grid  sizes  used  in  this  work  for  the 
finite  difference  representation  of  the  target  material  was  25  ym  for  a 
1  mm  water  drop  which  would  be  equivalent  to  50  ym  for  a  2  mm  water  drop. 

The  typical  grid  used  in  this  work  for  a  2  mm  water  drop  was  100  ym.  Sub¬ 
sequent  calculations  (Rosenblatt  and  Ito,  1978)  use  a  12.5  ym  grid  size 
for  a  2  mm  water  drop.  The  peak  radial  tensile  stresses  were  evaluated  as 
a  function  of  depth  into  the  target.  This  computation  would  only  supply 
the  values  of  the  stresses  averaged  over  a  distance  of  12.5  ym  from  the 
surface  and  applied  to  the  midpoint  of  this  cell.  The  available  analytic 
solution  indicates  there  is  a  significant  gradient  in  the  radial  stress 
with  distance  into  the  target  over  the  12.5  ym  cell  dimension,  so  a  con¬ 
siderably  finer  grid  would  be  required  to  provide  comprehensive  stress  data 
for  zinc  sulfide.  The  analytic  solution,  can  compute  the  values  of  the 
stress  components  in  increments  of  hundredths  or  thousandths  of  a  micrometer 
if  desired  with  computational  times  running  on  the  order  of  seconds  for 
each  depth  over  the  complete  range  of  radial  distances  that  these  stresses 
are  significant.  The  approach  used  by  Rosenblatt  and  Ito  for  the  surface 
crack  interaction  is  different  than  that  being  pursued  here. 

One  of  the  central  issues  pertaining  to  hydraulic  penetration  is  to 
understand  the  mechanism  by  which  hydraulic  penetration  of  cracks  leads  to 
crack  growth  and  erosion  pit  nucleation.  Obviously,  this  problem  is  com¬ 
plicated  in  general  because  of  the  combined  effects  of  the  hydrodynamics 
of  the  water  in  the  crack  and  the  interaction  with  the  locally  induced 
stress  waves  generated  during  the  penetration  process.  As  a  first  step 
toward  describing  this  process,  a  much  simpler  problem  has  been  considered: 
namely,  a  two-dimensional  pressurized  surface  crack  inc.lined  to  the  surface. 
The  result  of  this  calculation  is  the  stress  distribution  surrounding  the 
crack  with  particular  emphasis  on  the  effect  produced  by  the  proximity  of 
the  surface.  From  this  stress  distribution,  the  initicl  direct?on  and  pos¬ 
sibly  the  velocity  cf  crack  propagation  as  a  function  of  initial  crack 
angle  and  depth  will  be  determined  using  currently  available  frecture  models. 
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It  is  our  contention  that  the  fracture  problems  identified  and 
presently  being  analyzed  in  terms  of  a  continuum  mechanics  approach  are  the 
critical  problems  which  must  be  addressed  in  order  to  formulate  a  productive 
analytical  framework  for  describing  the  particle  impact  response  of  CVD  zinc 
sulfide.  These  analyses  have  been  carried  out  in  a  balanced  association 
with  observations  of  the  actual  effects  of  the  water  drop  collision  process. 
Experimental  procedures  have  been  utilized  to  gain  understanding  of  the  re¬ 
lationships  between  the  material’s  surface  condition  and  micros tructure  and 
the  nature  and  extent  of  the  resulting  crack  patterns. 

The  experimental  observations  of  the  fractures  produced  by  water  drop 
impingement  indicate  a  relative  insensitivity  to  the  grain  size  and  general 
grain  structure  in  CVD  ZnS.  On  the  other  hand,  Hackworth  and  Kocher  (1978) 
report  seeing  a  significant  reduction  in  damage  coincident  with  reduction 
of  grain  size  in  ZnSe.  This  result,  seemingly  in  conflict  with  the  present 
finding,  is  easily  reconciled  by  the  following  observations.  The  grain  size 
used  in  Hackworth  and  Kocher's  (1978)  work,  although  not  reported  numerically, 
appears  in  their  micrograph  to  be  a  large  grain  size  possibly  on  the  order 
of  100  ym  intercept  distance  in  the  plane  shown  and  unresolvable  at  the 
stated  magnification  of  135  (presiunably  <10  ym)  for  the  micrograph  of  the 
fine  grained  material.  This  grain  size  variation  represents  not  only  a 
much  larger  range  than  this  study,  but  the  larger  grain  size  is  cn  the 
order  of  the  impact  fracture  pattern  dimensions.  Combining  this  observation 
with  the  observations  made  on  ZnS  would  indicate  that  reductions  in  grain 
size  beyond  an  order  of  magnitude  below  the  fracture  size  will  not  yield  a 
signific  c.nt  erosion  resistance. 

Examination  of  the  micrographs  of  water  drop  impacts  presented  by 
Hackworth  and  Kocher  (1978)  to  illustrate  that  Che  fine  grain  ZnSe  has 
better  erosion  resistance  than  the  large  grain  ZnSe  reveals  that  although 
the  fractures  are  smaller  in  the  small  grain  material  there  are  also  many 
more  of  them.  This  indicates  that  most  probably  a  similar  amcunt  of  frac¬ 
ture  surface  area  was  created  for  each  impact;  i.e.,  a  similar  amount  of 
energy  dissipated  in  fracturing  the  specimen  but  the  form  of  the  fractures 
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changed.  Hackworth  and  Kocher  (1978)  also  report  that  although  the  resist¬ 
ance  of  ZnSe  improved  with  reduced  grain  size  it  did  not  equal  that  of  ZnS. 
This  indicates  that  erosion  resistance  is  governed  by  a  more  intrinsic 
material  variable  than  grain  size  or  hardness. 

The  effect  of  grain  size  and  microhardness  on  the  erosion  resistance 
of  zinc  sulfide  has  been  described  by  Adler  and  Hooker  (1978).  It  is  pointed 
out  that  substantial  strengthening  would  be  anticipated  when  the  grain  bound¬ 
aries  effectively  resist  fracture  propagation,  when  fracture  propagation  is 
confined  to  crystallographic  cleavage  planes,  and  when  the  grain  size  is 
much  smaller  than  the  drop  size.  With  the  exception  of  confirmation  of 
fracture  propagation  along  crystallographic  cleavage  planes  these  condi¬ 
tions  are  satisfied  by  the  CVD  ZnS  formulations  which  have  been  investi¬ 
gated.  Furthermore  Adler  and  Hooker  indicated  that  the  standard  mechanical 
fracture  tests  evaluate  the  response  to  the  largest  flaw  present,  however 
the  erosion  process  depends  upon  the  concurrent  propagation  of  a  surface 
flaw  population.  The  largest  flaw  may  not  be  representative  of  the  average 
behavior  of  the  population.  It  is  also  shovm  here  that  the  magnitude  of 
the  imposed  radial  stresses  for  2  mm  water  drops  impacting  at  velocities 
in  excess  of  300  ms  ^  is  many  times  the  calculated  fracture  stress  levels 
for  reasonable  flaw  dimensions  as  identified  from  the  I'VD  ZnS  specimens 
evaluatedc  In  view  of  the  overwhelming  amount  of  energy  available  and  the 
imposed  stress  levels,  it  is  unlikely  that  further  modification  of  the  bulk 
material  will  significantly  improve  the  erosion  resistance  of  CVD  ZnS. 

Considerably  greater  latitude  for  increased  resistance  to  rain 
erosion  can  be  gained  from  improved  polishing  procedures  and  subsequent 
maintenance  of  the  initial  surface  finish.  The  specimens  used  in  this 
program  were  polished  to  better  than  average  FLIR  window  specifications, 
but  the  capability  of  the  water  drop  to  interact  with  the  polishing  defects 
remaining  even  for  this  level  of  surface  preparation  is  clearly  evident  in 
the  water  drop  experiments.  Submicrometer  flaw  depths  will  have  to  be 
achieved  and  maintained.  The  application  of  the  antiref lectant  coating  over 
this  surface  quality  viill  certainly  have  a  positive  influence  on  improving 
the  erosion  resistance. 
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SECTION  IX 
CONCLUSIONS 

The  general  fracture  behavior  of  CVD  ZnS  due  to  water  drop  impingement 
has  been  investigated.  Three  zinc  sulfide  formulations  have  been  considered 
representing  a  range  of  grain  sizes,  grain  shapes,  and  inclusion  content. 

The  microscopic  observations  of  the  damage  produced  by  2  mm  diameter  water 
drops  impacting  from  222  to  684  ms  ^  indicate  that  there  is  no  discernable 
difference  in  the  fracture  trajectories  or  penetration  depths  for  the  micro- 
structural  variations  examined.  It  is  found  tlat  no  significant  improvement 
in  the  erosion  resistance  will  accrue  from  further  modification  of  the  bulk 
properties  of  the  CVD  processed  ZnS.  The  dominant  fracture  mode  is  trans- 
granular,  and  the  grain  size  is  small  enough  in  relation  to  the  loaded 
region  of  a  1  aim  or  larger  water  drop  that  no  benefit  will  be  gained  from 
further  grain  refinement. 

The  most  important  factor  in  relation  to  the  fracture  response  of  the 
currently  produced  CVD  ZnS  was  the  quality  of  the  surface  finish.  Calcula¬ 
tions  were  provided  which  indicated  the  imposed  tensile  stress  levels  for 
2  mm  water  drops  impacting  at  200  ms  were  comparable  to  the  fracture 

stress  for  reasonable  surface  scratch  dimensions  to  several  times  the  esti- 

-i 

mated  fracture  stress  as  the  impact  velocity  approached  700  ms  .  On  the 
basis  of  these  straightforward  estimates  of  the  critical  stress  levels,  a 
significant  increase  in  the  erosion  resistance  appears  to  be  possible  if 
the  surface  polishing  defects  do  not  exceed  a  depth  of  1  ym  and  preferably 
0.1  ym. 


The  need  for  evaluating  the  transient  stress  states  at  and  near  the 
surface  of  the  specimen  was  an  outgrowth  of  the  experimental  program  where 
it  appears  that  the  fracture  nucleating  surface  flaws  are  on  the  order  of 
5  ym  or  less.  The  available  analytical  and  numerical  stress  wave  analyses 
were  to  be  used  for  modeling  the  observed  fracture  response  of  CVD  ZnS,  how¬ 
ever  these  analyses  were  found  to  be  inadequate  for  various  reasons  which 
have  been  enumerated.  This  survey  of  available  approaches  motivated  the 
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development  of  the  mathematical  analysis  of  the  transient  stresses  due  to 
a  water  drop  collision  presented  in  Section  V.  Further  work  is  required 
to  obtain  numerical  results,  however  a  generalization  of  Blowers'  analyti¬ 
cal  approach  is  now  available  which  should  be  applicable  to  a  spatially 
and  temporally  varying  loading  function  applied  to  the  surface  of  the 
target . 

Two  fracture  conditions  have  been  identified  as  being  the  most 
relevant  to  the  general  degradation  of  infrared  transmitting  windows. 

The  first  is  the  interaction  of  the  transient  stresses  generated  by  the 
impact  with  surface  flaws,  and  the  second  is  hydraulic  penetration  of  the 
water  from  an  impacting  water  drop  into  a  pre-existing  surface  crack. 

Stress  analyses  for  both  of  these  critical  problems  have  been  initiated  as 
described  in  Section  VI. 

The  transmission  degradation  which  may  arise  due  to  cracks  in  the 
material  once  they  are  formed  also  has  been  treated  analytically.  The 
significance  of  the  crack  face  separation  distance  and  the  crack  orienta¬ 
tion  is  estimated  in  the  preliminary  calculations  presented  in  Section  VII, 
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